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2.  “Gamma-Ray  Bursts  and  the  Particle  Mass  Scale,”  Yang  Dai,  Alex  B.  Borisov,  James  W. 
Longworth,  Keith  Boyer,— and  Charles  K.  Rhodes.  Proceedings  of  the  International 
Conference  on  Electromagnetics  in  Advanced  Applications,  edited  by  R.  Graglia  (Politecnico  di 
Torino,  Torino,  Italy,  1999)  p.  3. 
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Strong-Electroweak  Coupling  Constant,  and  the  Higgs  Mass,”  Y.  Dai,  A.  B.  Borisov,  K.  Boyer, 
and  C.  K.  Rhodes,  Sandia  National  Laboratories,  Report  SAND2000-2043,  August,  2000. 

4.  “An  Efficient,  Selective  Collisional  Ejection  Mechanism  for  Inner-Shell  Population  Inversion 
in  Laser  Driven  Plasmas,”  W.  A.  Schroeder,  T.  R.  Nelson,  A  B.  Borisov,  J.  W.  Longworth,  K. 
Boyer,  and  C.  K.  Rhodes,  J.  Phys.  B,  in  press. 

5.  “L-shell  Emission  from  High  Z-solid  Targets  by  Selective  Inner-shell  Collisional  Electron 
Ejection,”  T.  R.  Nelson,  A.  B.  Borisov,  K.  Boyer,  S.  Cameron,  J.  W.  Longworth,  T.  S.  Luk,  A. 
McPherson,  W.  A  Schroeder,  J.  Santoro,  A  Van  Tassle,  and  C.  K.  Rhodes  submitted  to  J. 
Phys.  B,  December,  2000. 

(4)  Statement  of  Problem  Studied 

Three  computationally  intensive  problem  area  were  studied.  They  are  (A)  power  compression  in 
materials,  particularly  self-trapped  plasma  channels  [1],  (B)  the  phenomenon  of  y-ray  bursts  [2]  and  the 
relationship  to  the  particle  mass  spectrum  [3],  and  (C)  the  potential  role  of  a  new  organizing  principle  in 
the  regulation  of  biological  processes. 


(5)  Summary  of  the  Most  Important  Results 

A.  Power  Compression  (See  Appendix  1). 

Stable  self-channeling  of  ultra-powerful  (Po  ~  1  TW  -  1  PW)  laser  pulses  in  dense  plasmas  is  a  key 
process  for  many  applications  requiring  the  controlled  compression  of  power  at  high  levels.  Theoretical 
computations  predict  that  the  transition  zone  between  the  stable  and  highly  unstable  regimes  of 
relativistic/charge-displacement  self-channeling  is  well  characterized  by  a  form  of  weakly  unstable 
behavior  that  involves  bifurcation  of  the  propagating  energy  into  two  powerful  channels.  Recent 
observations  of  channel  instability  with  femtosecond  248  nm  pulses  reveal  a  mode  of  bifurcation  that 
corresponds  well  to  these  theoretical  predictions.  It  is  further  experimentally  shown  that  the  use  of  a 
suitable  longitudinal  gradient  in  the  plasma  density  can  eliminate  this  unstable  behavior  and  restore  the 
efficient  formation  of  stable  channels. 
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B.  Gamma-Ray  Bursts  and  Elementaiy  Particle  States 

1.  (ramma-Rav  Burst  Analysis  (See  Appendix  21 

A  mode  of  detonation  for  neutron  stars  is  derived  based  on  a  physically  defined  power 
residue  counting  system  used  to  construct  a  particle  mass  scale  that  reflects  supersymmetric 
fermion/boson  pairing.  Neutron  stars  with  specific  mass  values  (Mn)  can  become  unstable 
toward  decay  into  an  ensemble  of  highly  energetic  pairs  of  supersymmetric  prides  with  an 
energy  release  whose  upper  bound  can  approach  the  rest  mass  energy  of  the  original  body.  The 
subsequent  scattering  of  these  particles  provides  a  way  to  account  for  (1)  the  high  radiative 
yield  J)  of  gamma-ray  bursts,  (2)  the  large  variability  of  their  observed  energy,  and 

(3)  the  expected  collateral  production  of  ultrahigh  energy  cosmic  rays.  The  joint  action  of  (a) 
supersymmetric  classification  and  (b)  a  new  arithmetic  selection  rule  designated  as  congruence 
conservation  governs  the  condition  of  instability  and  the  dynamic  participation  of  the  magnetic 
monopole  obviates  baryon  number  conservation.  The  decay  patterns  of  certain  elernentary 
particles  and  putative  supermassive  forms  of  cosmological  dark  matter  may  be  subject  to 
similar  conditions,  -l^cifically,  the  analysis  predicts  the  existence  of  a  super-light  boson 
neutrino  Vi  that  plays  a  fundamental  role  in  weak  interactions  and  it  is  conjectured  that  the 
phenomenon  of  lepton  number  conservation  is  isomorphic  to  the  supersymmetric  conservation 
of  congruence.  The  vi  boson  also  provides  (i)  a  potential  basis  for  the  CP  violation  observed  in 
K£  decay  and  (ii)  an  alternative  description  of  the  attribute  of  strangeness.  These 
interpretations  support  the  conclusion  that  the  observations  of  gamma-ray  bursts  signal  the 
existence  of  a  radically  new  organizing  principle  of  physical  interactions  that  defines  a  general 
condition  of  stability  for  massive  particles  and  expresses  the  unstable  case  with  the  copious 
production  of  energetic  emissions. 

2.  Flementarv  Particle  States  (See  Appendix  3). 

The  construction  of  inverse  states  in  a  finite  field  F  p^enables  the  organiption  of  the  mass 
scale  with  fimdamental  octets  in  an  eight-dimensional  index  space  that  identifies  particle  states 
with  residue  class  designations.  Conformance  with  both  CPT  invariance  and  the  concept  of 
supersymmetry  follows  as  a  direct  consequence  of  this  formulation.  Based  on  two  parameters 
(P  and  gg)  that  are  anchored  on  a  concordance  of  physical  data,  this  treatment  leads  to  (1)  a 
prospective  mass  for  the  muon  neutrino  of  ~  27.68  meV,  (2)  a  value  of  the  unified  strong- 
electroweak  coupling  constant  a*  =  (34.26)"^  that  is  physically  defined  by  the  ratio  of  the 
electron  neutrino  and  muon  neutrino  masses,  and  (3)  a  seesaw  congruence  connecting  the 
Higgs,  the  electron  neutrino,  and  the  muon  neutrino  masses.  Specific  evaluation  of  the  masses 
of  the  corresponding  supersymmetric  Higgs  pair  reveals  that  both  particles  are  superheavy 
(>10^*  GeV).  No  renormalization  of  the  Higgs  masses  is  introduced,  since  the  calculational 
procedure  yielding  their  magnitudes  is  intrinsically  divergence-fi’ee.  Further,  the  Higgs  fulfills 
its  conjectured  role  through  the  seesaw  relation  as  the  particle  defining  the  origin  of  all  particle 
masses,  since  the  electron  and  muon  neutrino  systems,  together  with  their  supersymmetric 
partners,  are  the  generators  of  the  mass  scale  and  establish  the  corresponding  index  space. 
Finally,  since  the  computation  of  the  Higgs  masses  is  entirely  determined  by  the  modulus  of  the 
field  P^,  which  is  fully  defined  by  the  large-scale  parameters  of  the  universe  through  the  value 
of  the  universal  gravitational  constant  G  and  the  requirement  for  perfect  flatness  (Q  =1.0),  the 
seesaw  congruence  fuses  the  concepts  of  mass  and  space  and  creates  a  new  unified  archetype. 
On  the  basis  of  new  data  firom  the  Super-Kamiokande  Collaboration  [4],  this  analysis  has  been 
extended  to  give  prospective  values  for  the  electron,  muon,  and  tau  neutrinos  as  given  below. 
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The  previously  published  atmospheric  neutrino  data  did  not  distinguish  whether  muon  neutrinos  were 
oscillating  into  tau  neutrinos  or  sterile  neutrinos,  as  both  hypotheses  fit  the  data.  Using  data  recorded  in 
1 100  live  days  of  the  Super-Kamiokande  detector,  we  use  three  complementary  data  samples  to  study 
the  difference  in  zenith  angle  distribution  due  to  neutral  currents  and  matter  effects.  We  find  no  evidence 
favoring  sterile  neutrinos,  and  reject  the  hypothesis  at  the  99%  confidence  level.  On  the  other  hand,  we 
find  that  oscillation  between  muon  and  tau  neutrinos  suffices  to  explain  all  the  results  in  hand. 

PACS  numbers:  14.60.Pq,  %.40.Tv 

5m^  = 

6m^  =  3.2x  10”^  (eV^) 

(MEASUREMENT) 

mj,  =  27.68  X  10"^  eV 
=  62.98  meV 


Vx 

62.98  meV 

V, 

27.68  meV 

Ve 

0.808  meV 
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C.  Biological  Processes  • .  l-  i  • 

A  genetic  analogy  linking  the  physical  mass  states  of  elementary  systems  with  biological 

process  has  been  developed.  Specifically,  the  subgroup  structure  of  the  mass  ststes  [5]  suggests  a 
direct  link  with  the  large-scale  organization  of  metabolic  networks  [6].  A  connection  to  the 
dynamics  of  the  p53  tumor-suppressor  gene  [7]  is  also  seen.  A  statement  [5]  concerning  the  genetic 

analogy  is  given  below.  .  l  •  r  /in  u 

Arithmetic  conditions  between  particle  masses  can  be  defined  on  the  basis  of  (1)  the 
supersymmetric  conservation  of  congruence  and  (2)  the  observed  characteristics  of  particle 
reactions  and  stability.  Stated  in  the  form  of  common  divisors,  these  relations  can  be  interpreted  as 
expressions  of  genetic  elements  that  represent  elementary  particle  characteristics.  In  order  to 
illustrate  this  concept,  it  is  shown  that  the  pion  triplet  (tc*,  tc")  can  be  identified  by  the  existence  of  a 
common  divisor  d<rf:  in  a  way  that  can  also  account  for  the  observed  7t*/7i“  mass  splitting.  The 
results  indicate  that  a  corresponding  statement  may  also  hold  for  all  particle  multiplets. 
Classification  of  the  respective  physical  states  is  achieved  by  association  of  the  common  divisors 
with  residue  class  designations  in  a  finite  field  Fp^and  the  existence  of  the  multiplicative  group  of 
units  F*  enables  the  coircsponding  mass  parameters  to  be  associated  with  a  rich  subgroup 
structure!  It  is  further  shown  that  the  fundamental  source  of  supersymmetry  can  be  expressed  in 
terms  of  hierarchical  relationships  between  odd  and  even  order  subgroups  of  Fp*,  an  outcome  that 
reflects  itself  in  the  phenomenon  of  fermion^oson  pairing  of  individual  particle  systems. 
Accordingly,  supersymmetry  is  a  group  rather  than  a  particle  property.  The  status  of  the  Higgs 
subgroup  of  order  4  is  singular;  it  is  isolated  from  the  hierarchical  pattern  and  communicates 
globally  to  the  mass  scale  through  the  seesaw  congruence  by  specifying  the  generators  of  the 
physical  masses. 
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Abstract 

Stable  self-channeling  of  ultra-powerful  (Po  ~  1  TW  -  1  PW)  laser  pulses  in  dense  plasmas  is  a 
key  process  for  many  applications  requiring  the  controlled  compression  of  power  at  high  levels 
Theoretical  computations  predict  that  the  transition  zone  between  the  stable  and  highly  unstable  regimes 
of  relativistic/charge-displacement  self-channeling  is  well  characterized  by  a  form  of  weakly  unstable 
behavior  that  involves  bifurcation  of  the  propagating  energy  into  two  powerful  channels.  Recent 
observations  of  channel  instability  with  femtosecond  248  nm  pulses  reveal  a  mode  of  bifurcation  that 
corresponds  well  to  these  theoretical  predictions.  It  is  further  experimentally  shown  that  the  use  of  a 
suitable  longitudinal  gradient  in  the  plasma  density  can  eliminate  this  unstable  behavior  and  restore  the 
efficient  formation  of  stable  channels. 


I.  Introduction 

Robust  stabiUty  is  known  both  theoretically  [1-3]  and  experimentally  [1-4]  to  be  a  chief 
characteristic  of  relativistic/charge-displacement  self-channeling  [1-12].  Moreover,  it  is  known  that  the 
influence  of  the  ponderomotively  generated  charge-displacement  is  the  major  factor  [2]  governing  the 
stability  of  the  confined  modes  of  propagation  and  it  has  been  observed  that  this  feature  has  a 
counterpart  in  the  ponderomotive  stabilization  present  in  related  areas  of  plasma  physics  and  accelerator 
design  [2],  It  is  also  theoretically  clear  [2,3]  that  the  eigenmodes  play  a  key  role  in  the  dynamics  of  the 
stability,  since  they  function  as  attractors  and  efficiently  guide  the  energy  into  stable  states  of  high 
power  compression. 

A  regime  of  unstable  propagation  is  also  theoretically  predicted  [2]  to  exist,  although  a 
systematic  experimental  study  of  its  properties  is  lacking.  The  unstable  behavior  can  be  roughly 
classified  into  two  subcategories.  They  are  (1)  strongly  unstable  propagation  that  results  in  catastrophic 
filamentation,  an  outcome  associated  with  an  operating  point  located  deep  in  the  unstable  region  [2],  and 
(2)  weak  instability  associated  with  the  stable/unstable  boundary  zone  in  which  deviations  from  stable 
single  channel  propagation  initially  become  apparent.  The  former  regime  is  characterized  by  the 
development  of  a  group  of  peripheral  smaU-scale  filaments  each  containing  a  power  on  the  order  of  the 
critical  power  for  self-channeling.  Under  these  circumstances,  an  increase  in  the  incident  laser  power 
leads  primarily  to  a  disorganized  multiplication  of  filaments  with  the  outcome  that  the  trapped  power  in 
an  individual  filament  does  not  increase.  In  the  case  of  weak  instability,  however,  when  typically  two 
channels  are  formed,  the  power  of  the  strongest  channel  can  be  raised  with  an  increase  of  the  incident 
laser  power.  Overall,  the  principal  characteristics  of  the  weakly  unstable  behavior  are  (a)  bifurcation  of 
the  propagating  mode  and  (b)  the  development  of  lateral  displacements  of  the  confined  power. 

The  present  study  discusses  these  irregularities  in  the  propagation  by  examining  both  their 
experimental  observation  and  the  corresponding  theoretical  description.  The  main  findings  reported  in 
this  work  are  (i)  the  first  experimental  observations  of  bifurcation  of  the  propagating  channel  and  (ii) 
experimental  confirmation  of  the  previously  proposed  [1]  use  of  a  suitable  longitudinal  gradient  in  the 
plasma  density  for  both  the  restoration  of  stable  propagation  and  channel  optimization. 


II.  Discussion  of  Results 

Experimental  evidence  for  bifurcation  is  shown  in  Fig.  1(a),  an  image  that  displays  a  transverse 
view  of  the  spectrally  integrated  Xe(M)  radiation  (~lkeV)  produced  by  the  propagating  ultraviolet 
(248nm)  radiation  in  a  gaseous  xenon  cluster  target  [4,5]  having  a  Xe  density  of  ~  3x10  cm  .  The 
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methods  of  measurement  and  the  experimental  conditions  generally  pertaining  to  the  data  in  Fig.  1(a) 
have  been  previously  described  [4],  An  interpretation  of  the  observed  Xe(M)  spectrum  has  also  been 
given  [13], 

The  key  feature  shown  in  Fig.  1(a)  is  the  dual  structure  visible  at  the  longitudinal  position 
z=  1500  |im.  Since  x-rays  are  generated,  this  indicates  the  existence  of  two  separate  parallel  zones  in 
which  the  power  is  highly  compressed.  In  contrast,  at  the  position  corresponding  to  zs600  pm,  only  a 
single  feature  is  evident.  In  the  process  of  channel  formation,  the  incident  propagating  energy  collapses 
from  this  initial  region  to  a  thin  channel  whose  volume  is  sufficiraily  small  that  the  x-ray  emission  from 
it  is  not  visible  in  the  exposure  on  account  of  the  detection  limit  imposed  by  combined  influence  of 
(i)  the  relatively  low  spatial  resolution  (~  30  pm),  (ii)  the  bound  on  the  dynamic  range  of  the  detector, 
and  (iii)  the  quenching  of  the  Xe(M)  emission  in  the  high  intensifr  zone  of  the  channel  that  is  caused  by 
stripping  of  the  entire  3d-shell  by  field  ionization.  This  behavior  causes  the  longitudinal  gap  (800  pm  < 
z  <  1400  pm)  in  the  x-ray  emission  to  occur.  The  subsequent  bifiircation  of  the  thin  (-2-3  pm  radius) 
channel  [1-3]  causes  the  dual  feature  at  z=  1500  pm  to  appear.  We  note  that  other  studies  [14]  have 
examined  the  spectral  properties  of  the  Xe(L)  emission  from  the  intense  central  region  of  the  filament. 

The  observed  morphology  of  the  x-ray  signal  can  be  put  in  contact  with  computations  of  the 
dynamics  for  conditions  corresponding  to  those  of  the  experiment.  The  modeling  of  the 
relativistic/charge-displacement  self-channeling  is  based  on  a  relatively  simple  model  [5,6]  that  involves 
two  dimensionless  parameters;  one  represents  the  normalized  power  q,  while  the  other  describes  the 
normalized  radius  po  of  the  incident  radiation.  They  are  given  respectively  by 

q  =  Pq  /  P„  and  Po  =  rOp  o  /c .  (1) 

In  Eq.  (1)  Po  represents  the  incident  peak  power  associated  with  the  central  zone  of  the  pulse  and  Per 
denotes  the  critical  power  for  relativistic/charge-displacement  self-channeling  given  by  [5-7] 

Pcr=  (me.ocVe^)  J  g^(p)pdp  ((n/(0p,o)^  =1.6198x10'°  (©  /Op,,)^  W ,  (2) 

0 

in  which  me,o,  c,  and  e  have  their  standard  identifications,  go(p)  is  the  Townes  mode  [15],  and  o, 
Qpo  =  (47ie‘ and  ro,  respectively  denote  the  angular  frequency  corresponding  to  the 

propagating  radiation,  the  angular  frequency  of  the  unperturbed  plasma,  and  the  radius  of  the  incident 
transverse  intensity  profile.  A  chief  feature  of  the  self-channeling  process  is  the  stabilization  of  the 
transverse  beam  profile  near  one  of  the  z-independent  modes  of  propagation  designated  as  the  lowest 
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eigenmodes  [5,6]  Us,o(p)  (with  index  s,  0  <  s  <  1)  of  the  governing  nonlinear  Schrodinger  equation. 
These  eigenmodes  have  the  dimensionless  radius 

p.,.-[2ju;,.(p)dp/u;,,(o)r'=.  (3) 

0 

The  self-channeling  process  can  be  described  as  the  dynamic  adjustment  of  the  radius  ro  of  the 
incident  focused  beam  to  the  radius  rd,  of  the  laser  channel,  the  value  of  which  can  be  derived  from  the 
dimensionless  radius  pe.o  of  the  corresponding  lowest  eigenmode  [3,6]  as 


Teh  = 


(4) 


Previous  analysis  [6]  has  shown  that  in  the  case  of  relativistic  and  charge-displacement  self-channeling, 
the  dimensionless  radius  of  the  channel  pe,o  varies  extremely  slowly  over  a  large  range  of  normalized 
laser  power  ri,  pe,o  «  const  for  ri  >  1.5,  a  property  evident  in  Fig.  2.  Therefore,  from  the  expression  for 
the  plasma  frequency 

cOp=(37ie'N./m..o)'^%  (5) 

we  conclude  that  the  radius  of  the  laser  channel  formed  in  a  plasma  with  a  local  electron  density  N* 
scales  in  good  approximation  [3]  simply  as 


The  situation  involving  q  =  Po  /  Per  »  1  ^md  ro  /  rch  »  1  leads  to  a  strong  self-focusing  response 
which  can  be  either  stable  or  unstable  [1,2],  Unfortunately,  even  in  the  stable  case,  the  propagation  is 
generally  characterizing  by  a  substantial  loss  of  laser  power  and  the  eflBciency  of  the  process  is  grossly 
compromised  [3].  However,  the  efficiency  of  the  power  compression  can  be  optimized  when  ro  ~  rch  (or, 
equivalently,  po  ~  pe,o),  a  condition  that  can  be  dynamically  achieved  by  the  adjustment  of  the  local 
electron  density  Ne  (see  Eq.  (6)).  This  can  be  readily  achieved  with  the  use  of  an  appropriate 
longitudinal  gradient  in  the  electron  density  distribution  [3], 

Previous  work  [2]  has  shown  that  the  conditions  governing  the  stability  of 
relativistic/ponderomotive  self-channeling  are  well  described  by  two-dimensional  (q,po)  stability  maps 
of  the  form  illustrated  in  Fig.  (2)  and  that  specific  zones  corresponding  to  stable  and  strongly  unstable 
propagation  can  be  identified.  In  particular,  the  stability  and  efficiency  of  the  relativistic/charge- 
displacement  self-channeling  mechanism  have  been  analyzed  [1,2]  for  perturbed  Gaussian  beams 
having  incident  transverse  amplitude  distributions  given  by 

Uo(r,({))  =  C  exp(-0.5(r/ro)'  x(l-h(r/ro/2e,  cos(^(t))).  (7) 
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(8) 


In  the  profile  represented  by  Eq.  (7),  the  weak  azimuthal  perturbations  are  characterized  by 
g^=0.03,  q  =  l-4,  5  =  max|Io(r,({),)-Io(r,(i)2)|/Io  =0-21- 

The  information  presented  in  the  stability  (Ti,po)  map  illustrated  in  Fig.  (2)  can  be  used  to 
establish  a  direct  correspondence  between  the  experimental  data  shown  in  Fig.  1(a)  and  the  theoretical 
procedures  applicable  to  the  dynamics  of  pulses  with  amplitude  distributions  given  generally  by  the 
form  of  Eq.  (7).  The  gaseous  xenon  cluster  target  is  experimentally  [4,14]  produced  by  a  pulsed  valve 
having  a  ~  1.5  mm  orifice  that  is  equipped  with  a  vertical  barrier  on  the  side  facing  the  incident  248nm 
beam.  This  wall,  which  contains  an  aperture  having  a  diameter  of -100  pm,  blocks  the  lateral  expansion 
of  the  gas  except  in  the  local  zone  near  the  -100  pm  opening.  Accordingly,  a  focal  position  located 
sufficiently  inside  the  wall  corresponds  well  to  a  uniform  target  density.  In  contrast,  location  of  the 
focus  in  the  region  near  or  slightly  outside  of  the  100  pm  aperture  causes  the  focussed  pulse  to 
propagate  through  a  zone  in  which  a  substantial  longitudinal  gradient  in  target  density  is  present. 
Therefore,  adjustment  of  the  focal  position  along  the  direction  of  propagation  in  the  region  defined  by 
the  barrier  enables  the  alteration  of  the  density  gradient  experienced  by  the  focussed  radiation. 
Consequently,  there  exists  a  controllable  experimental  procedure  that  can  provide  either  a  uniform  target 

density  or  one  possessing  an  appreciable  density  gradient. 

The  ability  to  modify  the  longitudinal  density  profile  enables  the  dynamic  trajectory  of  the 
operating  point  [1,3]  of  the  system  to  be  grossly  altered.  The  principal  consequence  of  this  control,  for 
the  data  shown  in  Fig.  1(a),  can  be  understood  by  reference  to  Fig.  1(b)  and  Fig.  (2).  The  initial  condition 
corresponding  to  the  development  of  the  experimentally  observed  bifurcated  channel  illustrated  in 
Fig.  1(a)  is  designated  by  the  datum  Aumfom  in  Fig.(2),  a  point  that  lies  near  the  zone  in  which  strong 
filamentation  occurs  [2].  Although  in  the  stable  region,  it  is  significant  to  note  that  Aimfom,  is  clearly 
located  at  a  considerable  remove  fi'om  the  eigeninode  curve  of  the  system  pe,  ©(ty)  in  this  situation. 

The  computed  evolution  of  an  incident  pulse  corresponding  to  Auniform  is  illustrated  in  Fig.  1(b). 
The  salient  feature  of  the  theoretical  distribution  is  a  strong  spatial  bifurcation,  a  characteristic  that  is  in 
good  qualitative  correspondence  with  the  experimental  observation  shown  in  Fig.  1(a).  Indeed, 
bifurcation  into  two  roughly  equal  channels  is  shown.  Experimentally,  the  bifurcated  propagation 
becomes  visible  when  the  divergence  of  the  two  channels  is  sufficiently  great  that  they  can  be  separately 
resolved  by  the  x-ray  camera  [4].  A  continuation  of  the  computed  evolution  of  the  pulse  to  longitudinal 
distances  Z  >  400  pm  reveals  that  the  stronger  component  of  the  bifurcated  distribution  shown  in 
Fig.  1(b)  forms  a  stable  channel  that  is  denoted  by  point  Buniform  in  Fig.(2).  We  note  that  this  datum  falls 


5 


very  close  to  the  eigenmode  curve  [pe,o  (ri)],  a  basic  property  of  fully  evolved  stable  channels  [1,3]  and 

a  sharp  contrast  to  the  initial  condition  given  by  the  location  of  Auniform- 

Previous  studies  have  demonstrated  that  stable  channels  aggressively  seek  the  eigenmode 
distribution  [6],  Hence,  the  existence  of  a  channel  with  an  operating  point  corresponding  to  a  position 
far  from  the  eigenmode  curve  necessitates  a  major  reorganization  of  its  structure.  The  present  results 
indicate  that,  for  circumstances  requiring  a  sufficiently  large  adjustment,  the  process  of  bifurcation 
enables  the  operating  point  of  the  system  to  make  a  single  large  and  abrupt  step  toward  the  eigenmode 
curve.  This  mechanism,  however,  leads  to  very  high  losses,  as  the  relative  positions  of  Auniform  and 

Buniformin  Fig.  (2)  illustrate. 

The  use  of  a  plasma  gradient  greatly  alters  the  dynamics  of  propagation  [3].  Specifically,  it  was 
theoretically  predicated  [3]  that  the  process  of  bifurcation  could  be  eliminated  by  the  use  of  a  suitable 
plasma  gradient.  This  predication  has  now  been  experimentally  confirmed,  as  shown  by  the  data 
presented  in  Fig.  3(a).  The  morphology  of  the  channel  seen  in  Fig.  1(a)  is  transformed  into  that  shown  in 
Fig.  3(a)  by  a  relocation  of  the  incident  focal  zone  to  a  region  possessing  an  appreciable  density 
gradient.  It  is  apparent  that  the  bifurcated  region  (z  =  1500  pm)  of  the  former  has  been  replaced  by  an 
intense  single  channel  in  the  latter.  This  was  achieved  experimentally  solely  by  a  small  longitudinal  shift 
in  the  focal  position  of  the  incident  pulse  to  the  region  near  the  100  pm  aperture  where  a  significant 
gradient  in  the  target  density  exists.  This  adjustment  places  the  operating  point  of  the  incident  radiation 
in  the  proximity  of  the  eigenmode. 

A  corresponding  calculation  of  the  pulse  evolution  which  includes  a  suitable  density  gradient, 
shown  in  Fig.  3(b),  illustrates  the  smooth  development  of  a  single  channel.  In  this  case,  the  imtial 
condition  has  been  remapped  to  point  Agradientin  Fig.  (2)  and  the  trajectory  of  the  operating  point  leads  to 
the  stable  configuration  by  Bgradient  Since  the  loss  caused  by  the  formation  of  the  second  subsidiary 
channel  in  the  bifurcated  case  has  been  eliminated,  the  efficiency  of  the  transport  into  the  channel 
formed  has  been  significantly  increased.  The  existence  of  the  plasma  gradient  enables  the  pulse  to 
reorganize  in  a  less  abrupt  manner,  thereby  eliminating  the  need  for  the  bifurcation. 

m.  Conclusions 

Experimental  observations  of  channel  instability  show  good  correspondence  to  computations  of 
the  channel  dynamics.  The  leading  instability  is  seen  to  be  bifurcation  of  the  propagating  mode.  The  use 
of  an  appropriate  density  gradient  can  control  this  instability  and  lead  to  the  efficient  development  of 

stable  channels. 
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Figure  Captions 

Figure  1:  . 

(a)  Experimentally  observed  morphology  of  channelled  propagation  obtained  by  imaging  the 

Xe(M)  radiation  (~  IkeV")  transverse  to  the  channel  with  an  x-ray  pinhole  camera  ha\ing  a 
spatial  resolution  of  -  30  pm.  The  procedures  of  the  measurement  are  described  in  Ref  [4}. 
The  incident  248  nm  pulse  had  an  energy  of  -  350  ml  and  a  pulse  width  of  -  270fs.  The  initial 
focusing  action  begins  at  z  =  0.5  mm  and  a  clearly  bifurcated  region  is  visible  at  z  =  1.5  mm. 
See  text  for  discussion. 

(b)  Computed  pulse  evolution  corresponding  to  the  experimental  circumstances  yielding  the 
observed  pattern  shown  in  Fig.  1  (a).  A  strong  bifurcation  of  the  propagating  node  is  seen,  an 
outcome  matching  the  experimental  result. 

Figure  2: 

Stability  map  with  coordinates  (rj,  po).  The  eigenmode  curve,  pe.o('n)  is  shown.  A  full  description 
appears  in  Ref  [1]  .  See  text  for  discussion  of  designated  loci. 

Figure  3:  •  .  • 

(a)  Morphology  of  channel  seen  with  adjustment  of  the  incident  focal  volume  to  a  region  havmg  an 

appreciable  longitudinal  plasma  density  gradient.  The  procedures  of  measurement  were  identical 
to  those  used  in  Fig.  1(a).  A  single  channel  is  seen.  See  text  for  discussion. 

(b)  Computed  channel  evolution  for  a  pulse  corresponding  to  the  experimental  conditions  yielding 
the  result  in  Fig.  3(a).  A  single  channel  is  observed  in  conformance  with  observation. 
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ABSTRACT 


A  mode  of  detonation  for  neutron  stars  is  derived  based  on  a  physically  defined  power  residue  counting  system  used 
to  construct  a  particle  mass  scale  that  reflects  supersymmetric  fermionyboson  pairing.  Neutron  stars  with  specific  mass 
values  (Mg)  can  become  unstable  toward  decay  into  an  ensemble  of  highly  energetic  pairs  of  supersymmetric  particles 
with  an  energy  release  whose  upper  bound  can  approach  the  rest  mass  energy  of  the  original  bo(fy.  The  subsequent 
scattering  of  these  particles  provides  a  way  to  account  for  (1)  the  high  radiative  yield  (~10^^-10^’  J)  of  gamma-ray  bursts, 
(2)  the  large  variability  of  their  observed  energy,  and  (3)  the  expected  collateral  production  of  ultrahigh  energy  cosmic 
rays.  The  joint  action  of  (a)  supersymmetric  classification  and  (b)  a  new  arithmetic  selection  rule  designated  as 
congruence  conservation  governs  the  condition  of  instability  and  the  dynanuc  participation  of  the  magnetic  monopole 
obviates  baiyon  number  conservation.  The  decay  patterns  of  certain  elementary  particles  and  putative  supermassive 
forms  of  cosmological  dark  matter  may  be  subject  to  similar  conditions.  Specifically,  the  analysis  predicts  the  existence 
of  a  super-light  boson  neutrino  vj  that  plays  a  fundamental  role  in  weak  interactions  and  it  is  conjectured  that  the 
phenomenon  of  lepton  number  conservation  is  isomorphic  to  the  supersynunetric  conservation  of  congruence.  The  vi 
boson  also  provides  (i)  a  potential  basis  for  the  CP  violation  observed  in  dec^  and  (ii)  an  alternative  description  of 
the  attribute  of  strangeness.  These  interpretations  support  the  conclusion  that  the  observations  of  gamma-tt^  bursts 
signal  the  existence  of  a  radically  new  organizing  principle  of  physical  interactions  that  defines  a  general  condition  of 
stability  for  massive  particles  and  expresses  the  unstable  case  with  the  copious  production  of  energetic  emissions. 


I.  Introduction 

The  spectacular  energy  release  (~10^’-10^’  J),  flequency 
of  occurrence,  high  variability  in  outyut  (~10*-10^,  and 
spectral  characteristics  (~10'‘-1  McV)  rf  gamma-ray  bursts 
currently  lack  a  unified  description  [1-3].  Physical  evidence 
is  also  extant  which  suggests  the  correlation  of  these 
observations  with  other  energetic  cosmic  phenomena  [4,5]. 
Various  pictures  of  these  enormously  powerful  radiative 
events  plausibly  involve  neutron  stars  in  processes  that  give 
rise  to  the  production  of  super-relativistic  charged  particles 
whose  collisional  and  magnetic  Bremsstrahlung  is  the 
observed  emission.  The  present  stutty  examines  a  specific 
mechanism  which  enables  an  individml  neutron  star  (NS) 


with  mass  Mg  to  undergo  spontaneous  detonation  and 
produce  an  array  of  superenergetic  particles  with  an 
aggregate  yield  comparable  to  the  rest  mass  energy  (MgC^) 
of  the  original  body. 

Basically,  it  is  found  that  the  system  can  become 
unstable  if  appropriate  conditions  governing  its  structure  are 
satisfied.  These  conditions  are  defined  by  the  existence  of  a 
new  form  of  numerical  (arithmetic)  selection  rule  that 
reflects  both  the  specific  magnitude  of  the  mass  and  its 
supersymmetric  (fermi/bose)  classification.  An  isolated 
elementary  event  of  this  nature  considerably  simplifies  the 
circtunstances  of  a  gamma-ttty  burst  in  comparison  to  other 
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descriptions  [1-4],  since  the  participation  of  another 
exceptional  astrophysical  body  (e.g.  neutron  star  or  black 
hole)  is  not  necessary.  Furthermore,  with  the  customary 
mass  limits  [6,7]  for  neutron  stars,  it  also  produces  an 
energy  yield  which  conforms  to  the  observed  range.  As 
indicated  below,  these  considerations  may  also  be  relevant 
to  the  imderstanding  of  the  phenomenon  of  lepton  number 
conservation,  CP  violation  in  K®  decay  [8,9],  the  attribute 
of  strangeness,  and  reactions  involving  putative 
supermassive  forms  of  cosmological  dark  matter  [10-12]. 

II.  Discussion 

A  physically  anchored  scale  invariant  system  of 
enumeration  has  been  developed  [11,12]  which  is  capable  of 
representing  the  particle  mass  spectrum  in  a  way  that 
explicitly  respects  supersynunetric  fermion/boson  pairing 
[12].  The  manifest  scide  invariance  of  the  counting 
procure  [11,12]  confers  on  it  the  capacity  to  represent  a 
particle  with  an  arbitrarily  large  mass,  therd>y  enabling  a 
description  of  neutron  stars  or  even  more  massive  objects. 
A  specific  result  of  the  analysis  of  the  mass  scale  [12]  was 
the  prospective  identification  of  a  fimdamental 
supersynunetric  pair  composed  of  a  light  (rUa  ~0.808  tneV) 
fermion,  presumably  representing  the  electron  neutrino,  and 
its  supersymmetric  partner,  a  heavy  (mp~1.22  x  10“  eV) 
sneutrino  which  serves  as  a  candidate  for  an  irtunobile  form 
of  cosmological  dark  matter  [10-12].  The  predicted 
characteristics  [1 1,12]  of  the  sneutrino  (m^)  are  such  that  it 
can  account  for  the  existence  [12]  of  a  cosmological 
constant  with  a  value  of  A  =  0.7.  Since  the  masses  nia  and 
trip  of  these  two  particles  physically  define  the  mathematical 
bases  of  a  dual  counting  operation,  they  were  nominated  as 
the  basis  of  a  fimdamental  organizing  principle  through  their 
roles  as  the  twin  generators  of  tte  mass  scale.  The 
supersymmetric  fermion/boson  classification  that  emerged 
[12]  was  also  found  to  be  consistent  with  the  values  of  the 
masses  computed  in  prior  work  [11]  for  both  the  electron 
and  the  protoa 

The  physically  established  counting  system  [12],  which 
represents  particle  masses  with  a  normalized  integer 
amplitude  B,  is  expressed  with  power  residues  and 
congruences  [13]  (tf  the  form 

gi«  sB(modPa).  (1) 

in  which  P.  is  a  prime  modulus,  g.  is  a  primitive  root  [14- 
17]  of  Pa,  and  t.  is  the  integral  index  parameter  (power) 
whose  range  is  determined  1^  the  Euler  function  ^o)  as 
0<ta<  Pa-1.  Furthermore,  since  the  ^>ecific  magnitudes  of 
Pa  and  ga  are  independently  and  directly  founded  [12]  on  an 
interlocking  set  of  observational  data  involving  the  Hubble 
dimension,  the  mass  of  the  universe,  chief  characteristics  of 
the  cosmic  microwave  background,  the  magnitude  of  the 
gravitational  constant  G,  and  the  value  of  the  fine-structure 
constant  a,  the  power  residue  formulation  combines  a 


physically  unique  definition,  based  on  both  cosmic 
phenomena  and  microscale  couplings,  with  a  strict 
mathematical  imperative,  the  primitive  root  relation 
between  Pa  and  g*.  The  integer  also  physically 
represents  the  mass  of  the  magnetic  monopole  [11],  a 
crucial  condition  in  relation  to  the  stability  of  a  neutron  star, 
since  the  dynamic  presence  of  that  particle  [18]  enables  the 
occurrence  of  processes  which  violate  baiyon  number 
conservation  [19]. 

Consider  an  isolated  uncharged  rotationally  static  mass 
Mo  and  proportionally  corresponding  [11,12]  normalized 
mass  parameter  Bo.  From  Eq.  (1),  we  can  write  the  explicit 
relation 

®o  ~  ^ao^a  (2) 

for  appropriate  integer  indices  (too,  kao).  Since  previous 
studies  [11,12]  have  identified  ga  as  an  even  number  and  Pa 
as  an  odd  prime,  the  correspondence  Bo  (odd,  even)  and  koo 
(odd,  even)  for  the  modulus  index  holds.  Furthermore, 
since  the  supersymmetric  classification  [12]  generally 
partitions  fermions  and  bosons  with  the  parity  of  ko,  any 
mass,  including  a  neutron  star,  can  have  either  fermi  orbose 
character  depending  upon  the  corresponding  parity  of  the 
mass  parameter  Bo.  llie  previous  study  of  the  mass  scale 
[12]  1^  to  the  conclusion  that  ka  (even/odd)  corresponds  to 
(fermi/bose),  respectively. 

We  now  derive  ^e  conditions  necessary  for  the 
elementary  process  of  decay  of  a  system  with  mass  Mo  into 
two  energetic  particles  with  masses  m,  and  and 
corresponding  integer  mass  parameters  B^  and  By  in  the 
reaction 

Mq  m^  +  my .  (3) 

Using  the  normalized  momenta  (A^  A,,  Ao),  masses 
(Bx,  By),  and  total  energies  (Cx,  Cy,  Co)  defined  generally  in 
earlier  woik  [11],  the  corresponding  relativistic 
energy/momentum  relations 

Bo=Co.  (4) 

®x  ~  ^x » 
and 

Ay  +  By  =  Cy  (6) 

hold  in  the  rest  firame  (Ao=0)  of  the  particle  M..  With  the 
conditions  for  conservation  of  linear  momentum  and  energy 
given  by 

Ajj  +  Ay  —  0  (7) 
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and 

Cx  +  Cy  =  Cq  =  Bq  ,  (8) 


The  restriction  of  Eq.(14)  to  decay  pathways  involving 
supersymmetric  pairs  defines  an  important  subclass  of 
solutions.  The  properties  of  these  solutions  can  be  examined 
by  utilizing  the  conjecture  that  the  sum  I,  of  the  masses  of 
all  supersymmetric  pairs  is  governed  by  the  statement 


the  solution  of  Eqs.  (5)  and  (6)  yields  the  individual  particle 
energies 


Zs  =0(niodPa),  (15) 

the  outcome  suggested  by  the  study  [12]  of  the  systematics 
of  the  power  residue  mass  scale.  Therrfore,  if  we  assume 
that  Bx  and  By  represent  a  supersymmetric  pair,  w  have 


and 


Bx  +By  s  0(modP(x), 


or  equivalently, 

(10) 

®y  =  ^^xyPa  “  Bx 


(16) 

(17) 


Since  Cx  and  Cy  are  constrained  to  have  integral  values  [11], 
we  require  that 


for  some  integers  k,.  and  k_,  a  condition  that  yields 

±^B|-Bj^=(2k±-Bo)Bo.  (12) 

With  the  further  definition  of  the  effective  modulus 
parameters 

ko±=2k±-Bo,  (13) 

Eq.  (12)  can  be  expressed  as  the  quadratic  congruence 

B2.B;(mod  B,).  (14) 

Equation  (14)  iqtresents  an  arithmetic  selection  rule 
relating  the  thm  mass  parameters  Bx,  By.  and  Bo  for  the 
two-body  decay  given  by  Eq.(3).  Consequently,  the 
permissible  solutions  of  Eqs.(9)  and  (10)  are  subject  to 
constraints  that  are  dependent  upon  the  detailed  arithmetic 
structures  of  Bx,  By,  and  Bo.  FromEq.  (13)  it  is  immediately 
seen  that  the  supersymmetric  dassification  of  mass 
parameters  causes  two  classes  of  solutions  corresponding  to 
bose  and  fermi  species  to  exist  Specifically,  we  have  the 
effective  modulus  parameters  kot  (even/odd)  for  Bo 
(fermi/bose),  respectively.  Many  characteristics  of  the 
solutions  of  Eq.(14)  are  known  [13,14,20,21).  For  example, 
any  integer  [14,20]  represented  by  the  form  ^*-By*  is  either 
odd  or  a  multiple  of  4,  and  general  classes  of  solutions  of 
the  congruence  are  available  in  both  cases  [21]. 


for  some  odd  integer  k^,  since  Bx  +  By  must  be  odd. 
Explicitly  forming  the  quadratic  difference 

Bx  “By  =  2kxyBxPa  “kxyP<x,  (18) 

it  follows  that 

b|  s  B^(modPa).  (19) 

Hence,  the  combination  of  Eqs.  (14)  and  (19)  yields  the 
additional  relation  [17] 

Bx  sBy(modNao)  (20) 

in  which  the  modulus  Noo  is  the  least  common  imiltiple  of 
Bo  and  Since  P.  is  a  prime,  the  minimum  value  of  Noo 
attainable  is  P,.  Physical^,  of  course,  this  decay  can  only 
occur  if  Bo^  P«. 

We  now  summarize  the  development  given  above.  The 
conditions  stated  in  Eqs.  (13)  and  (14)  define  a  general 
arithmetic  selection  rule  correspond^g  to  amplitudes  for 
two-bo(fy  decay  of  fermi  and  bose  partides  that  is  consistent 
with  supersymmetric  classification.  The  restriction  of  this 
selection  rule  to  pathways  involving  the  production  of 
supersymmetric  fermion/boson  pairs  is  represented  by  Eqs. 
(14),  (19),  and  (20).  We  note  that  fulfillment  ai  these 
conditions  is  explicitly  consistent  with  integral  solutions  for 
Ax  and  Ay  and  the  condition  Ax+Ay=0,  a  statement  whose 
proof  is  readily  demonstrated  by  di^  subriitmioo  of  Eqs. 
(9)  and  (10)  into  the  coneqxinding  Eqs.  (S)  and  (6).  Hence, 
the  identical  selection  rules  govern  both  energy  and 
momentum.  We  designate  this  selection  rule  generally  as 
the  conservation  of  congruence. 
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A.  Conservation  of  Congruence  and  Angular 
Momentum 

The  consequences  of  the  joint  conservation  of 
congruence  and  angular  momentum  on  the  decay  pathways 
of  an  isolated  massive  particle  are  now  evaluate.  At  a 
further  point,  the  role  of  baryon  number  conservation  and 
related  implications  concerning  lepton  number  conservation 
[19]  will  be  explored. 

1.  Fermi  Mass  Parameter 

In  this  case,  since  Bo  and  the  corresponding  modulus 
parameters  ko±  are  even,  is  perforce  even  on  the 

basis  of  the  supersymmetric  classification  [12].  Indeed,  as 
noted  above  [14],  is  generally  a  multiple  of  4. 

Hence,  B*  and  By  are  either  (a)  both  even  (fermi)  or  (b)  both 
odd  (^se).  Since  two  fermions  or  two  bosons  cannot 
represent  the  angular  momentum  state  of  a  single  fermion, 
the  two-body  decay  paths  into  either  fermions  or  bosons  is 
strictly  forbidden  by  the  conservation  of  angular 
momentum.  The  remaining  possibility  for  two-body  decay, 
which  leads  to  a  final  state  composed  of  a  fermion/boson 
pair,  is  barred  bv  supersymmetry  fi'om  Eqs.  (13)  and  (14). 
Hence,  all  two-bodv  pathways  for  the  decay  of  fermi  species 
(Bo)  are  forbidden. 

It  is  important  to  recognize  that  this  result  is  entirely 
dq>endent  upon  the  classification  of  fermions  with  even 
ma«  parameters  in  the  construction  of  the  mass  scale  [12], 
an  outcome  that  arose  perforce  fi’om  consideration  of 
supersymmetric  fermion/boson  pairing  with  respect  to  the 
identification  of  the  mass  of  the  electron  neutrino  (v,).  Had 
that  assignment  been  the  reverse,  with  the  specification  of 
odd  mass  parameters  for  fermions,  the  reaction  shown  by 
Eq.  (3)  would  be  allowed  for  a  fermion/boson  pair,  since  the 
integer  parity  condition  demanded  Eqs.  (11)  through  (14) 
and  the  angular  momentum  requirement  could  be 
simultaneously  satisfied.  This  alternative  assignment  would 
have  Lifted  the  bar  imposed  by  supersymmetry.  Therefore, 
the  conclusion  is  reached  that  all  two-bodv  amplitudes  are 
forbidden  for  fermi  character  Bo  by  the  combined  demands 
of  conservation  of  supersymmetric  congruence  and  angular 
momentunL  Specifically,  we  arrive  at  the  import 
outcome  that  the  sinele  fermion/boson  pair  charnel  is 
blocked  by  supersymmetry. 

2.  Bose  Mass  Parameter 

Boson  decay  gives  B,  odd  and  Bx^-By’  odd  with  the 
direct  consequence  that  B.^  and  By^  are  necessarily  of 
opposite  parities.  It  follows  that  the  only  supersymmetry 
allowed  two-body  final  state  is  a  fermion/boson  pair. 
Ifowever,  since  tUs  configuration  carmot  satisfy  the  an^ar 
momentum  requirements  for  a  bose  system,  the  amplitude  is 
forbiddea  Again,  the  joint  action  of  the  conservation  of 
supersymmetric  congruence  and  angular  momentum  rules 
out  an  amplitude  for  two-body  decay.  It  is  important  to  note 
fiom  Eqs.  (14)  and  (20),  however,  that  for  a  system  of  bose 


character,  the  single  fermion/boson  channel  is  allowed  bv 
supersymmetry  in  contrast  to  the  fermi  case  with  B,,  even. 

Table  I  smnmarizes  the  conclusions  governing  two- 
body  amplitudes  for  both  fermi  and  bose  particles.  The 
astonishing  result  is  that  aU  two-body  chatmels  are 
forbidden  by  the  joint  fulfillment  of  congruence  and  angular 
momentum  conservation.  Obviously,  a  bold  contradiction  to 
the  clear  e.>q>erimental  observation  [22]  of  two-body  decay 
charmels  for  a  wide  range  of  elementary  systems 

(e.g.  n*  — >  )  is  presented  by  this  highly  unexpected 

outcome.  This  important  conflict  is  fuUy-resolved  below  in 
Section  n.A.3.b  through  the  development  of  a  specific 
physically  motivated  modification  of  Eq.(14). 


B, 

Two-Body 

Ampbude 

SdectioQ  Fries 

Angular  MomeG&m 

Si^ersymmetiy/ 

CoQgruence 

Fcnrri 

fi+f. 

Forbidden 

/dlowe-i 

Fenni 

f«+b. 

AIlow«d 

Forbidden 

Fermi 

bx  +  b. 

Forbidden 

Allowed 

Bose 

Allowed 

Forbidden 

Bose 

f,+  b. 

Forbidden 

Allowed 

Bose 

b,  +  b. 

.Allowed 

Forbidden 

Table  I:  Summary  of  two-body  amplitudes  for  the  decay  of 
a  fermi/bose  particle  with  mass  parameter  Bo.  The  notation 
fi/b)  designates  fermi  or  bose  particles  for  the  final  state 
m*-t-my  given  by  Eq.  (3).  The  joint 

requirements  of  congruence  and  angular  momentum 
conservation  block  all  two-body  chatmels. 


3.  Supersymmetric  Pair  Pathways  and  Lepton/Baryon 
Conservation 

a.  Three-Body  Amplitudes 

The  addition  of  a  third  body  (m  J  consistent  with  the 
requirements  of  angular  momentum  conservation,  which 
converts  Eq.(3)  to  the  three-body  process 

my  +  m^ ,  (21) 

enables  all  six  chatmels  in  Table  I  to  become  allowed.  In 
interesting  contrast  to  the  fundamentally  unsolvable  nature 
of  the  standard  3-body  problem  [23]  and  the  complexities  of 
radiative  corrections  to  normal  channeb  of  P-decay  [24],  is 
easily  seen  that  the  energies  and  momenta  characterizing  the 
three-body  final  state  can  be  exactly  determined  on  the 
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integral  cnergy/momcntum  lattice.  This  fact  can  be 
demonstrated  from  known  properties  of  the  solution  [11]  of 
the  elementary  case  which  involves  values  of  Bx,  By,  and  B, 
given  by  odd  primes.  Subject  to  the  physically  necessary 
condition  Bx+]^+Bx<  Bo.  this  case  admits  the  possibility  of 
only  a  single  solution  for  the  three-body  Hn^  state.  It  is 
readily  shown  that  the  specific  solution  for  B^  By,  and  B, 
prime  corresponds  to 


and 


B^+1 

2 

(22) 

B^  +1 

2  ’ 

(23) 

Bi+1 
’  2  ’ 

(24) 

under  the  condition  that  the  subsidiary  statements 


2Bo=B|-hB^+B|-h3  (25) 

and 

Bl+Bj^Bl+l^Bl-Bj+2  (26) 

are  fulfilled.  Equation  (23)  represents  energy  conservation 
while  the  inequalities  [23]  expressed  in  Eq.(26)  enable 
satisfaction  of  the  conservation  of  linear  momentum.  For 
explicit  confirmation  of  solvability,  the  values  (Bx,,  By,  BJ 
=  (29,  23,  31)  uniquely  satisfy  the  conditions  of  Eqs.  (25) 
and  (26)  for  Bo=1167.  Since  all  four  mass  parameters  are 
odd,  this  rmmerical  example  would  rqrresent  the  decay  of  a 
boson  into  three  bosons,  an  outcome  that  also  satisfies  the 
demand  fw  angular  momentum  conservatiotL  This  simple 
case  demonstrates  the  existence  of  arithmetic  constraints  on 
the  set  mass  parameters  {B,^  By,  Bx,  Bo}  which 
determine  the  condition  for  solvabilify. 

The  generalization  to  composite  mass  parameters, 
which  gives  the  possibilify  of  a  multiplicity  of  solutions,  is 
readily  foanulat^  In  o^r  to  obtain  a  high  densify  of 
solutions  per  unit  normalized  energy  [11]  of  the  individual 
particles  m.,  niy,  and  nix,  the  corresponding  mass  parameters 
Bx,  By,  and  Bx  must  be  composite  high  k-factorable  low  z* 
smooth  integers  [26]  whose  arithmetic  factor  structures  are 
mutually  sifoject  to  numerical  constraints.  Basically,  these 
constraints  must  be  such  that  the  number  of  representations 
N(Bo)  of  the  integer  Bo,  which  denotes  the  total  energy  of 
the  system,  through  the  expression 

+  Cy  +  C2  ~  Bo  >  (^T) 

is  high.  Formally,  we  can  express  Eq.(27)  as 


2Bo—  dx+dy+dx+ 


(28) 


in  which  dx,  dy,  and  d,  respectively  represent  the  appropriate 
divisors  of  Bx^  By^  and  Bx*.  Satisfaction  of  linear 
momentum  conservation  requires  that  the  subsidiaiy 
conditions 


+  dy  — 


+  dv 


(29) 


also  be  fulfilled.  We  observe  that  Eqs.(28)  and  (29)  coincide 
with  the  corresponding  statements  in  Eqs.(25)  and  (26)  for 
the  case 


dx=dy=dx=l.  (30) 


The  existence  of  these  arithmetic  conditions  obviously 
carries  important  implications  for  the  systematics  of  the 
particle  mass  scale  [12]  and  the  observed  channels  of  decay. 
For  example,  if  the  sets  of  mass  parameters  (Bo,  B*,  By,  B,} 
are  such  that  Eqs.  (28)  and  (29)  do  not  possess  a  solution  for 
all  physically  possible  choices  of  Bx,  By,  and  Bx,  the  particle 
M«  in  Eq.  (21)  would  behave  as  a  stable  system  like  the 
proton.  In  this  wav.  Eos.  (28)  and  (291  define  the  conditions 
Eoveming  particle  stability. 

Since  the  squares  of  suitably  smooth  integer  mass 
parameters  for  typical  particles  like  the  electron  and  proton 

can  have  more  than  10*^  divisors  [11],  the  requirement  for 
a  high  representation  number  N(Bo)  involves  a  complex  set 
of  relations  concerning  the  arithmetic  structures  of  the  mass 
parameters  whose  precise  nature  remains  to  be  formulated 
[21].  It  may  be  possible  to  consider  Eq.  (28)  as  a  restricted 
class  of  partitions  [27]  of  the  even  integer  2Bo. 

A  physical  example  of  the  three-body  final  state  for  the 
decay  of  a  fermion  is  the  commonly  observed  conversion  of 
a  free  neutron  to  a  proton,  an  electron,  and  an  electron  anti¬ 
neutrino  by  the  3-b^  amplitude 


n->p+e  -t-Vg.  (31) 

This  process  explicitly  exhibits  the  conservation  of  lepton 
number,  an  empiri^y  established  quantity  which  is 
believed  to  be  separately  conserved  for  electron,  muon,  and 
tau  species.  In  principle,  of  course,  the  third  body  (m,)  in 
Eq.(21)  could  also  be  a  photon,  as  it  is  in  the  Dalitz  decay 
[22]  of  the  neutral  pion  (n*^.  In  light  of  the  discussion  given 
above,  the  comparison  between  the  character  of  the 
theoretical  results  illustrated  in  Table  I  and  the 
experimentally  observed  behavior  of  particle  decay  channels 
(e.g.  the  neutron)  suggests  that  supersymmetric  congruence 
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conserv'ation  may  potentially  provide  a  mathematical  basis 
for  the  phenomenon  of  lepton  number  conservation. 
Separate  conservation  of  lepton  number  for  electron,  muon, 
and  tau  amplitudes  would  then  be  direct  consequences  of 
the  detailed  arithmetic  structures  of  the  mass  parameters  of 
those  particles  (e.g.  c",  p',  x",  ,  v^,  Wj )  and  the  ability  of 

the  specific  factors  of  these  parameters  to  satisfy  the 
arithmetic  demands  expressed  by  Eqs.(28)  and  (29)  with 
integers  of  the  conect  parity,  in  addition  to  the  requirement 
for  the  conservation  of  angular  momentum.  This  hypothesis 
is  currently  under  computational  evaluation  with  respect  to 
the  p-decay  of  the  neutron  (Eq.  (31)],  since  the 
neutron/proton  mass  ratio  is  experimentally  Imown  to  high 
precision  [22]  and  candidate  values  for  the  mass  parameters 
of  the  p,  e',  and  v^  are  available  [11,12]. 


b.  Analysis  of  Two-Body  Channels 


Two-body  decay  channels  for  both  fermi  and  bose 
systems  are  well  established  experimentally  for  many 
particles  (e.g  ,  and  A^-^p  +  n").  It  is 

obvious  that  these  observational  facts  stand  in  direct 
contradiction  to  the  conclusions  stated  in  Table  I.  The 
manner  in  which  this  conflict  can  be  resolved  is  now 
described. 

The  resolution  of  this  issue  requires  the  participation  of 
an  additional  interacting  particle.  The  influence  of  a  third 
body  in  the  final  state  of  the  system  with  a  mass  parameter 
of  unity  (B*  =  1),  as  expressed  through  Eqs.(28)  and  (29), 
has  a  magnitude  sufficient  to  nullify  fully  the  arithmetic 
restrictions  given  by  Eq.(14).  Since  this  conclusion  is  valid 
for  Bx,  By,  and  Bo  of  any  numerical  size,  the  minimal  value 
Bx  =1  has  scale  invariant  consequences.  Hereafter,  a  bose 
particle  with  the  mass  parameter  of  unity  will  be  designated 
with  the  symbol  vj. 

The  physical  consequences  of  this  super-light  boson 
can  be  readily  evaluated.  Specifically,  with  Bx  -  1, 
Eqs.(28)  and  (29)  are  respectively  modified  to  read 


and 


=  2(Bo-l) 


(32) 


(33) 


results  which  lead  immediately  to  the  statement 


Bi-B5=(dx-dy)(Bo-l)  (34) 


With  the  identification  of  the  fyctor  (dx-dy)  with  the 
modulus  index 


k  =  dx-dy.  (35) 

Eq.  (34)  can  be  written  as  the  congruence 

®x  —  By(mod(Bo  -1)) ,  (36) 

a  form  which  is  identical  to  Eq.(14)  under  the 
transformation  of  the  modulus  B©  Bq  - 1 .  It  also 
follows  fi'om  Eqs.  (32)  and  (33)  that 


Bo  =dxo  +dyo  +  1.  (37) 

in  which  dxoand  dyo  respectively  represent  specific 

members  of  the  sets  of  divisors  of  the  integers  Bx^  and  By^. 

The  shift  of  unity  in  the  value  of  the  modulus 
completely  eliminates  the  restrictions  listed  in  Table  I 
arising  from  violation  of  the  supersymmetiy/congruence 
selection  rule.  G>nsider  the  known  [22]  weak  decay  process 

(38) 


a  reaction  that  is  forbidden  by  supersymmetry  in  Table  I. 
However,  if  the  reaction  given  by  Eq.(38)  is  physically 
reinterpreted  to  be 


x+ -►p++v^+vi,  (39) 

the  amplitude  becomes  fully  allowed  as  a  three-body 
process  in  accord  with  Eq.(36).  Furthermore,  since  Eq.(24) 
states  that  a  particle  with  a  mass  parameter  of  unity  (Bx  =  I) 
can  only  have  a  vanishing  momentum  (Ax  =  0),  the 
vi  particle  plays  no  role  in  the  distribution  of  momentum. 
Therefore,  the  kinematic  conservation  of  momentum  is 
bbservationally  identical  in  both  Eqs.  (3  8)  and  (39).  The 
particle  V)  acts  as  an  experimentally  invisible  boson 
neutrino  whose  rest  mass  energy  is  given  by  the  quantify 
(Eo),  the  quantify  which  defines  the  minimum  energy 
uncertainty  of  the  system  [11,12].  Hence,  energy 
conservation  is  also  observationally  unaffected  by  the 
coupling  to  the  vi  particle.  Cktnsequently,  the  three-body 
reaction  for  %*  decay  shown  in  Eq.(39)  behaves 
kinematically  as  a  two-botfy^  amplitude  that  would  not  be 
distinguished  from  Eq.  (38)  by  standard  procedures  of 
measurement  In  contrast  the  decay  of  the  neutron  given  in 
Eq.(31)  behaves  kinetically  as  a  conventional  three-body 
process  and  yields  a  continuous  spectrum  of  particle 
energies.  The  amplitudes  governing  the  decay  patterns  of  all 
other  known  particles  [22]  appear  consistent  with  this 
formulatioa  For  example,  the  principal  decay  mode  of  the 
A°hyperon  (A®->p  -i-  «"),  a  fermion,  can  be  likewise 
considered  as  the  reaction 

A®  -♦p-nc“-t-v, . 


(40) 
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This  reassignment  of  the  amplitude  for  A°  decay  in  Eq. 
(40)  also  suggests  an  alternative  interpretation  of  the 
attribute  of  strangeness.  Specifically,  the  low  rates  generally 
exhibited  by  nonleptonic  hyperon  (e.g.  A®,  V,  CT)  decay 
can  be  regarded  to  be  the  result  of  a  suitably  weak  coupling 
to  the  vi  boson  in  processes  of  the  form  shown  in  Eq.  (40) 
for  the  A®  particle.  At  the  same  time,  it  is  important  to  note 
that  the  corresponding  production  of  the  A®  system  by  the 
reaction 

x~+p->A°+K°  (41) 

can  occur  unsuopressed  directly  through  the  strong 
interaction.  This  follows,  since  the  participation  of  four 
strongly  interacting  particles  (n”,  p.  A®,  K®)  allows  the 
arithmetic  conditions  specified  by  ^s.  (28)  and  (29)  to  be 
satisfied.  The  conundrum  of  strong  production  joined  with 
weak  decay  in  the  n  +  p  channel  is  thereby  resolved. 
Further,  the  nonoccurence  of  the  ampUtude 

x-+p->n  +  KO  (42) 

can  also  be  understood  as  a  case  for  which  the  divisors  of 
the  corresponding  mass  numbers  of  the  particles  do  not 
admit  a  solution  of  Eqs.  (28)  and  (29).  Therefore,  the 
behavior  associated  with  the  assignment  of  strangeness  can 
be  alternatively  considered  to  be  the  natural  consequence  of 
(a)  weak  coupling  to  the  V|  boson  and  (b)  the  arithmetic 
requirements  of  Eqs.  (28)  and  (29). 

The  vi  particle  can  be  given  a  supersymmetric 
classification  [12]  in  relation  to  the  mass  scale  expressed  by 
Eqs.  (1)  and  (2).  Specifically,  in  accord  with  the  previous 
development  [12]  of  Eq.(15),  the  corresponding 
supersymmetric  partner  of  the  V]  particle  has  a  mass 
parameter  with  the  value  P«-l.  Therefore,  this  assignment 
predicts  the  existence  of  a  supermassive  fermion  4>i  with  a 
magg  infinitesimally  shifted  firom  that  of  the  monopole 
[11,12],  Overall,  the  characteristics  of  the  particle  mark  it 
as  an  attractive  candidate  for  a  supermassive  component  of 
cosmological  dark  matter  [10-12]. 

Since  the  construction  of  the  supersymmetric  mass 
scale  led  to  the  definition  of  a  four  dimensional  index  space 
[12],  the  corresponding  indices  of  the  vi  particle  associated 
with  the  power  residue  formulation  (cf.  ^s.  (1)  and  (2))  can 
be  given.  They  are  t«=tp=0  for  the  power  indices  and 
ka=kp=0  for  the  modulus  indices.  Hence,  the  vi  boson 
physically  defines  the  origin  of  the  index  space  governing 
the  classification  of  particle  masses.  It  is  also  simply  shown 
[12]  that  the  twn  power  indices  corresponding  to  the 
supersymmetric  partner  ^  are  degenerate  with  values  of 
t.=tp=(P.-l)/2. 

The  presence  of  the  vi  boson  enables  a  three-body 
amplitude  to  behave  obseivationally  as  a  two-body  process. 
For  example,  if  a  neutral  particle  were  to  have  two  bosonic 
modes  of  decay,  one  involving  three  bosons  (e.g. 
and  another  yielding  two  bosons  and  an  unobserved 


Vi  particle  (e.g.  7t*,Ji‘,vi),  experimentally  distinct  3n  and 
2n  chaimels  with  different  lifetimes  would  be  identified.  In 
this  case,  the  two  decay  modes  are  simply  related  by  a  n/  vi 
e.\change.  Interestingly,  for  K  mesons,  this  experimental 
characteristic  is  precisely  the  signature  of  direct  CP 

violation  in  K®  decay  [8,9].  Therefore,  the  existence  of  a  vi 

boson  could  provide  a  mathematical  basis  for  CP  violation 
which  is  founded  on  the  power  residue  formulation  of  the 
particle  mass  scale  [12]. 

c.  Baryon  Conservation 

Baryon  conservation  naturally  forbids  the  detonation  of 
a  neutron  star  in  a  process  which  calls  for  a  massive 
decrease  in  baryon  number  to  occur.  However,  this 
customary  restriction  can  be  obviated  by  the  dynamic 
participation  of  the  magnetic  monopole  [18,19],  foe  key 
particle  promoting  baryon  nonconservation.  Since  foe 
theoretic^  picture  underlying  foe  construction  of  foe  mass 
scale  [12]  suggests  that  foe  mass  numbers  associated  with 
supers>mmetric  fennion/boson  pairs  generally  satisfy  foe 
congruence  given  by  Eq.  (16),  we  can  generally  write 

Bf  +  Bb  s  0(modPa  ) ,  (43) 

in  which  the  subscripts  f  and  b  denote  fermi  and  bose 
particles,  respectively.  We  recall,  since  Bft-Bb  must  be  odd, 
this  statement  specifies  that  foe  masses  physically  sum  to  an 
odd  multiple  (nPa,  n=l,  3,  S, ...)  of  foe  monopole  mass  [12]. 
Hence,  for  foe  case  n=l,  foe  sum  of  foe  masses  of  an 
arbitrary  supersymmetric  pair  would  corresiwnd  exactly  to 
foe  monoDole  mass  shell  fP-).  In  turn,  this  suggests  that  foe 
barvon  nonconserving  decay  of  a  massive  object  can  occur 
thrxHigh  the  production  of  an  ensemble  of  supersymmetric 
pairs,  or  equivalently,  by  a  process  that  goes  directly 
through  the  monopole  state  for  every  particle/pair.  The 
monopole  would  then  pl^  foe  role  of  the  fiindamental 
modulus  connecting  foe  individual  decay  channels  through  a 
barvon  nonconserving  amplitude.  Therefore,  with  this 
interpretation,  foe  ability  to  involve  dynamically  foe 
monopole  state  through  the  satisfaction  of  (43)  obviates 
foe  restriction  of  barvon  conservation. 

We  conclude  that  two  conditions  must  be  fulfilled  to 
permit  the  existence  of  a  baryon  nonconserving  chaimei  for 
the  decay  of  an  isolated  body  with  mass  They  are  (1) 
the  presence  of  three  or  more  particles  in  foe  final  slate  and 
(2)  the  availability  of  an  amplitude  invohing  the  production 
at  supersymmetric  pairs  which  couple  directly  through  foe 
monopole  mass  shell  P,. 

Allowed  decay  chaimels  for  both  bose  and  fermi 
character  systems  can  now  be  classified.  Under  foe 
combined  requirements  of  Eqs.  (13),  (14),  (19),  (20),  (28), 
and  (29)  it  is  immediately  seen  that  bose  species  can 
undergo  baryon  nonconserving  decay  into  chaimels 
involving  supersymmetric  pairs  (fx  +  by)  m  reactions  of  foe 
general  form 
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Bo(bose)  ->  (fx  +  by  )2tn-l  +  (fz)2n-l 

and 

Bo(bose)->(fx +by)2m+(bz)n-l.  ('»5) 

in  which  m  and  n  are  nonzero  integers  (m,n=l,2,3, ...),  fx 
and  by  denote  respectively  a  supersynunetric  particle  pair 
with  corresponding  mass  parameters  Bx=Bf  and  By=Bb,  and 
fjbt  designates  an  appropriate  fermion/boson  companion 
with  mass  parameter  Bi.  The  two  decay  pathways  given  by 
Eqs.  (44)  and  (45)  can  satisfy  (a)  the  specific  requirements 
for  conservation  of  angular  momentiun,  linear  momentum, 
and  energy,  (b)  the  general  consequences  of  the  arithmetic 
constraints,  and  (c)  the  need  to  proceed  on  the  monopole 
mass  shell  for  the  (f*  +by)  supersymmetric  pairs. 

The  corresponding  conditions  are  now  determined  for 
the  decay  of  a  system  with  fermi  character.  A  final  state, 
mntaining  a  Single  pair  is  forbidden  by  supersymmetry,  as 
shown  in  Table  I,  and  the  two-pair  chaimel  is  blocked  by  the 
conservation  of  angular  momentum..  However,  the  presence 
of  an  additional  bose  or  fermi  particle  (b,  or  f*,  respectively) 
in  the  final  state,  as  discussed  above  in  Sections  n.A.3..a  and 
n.A.3.b,  obviates  these  restrictions  and  reactions  involving 
the  production  of  supersymmetric  pairs  (fx+by)  conforming 
to  Eq.  (43)  of  the  general  form 

Bo(fermi)  (fx  +by)2m-l  +  (bz)n 
and 

Bq  (fermr)  — >  (fx  +  by  )2m  (fz)2n— 1 

with  m  and  n  integral  (m,n=l,2,3,...)  become  allowed.  In 
parallel  with  the  bose  case  described  above,  these  two 
rhanm-ls  allow  the  fill!  Satisfaction  of  all  the  corresponding 
constraints  for  fermi  species  Bg. 

The  decay  charmels  given  by  Eqs.(44),  (45),  (46),  and 
(47)  represent  amplitudes  for  baryon  nonconserving 
reactions  which  yield  supersymmetric  (fx+by)  pairs  for  both 
fermi  and  bose  character  systems.  Spe^cally,  for  the  case 
nFrF=l,  Eqs.  (44)  and  (46)  can  be  put  in  direct 
correspondence  with  the  statements  expressed  in  Eqs.  (28), 
(29),  and  (43).  Consequently,  in  the  discussion  below,  we 

will  assume  t^  the  pathways  for  the  decay  of  a  neutron  star 
involve  the  disintegration  into  an  ensemble  of  energetic 
an^fsvmmetric  pairs  and  associated  particle  companions 
whose  subsequent  electromagnetic  scattering  gives  rise  to 
the  phenomenon  of  gamma-ray  bursts. 


B.  Detonation  of  a  Neutron  Star 

We  now  direct  our  attention  to  the  general  systematics 
that  ai^ly  to  the  detonation  of  a  neutron  star  with  mass  Mg 
and  corresponding  mass  parameter  Bg.  Since  the 
congruence  expressed  in  Eq.  (1)  can  represent  any 
normalized  physical  magnitude  [1 1,12],  it  can  express  both 
the  energies  {Q}  and  mass  parameters  {Bj}  the  set  of 


particles  {X;}  involved  in  the  neutron  star  decay.  The 
normalized  integer  energies  of  these  particles  (Q),  which 
are  obtained  by  the  calculational  procedures  previously 
described  [1 1],  are  assumed  to  give  the  sum 

Co=ZC.  (48) 

under  the  constraint  of  momentum  conservation,  namely, 

ZA.  =0  (49) 

1 

Therefore,  the  parameter  C<i  represents  the  total  energy  of 
the  energetic  decay  products  in  the  process 

NS^{Xi}  (50) 

and  the  individual  energies  Q  can  assume  even  and  odd 
values  according  to  the  factor  structures  of  the  mass 
parameters  of  the  particles  involved  in  the  decay  (11]. 
Table  n  summarizes  the  parities  of  the  energy  parameters  Q 
that  can  be  attained  for  different  individual  bose  and  fermi 
particles  based  on  the  possible  arithmetic  structures  of  the 
conesponding  mass  parameters. 


Particle  Mass 

Parameter  B 

Energy  C:  Parity  | 

even 

odd 

Bose  (odd)  (2°) 

X 

/ 

•/ 

Fermi  (even)  (2^) 

.  /■' 

X 

Fermi  (even)  (2“,ni2) 

,  /’ 

/' 

V 

X!  Forbidden  parity  of  energy 


/  Allowed  paii^  of  energy 

Table  II;  Parities  of  the  energy  parameters  C  for  different 
arithmetic  structures  of  the  mass  parameter  B.  Bose  species 
always  yield  odd  values;  fermi  systems  give  rise  to  two 
classes  depending  upon  whether  the  necessary  fktor  of  2 
occurs  eitter  to  the  first  (n=l)  or  to  a  higher  (n22)  power  in 
the  corresponding  mass  parameter. 


In  order  to  conserve  energy  for  mass  paiameiefs  B.  of 
either  fermi  or  bose  character,  it  follows  that  the  energy  sum 
Ca  given  1^  Eq.(48)  must  be  able  to  assume  a  cone^Moding 
even  or  odd  parity,  respectively.  Inspection  of  the  entries  in 
Table  n  demonstrate  that  this  requiremeiM  can  be  satisfied, 
but  that  its  fulfillment  will  impose  a  complex  conanint  on 
the  munber  of  pairs  involved  and  the  detailed  factor 
structure  of  the  mass  parameters  of  the  participating 
particles.  We  encounter  a  basic  arithmcsic  question 
concerning  the  representation  of  a  given  number  B,  by  a 
sum  of  terms  whose  forms  are  determined  by  the  teiaiivistic 
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energy/moroentum  relation  and  the  specific  factor  structures 
of  the  particles  produced  in  the  final  state.  Since  the  gross 
physic^  properties  of  the  participating  particles  are 
probably  unknown,  since  the  arithmetic  structures  of  the 
corresponding  mass  parameters  are  certainly  unknown,  and 
since  consideration  of  a  neutron  star  requires  the  evaluation 
of  this  issue  at  an  enormous  number  scale,  no 

general  analytic  solution  to  this  problem  is  presently  at 
hand.  This  matter  clearly  defines  a  central  question  for 
future  study  that  is  intimately  connected  with  the  detailed 
character  of  the  particle  mass  scale. 


R2mo  =Bn(modPa)  (54) 

for  a  suitable  modulus  parameter 

kCT  =  -(2m  +  s).  (55) 

Three  basic  characteristics  of  the  bose  system  for  the 
pure  supersymmetric  chaimel  (m,l)  are  inunediately 
apparent.  Consider  the  system  for  s=0.  (A)  For  all  mass 
parameters 


C.  Instability  of  a  Neutron  Star 


Bn  S  2P(x 


(56) 


Figure  (1)  schematically  illustrates  the  energetic 
considerations  associated  with  the  decay  of  a  neutron  star 
with  bose  character  Bn  into  supersymmetric  pairs  in  accord 
nith  the  final  states  specified  by  Eq.  (45).  This  process  is 
the  sole  possibility  for  disintegration  purely  into  ensembles 
of  supersymmetric  pairs,  the  condition  defined  by  (m,n)  = 
(m,l).  For  this  channel,  the  total  energies  Cimo  for 
ensembles  of  zero-momentum  pairs  are  illustrated  (e.g. 
C2a,  C40).  With  the  assumption  of  the  pairwise  conservation 
of  momentum  in  Eq.(47)  that  expression  is  simplified  and 
the  corresponding  energies  can  be  readily  computed  [11]. 
An  exact  energy  resonance  €‘3^1  modulo  (s)  the  monopole 
mass  (Pa)  for  the  general  boson  case  corresponding  to 

®n  ~  ^2mo  (5^) 


is  also  shown.  Taking  into  account  the  rest  energy  of  the 
2m-pairs  (2mPa)  and  defining  the  kinetic  energy  release  as 


R*  =r* 
xmcj  ^Imcy 


-2inPa , 


(52) 


we  can  write  the  relation 


R5,„,=Bn-(2m+s)Pa.  (53) 

With  the  parities  fi>r  the  energy  parameters  illustrate  in 
Table  n.  the  integers  s  can  assume  all  non-negative  values 
(s=0,  1,  2,  ...).  It  is  important  to  note  that  the  known 
properties  [11]  of  the  solutions  of  the  energy  parameters 
C*2bw  are  clearly  such  that  the  corresponding  kinetic 
energies  R'^  can  be  super-relativistic. 

Equation  (53)  represents  the  condition  of  instability  for 
a  neutron  star  with  bose  character  for  the  decay  clu^el 
(m,n)  =  (m,l)  defined  by  Eq.(45)  in  terms  of  (1)  the  neutron 
star  mass  B,,  (2)  the  total  kinetic  energy  of  the  2m 
supersymmetric  pairs,  the  value  of  which  is  dependent  upon 
the  exact  arithmetic  structure  of  the  mass  parameters  of  the 
particles  constituting  these  pairs  [11,12],  and  (3)  the 
monopole  mass  represented  by  the  prime  Pa.  This  relation 
can  be  summarized  by  the  statement 


the  system  is  stable.  The  value  2Pa  is  excluded  in  Eq.  (56), 
since  B^  for  bose  character  is  odd  and  2Pa  is  even.  Hence, 
R*2o>0  with  a  minimum  numerical  value  determined  by  the 
factor  structures  of  particles  involved,  a  condition  consistent 
with  the  data  appearing  in  Table  H.  (B)  Small  changes  in  the 
value  of  Bn  can  lead  to  relatively  large  variations  in  R*2iiki, 
the  integer  representing  the  energy  yield  of  the  decay.  This 
arises  naturally,  since  the  modulus  parameter  ka  is  tjpically 
expected  to  have  a  very  large  magnitude  (ko~10^*).  As 
shown  in  Fig.  (1),  this  situation  produces  a  dense  set  of 
overlapping  sequences  of  points  (C'2n«,)  arising  from  all 
orders  up  to  2m,  a  condition  that  enables  k,,  to  experience  an 
extremely  wide  variation  over  a  small  range  of  Bn. 
Physically,  this  feature  would  be  revealed  by  an 
exceptionally  large  variance  in  the  observed  yield  of  ganuna 
ray  bursts,  a  leading  characteristic  of  this  phenomena  [1-3]. 
Indeed,  a  variance  on  the  order  of  the  average  yield  is 
consistent  with  both  observation  and  the  conditions  given  in 
Eqs.  (54)  and  (55).  (C)  Finally,  sufficiently  low  magnitudes 
of  ka  enable  the  yield  to  approach  the  full  rest  mass  energy 
of  the  original  bod^  provided  B„»Pa,  a  condition  that  is 
obviously  very  well  satisfied  for  steUar  systems. 

We  now  summarize  the  conditions  of  instability 
determined  above.  The  stability  of  a  neutron  star  having 
bose  character  against  disintegration  solely  into 
supersymmetric  pairs  is  governed  by  Eqs.  (54)  and  (55); 
stability  is  maintained  unless  the  mass  of  the  tystem  ^ 
satisfies  an  exact  arithmetic  condition.  Moreover,  the 
precision  involved  is  exceptionally  high,  since  P.  and  B, 
respectively  represent  number  scales  [11,12]  of  ~10*  and 
~10^,  and  a  shift  of  unity  finm  the  condition  specified  is 
governing.  In  addition,  since  the  exact  values  of  R’2nio  and 
the  corresponding  density  of  the  solutions  are  determined  by 
the  detailed  arithmetic  structure  of  the  maw  parameters 
[11],  there  is  a  high  sensitivity  to  the  physical  characteristics 
of  the  daughter  particles.  Finally,  the  key  role  of  the 
monopole  in  relation  to  both  supertymmetric  pairing  and 
baryon  nonconservation  is  expressly  communicated  by  the 
modulus  Pa  in  Eq.  (54). 

The  corresponding  conditions  of  stability  for  the  boson 
decay  channel  given  by  Eq.(44)  and  the  two  fenni 
amplitudes  expressed  Ity  Eqs.(46)  and  (47)  can  be 
understood  mutatis  mutandis.  Ety  inspection  of  Eqs.(44)  and 
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(46),  s>'stems  with  either  hose  or  fenni  character  are  stable 
for 

Bn^Pa+Bz  (57) 

in  which  B^  denotes  the  mass  parameter  of  the  companion 
particle.  Correspondingly,  fermi  species  Bn  is  stable  against 
decay  though  the  amplitude  defined  by  Eq.(47)  for 

Bn  ^  2P(x  +  Bz .  (58) 

From  the  discussion  in  Section  n.A.3.b  above,  the  minimum 
value  of  Bz  is  unity.  It  follows  that  unconditional  stability 
against  all  channels  involving  the  production  of 
supersymmetric  pairs  is  present  for 

Bn^Pa>  (59) 

the  monopole  mass. 


HL  Conclusions 

The  phenomenon  described  above  represents  the 
conversion  of  cold  degenerate  matter  at  the  physical  scale  of 
a  neutron  star  directly  to  energetic  ra^ation  by  a 
fundamental  instability  of  matter  that  is  governed  by  an 
exact  physically  based  counting  mechanism.  Speciii^ly, 
neutron  stars  can  develop  an  instability  toward  spontaneous 
detonation  with  the  production  of  super-relativistic 
supersymmetric  pairs  and  particles  at  specific  mass  values. 
These  conditions  enable  the  violation  of  baiyon  number 
conservation  and  the  resulting  energy  yield  can  approach  the 
rest  mass  energy  of  the  original  body,  an  upper  limit  that 
incorporates  the  observed  output  of  cosmic  gamma-ray 

bursts  including  the  very  large  (~10^’j)  event  (GRB 
990123)  recently  detected  [2,31.  Furthermore,  the 
characteristics  of  the  counting  operation,  which  define  the 
condition  of  instability,  provide  a  simple  explanation  for  the 
large  variance  in  the  observed  energies.  Further,  a  high 
neutrino  yield  is  not  mandatory  in  this  model,  a  feature  that 
both  contrasts  significantly  with  other  descriptions  of 
gamma-ray  bursts  and  is  consistent  with  recent  experimental 
[28]  findings  indicating  a  low  neutrino  emission.  In 
addition,  the  characteristics  of  the  particle  energy  spectrum 
[11]  given  by  Eqs.(9)  and  (10)  can  readily  account  for  the 
production  of  the  very  high  energy  quanta  [~18  GeV] 
reported  [28]  in  cormection  with  GRB  940217.  The 
findings  predia  the  existence  of  a  maximally 
supersymmetry-breaking  particle  pair  consisting  of  a  super¬ 
light  boson  neutrino  vi,  which  plajs  a  fundamental  role  in 
weak  interactions,  and  a  supermassive  form  of  fermionic 
dark  matter  [10].  The  results  also  sun>ort  the  conjecture  that 
the  empiric^y  established  phenomenon  of  lepton  number 
conservation  is  isomorphic  to  the  joint  requirements  of 
angular  momentum  conservation  and  the  supersymmetric 
conservation  of  congruence.  It  is  further  hypothesized  that 


the  existence  of  the  Vi  particle  may  provide  an  alternative 
theoretical  basis  for  (1)  CP  violation  in  decay  [8,9]  and 

(2)  the  attribute  of  strangeness,  a  characteristic  example  of 
which  is  given  by  the  properties  of  hyperon  decay  into 
nonleptonic  channels.  C>vet^l,  the  pattern  of  findings  that 
has  emerged  points  strongly  to  the  existence  of  a  highly 
energetic  realm  of  physical  interactions  that  is  based  on 
supersymmetry  and  a  new  arithmetic  organizing  principle. 
Its  presence  imposes  a  form  of  scale  invariance  that  enables 
sm^  particles  and  large  bodies  to  have  similar 
characteristics  of  stability.  Hence,  both  neutrons  and 
neutron  stars  can  exhibit  unstable  behavior. 
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A  zone  of  stability  given  by  £q.  (56)  is  indicated.  Unconditional  stability  exists  against  all 
decay  channels  involving  supersyxnmetric  pairs  for  both  fertni  and  bose  character  systems  for 
Bn  i  Po(  See  text  for  discussion  and  definitions  of  the  corresponding  parameters. 
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ABSTRACT 

The  construction  of  inverse  states  in  a  finite  field  enables  the  organization  of  the  mass  scale  with 
fi.nHamAntQi  octcts  in  sn  eigjit-dimensional  index  space‘that  identifies  particle  states  with  residue  class 
designations.  Conformance  with  both  CPT  invariance  and  the  concept  of  supersymmetry  follows  as  a 
direct  consequence  of  this  formulation.  Based  on  two  parameters  (P„  and  gj  that  are  anchored  on  a 
concordance  of  physical  data,  this  treatment  leads  to  (1)  a  prospective  mass  for  the  muon  ne^o  of 
~  27.68  meV,  (2)  a  value  of  the  unified  strong-electiowcak  coupling  constant  a*  =  (34.26)"  that  is 
physically  by  the  ratio  of  the  electron  neutrino  and  muon  neutrino  masses,  and  (3)  a  seesaw 

congruence  ^nnArring  the  ffiggs,  the  electron  neutrino,  and  the  union  neutrino  masses.  Specific 
evaluation  of  the  masses  of  the  corresponding  supersymmetric  Higgs  pair  reveals  that  both  particles  arc 
superheavy  (>10“  GeV).  No  renormalization  of  the  Higp  masses  is  introduced,  since  the  calculational 
procedure  yielding  their  magnitudes  is  intrinsically  divergence-free.  Fur^,  the  Higgs  fulfills  its 
conjectured  role  through  the  seesaw  relation  as  the  particle  defining  the  origin  of  all  particle  masses, 
since  the  electron  and  muon  neutrino  systems,  together  with  their  siqwrsymmemc  partners,  are  the 
gfnAratnre  of  the  maw  scalc  and  establish  the  corresponding  index  space.  Finally,  since  the  computauon 
of  the  wiggg  "<«««■«  is  entirely  determined  by  the  modulus  of  the  field  P^,  which  is  fully  defined  by  the 
large-scale  parameters  of  the  universe  through  the  value  of  the  universal  gravitational  constant  G  and  the 
requirement  for  perfect  flatness  (Q  =  1.0),  the  seesaw  congraence  fuses  the  concepts  of  mass  and  space 
and  creates  a  new  unified  ardietype. 
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L  Introduction 

The  concept  of  mass  [1]  has  been  an  everlasting  challenge  to  physi^  theory  [2].  In  general,  an 
organization  of  the  particle  mass  scale  is  sought  that  exhibits  certain  salient  features.  Four  prominent 

desired  characteristics  are  (A)  the  equivalence  of  particle  P  and  antiparticle  P  masses  (mp  =*  mp),  a 

demand  of  CPT  invariance  [3,4],  (B)  an  expression  of  supersymmetry  [5],  a  principle  invoking  fermion 
/boson  (b)  pairing  that  presently  exists  as  a  very  attractive  theoretical  hjqiothesis,  (C)  the  existence  of 
the  Higgs  particle  [6],  the  entity  which  introduces  mass  into  the  Standard  Model  [7],  and  (D)  a  basis  for 
the  value  of  the  unified  strong-electroweak  coupling  constant  a*,  the  physical  parameter  regulating  non- 
gravitational  interactions  [8]. 

Taken  together,  the  two  features  (A)  and  (B)  of  the  mass  spectrum  noted  above  incorporate  three 
forms  of  two-particle  associations.  The  requirement  of  CPT  separately  pairs  both  fermi  (e.g.  e  /e~)  and 
bose  (e.g.  tCItC  and  7i®/ic*)  species.  Supersymmetry  connects  the  fermi  and  bose  particle  genres  through 
the  formation  of  corresponding  supersymmetric  pairs  (?/?*.)•  Therefore,  these  two  classes  of 
relationships  produce  three  pairings  b/b,  and  b/Q  yielding  a  system  in  which  each  particle  P  has 
both  an  antiparticle  P  and  a  supersymmetric  partner  Psi.  The  overall  result  is  the  construction  of  the 

particle  mass  scale  [9]  with  fundamental  quartets  (P,  P,  P»,  P«).  Two  previous  studies  [9,10]  explicitly 
illustrate  the  manner  in  which  this  mass  scale  can  be  organized.  A  leading  prediction  of  these  analyses  is 
a  prospective  value  of  0.808  meV  for  the  mass  mve  of  the  electron  neutrino  v*.  a  magnitude  that  is 
consistent  with  currently  avmlable  data  [1 1]. . 

The  present  work  demonstrates  that  the  three  classes  of  pairings  associated  with  CPT  and 
supersymmetry  given  above  follow  as  immediate  consequences  of  the  defimtion  of  a  new  class  of 
particle  states  designated  as  inverse  states.  This  description  additionally  associates  each  particle  state  P 
with  a  corresponding  inverse  state  Pm  through  a  specific  relationship  of  the  particle  masses.  The 
outcome  is  an  affiliation  of  four  particles  (P,  Pm,  Pin,  (Pin)«),  normally,  but  not  necessarily,  ^h 
nondegenerate  mass  values,  which  are  further  related  to  a  corresponding  quartet  of  antiparticles.  The 
final  result  is  (1)  the  formation  of  an  octet  of  states  that  generally  represents  four  distinct  mass  values 
and  (2)  the  construction  of  a  corresponding  eight-dimensional  index  space  that  provides  a  system  of 
particle  classification. 

This  theoretical  picture  can  be  implemented  by  incorporating  the  previous  analyses  [9,10]  into  the 
mathematical  structure  of  a  finite  field  [12].  Since  an  algebraic  field  possesses  the  operation  of  division, 
this  extension  permits  the  definition  of  inverses  and  the  corresponding  construction  of  a  set  of  umts 
{U}.  An  important  consequence  is  the  fiict  that  these  units  comprise  a  group  structure  [12].  As  shown  in 
the  development  below,  we  will  be  primarily  concerned  with  the  subset  of  fields  fy  in  which 

P  3  l(mod  4),  the  physically  motivated  case  [9,10,13]. 

The  ability  to  define  inverse  states  leads  to  several  new  statements  concerning  the  organization  of 
the  mass  scale  and  satisfaction  of  the  four  desired  characteristics  (A-D)  listed  above.  Chief  among  the 
specific  findings  are  (1)  the  muon  neutrino  mass  with  a  predicted  value  of 

m  =27.68  me (1) 

(2)  the  conclusion  that  the  unified  strong-electroweak  coupling  constant  a*  is  physically  defined 
through  the  electron/muon  neutrino  mass  ratio  by 
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(2) 


niv  , 

a*=  —S-=(34.26)-\ 

niy 

»» 

(3)  the  identification  of  uniquely  specified  candidate  values  for  the  masses  of  a  supersymmetric  doublet 
that  corresponds  to  the  Higgs  particle,  and  (4)  a  seesaw  congruence  that  connects  the  electron  neutrino 
mass,  the  muon  neutrino  mass,  and  the  masses  of  the  supersymmetnc  Higgs  pair  [6-8,14].  Explicit 
evaluation  of  the  pair  of  Higgs  masses  indicates  that  both  particles  are  in  the  wperheavy  regime 
(>10^*  GeV).  Since  the  computation  of  the  Higgs  masses  is  intrinsically  divergence-fi-ee,  no 
renormalization  [6]  is  necessary.  The  stated  values  of  this  set  of  fundamental  parameters  are  consistent 
with  the  ranges  of  their  potentially  expected  magnitudes  [6-8,11,14-17]  and  all  currently  available 
experimental  data.  The  overarching  conclusion  that  follows  firom  the  form  of  the  seesaw  relation  derived 
in  this  analysis  is  that  the  concepts  of  space  [18.191  and  mass  [1.2]  are  fundamentally  inseparable  The 
r«.giiits  directly  that  these  two  phvsical  archetypes  constitute  an  identity. 


n.  Theoretical  Development 

A.  Elementary  Properties  of  Finite  Fields 

A  commutative  division  ring  is  designated  as  a  field  [12].  For  any  prime  integer  P,  the  field 
consists  of  the  set  of  residue  classes 


{[0]p,[l]f,[2]p,...,[P-l]p}.  (3) 

Computations  performed  with  elements  [x]p  of  Fpare  conducted  with  the  customary  rules  of  arithmetic 
and  reduction  modulo  the  prime  P.  Since  all  primes  possess  primitive  roots  [20],  power  residue  systems 
of  enumeration  in  FpCan  be  established  and  the  properties  of  congruences  [21]  are  central.  Physically 

anchored  procedures  of  counting  based  on  these  concepts  were  used  in  earlier  studies  to  describe  the 
organization  of  the  particle  mass  scale  [9,10]  and  the  phenomenon  ofy-ray  bursts  [22]. 

The  existence  of  the  operation  of  division  enables  the  definition  of  an  inverse  [x]p  for  every 
element  [x]p  of  Fp  except  [0]p.  Hence,  a  set  of  units  {U},  defined  by  the  statement 

MpMp^  =  [l]p,  (4) 

exists  and  it  is  known  that  these  units  (U)  form  a  multiplicative  group  [12].  Furthermore,  a  simple 
general  procedure  exists  for  the  calculation  of  the  inverses  [x]p‘  based  on  Fermat’s  Little  Theorem  [23], 
Specifically,  Fermat’s  result  states  that 

V  X,  gcd(x,  P)  =  1,  [x]J  =  [l]p,  (5) 

fi-om  which  it  follows  directly  that 

[x]p-‘=[x]J-*.  (6) 
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Eq.(6)  explicitly  shows  that  an  inverse  [x]'*  generally  depends  on  both  the  value  of  x  and  the 
modulus  P. 


B.  Specific  Mathematical  Properties  of  Inverses 

All  inverses  can  be  computed  directly  through  Eq.(6).  However,  a  restricted  subset  of  inverses  has  a 
special  elementary  formulation.  Anticipating  the  need  for  this  result  below,  we  state  the  lemma 
specifying  this  set  of  residues  and  give  the  corresponding  proof 

Lemma  (I):  Let  Pa  be  prime,  B*  be  a  divisor  of  Pa  -  1,  and  (Bx)«  =  Pa  -  . 

Then,  [(Bx)ss]p  =  [(Pa  ~  1)/Bx]p 


Proof:  By  definition,  [(Bx)«]  p^=  [(Bx)„]  p 


,P.-2 


(7) 

(8) 
(9) 


P  —  2 

Since  (Pa  -  Bx/‘"  ^  =  Z  (-1)‘(^“7 ^  )pS“" ^  B* 

i  =  0  * 

P--  2 

~  “Bx  *  +  k  •  Pa  for  some  integer  k, 

we  have  (Bx)S“~  -bJ““  ^(mod  Pa)  . 

P  -  2 

Suppose  Bx*  s  A(modPa), 
then  -Bx*~  p.  -  A(mod  Pa)  . 

In  order  to  evaluate  the  constant  A  fi-om  Eq.  (8),  we  write 

^  =  Bx^®  ^  B*  s  ABx(mod  Pa),  if gcd  (B* ,  Pa)  =  1. 

Therefore,  ABx  =  1  +  k*  •  Pa  =  (kx+  1)  +  kx(Pa  -  1)  for  some  integer  k*. 

Then  A  =  +  k*  .  Since  —  is  an  integer  by  assumption,  ^  must  also 

*  X  °x  Bx 

be  an  integer. 

Let  — — — = N ,  where  N  is  an  integer. 


Then  we  have  A  =  N  +  (NBx  -  1) 


=  N  +  N(Pa-l)- 


Pg-l 

Bx 

Pg-l 

B. 


=  NPa- 


Pg-  1 
B„ 
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From  Eq.(9), 


-B5‘"^  =  Pa-NPa  +  ^5^  s  Za_i(modPa), 

®x  ®x 

and  from  Eq.(7) 


[(Bx)ss5‘  =  ^  =  [(Pc  -  1)/B  . 


(10) 


Q.E.D. 


Hence,  B*  must  be  a  divisor  of  Pa  -1. 

Lemma  (1)  specifies  a  subset  of  residue  classes  whose  corresponding  supersymmetric  [9,10] 
partners  (B^^  possess  inverse  states  of  the  special  form  given  by  Eq.(lO).  The  number  of  elements  of 
the  set  {Bjp}  is  given  by  the  arithmetic  function  d(Pa  -  1),  the  function  [24]  which  represents  the  number 
of  divisors  of  the  integer  (Pa  -  1).  On  the  basis  of  an  earlier  study  [13],  d(Pa  -  1)  =  2.32  x  10‘*.  By 
inspection  of  Eq.(lO),  all  members  of  {B,j}  obey  the  congruence 

pi  KB,).];'- -l(modP„).  (11) 

a  o 

The  subset  {g}  of  elements  of  whose  members  are  primitive  roots  of  the  modulus  Pa  also  enjoy  a 
select  status  [9,10,21,22],  since  these  integers  serve  a  generating  function.  Of  particular  significance  will 
be  the  restricted  subset  {g,}  of  primitive  roots  such  that 

{g,}  =  {g}  n  {BJ.  (12) 


Specifically,  it  will  be  shown  below  that  the  subset  {g^}  relates  directly  to  both  the  magnitude  of  the 

unified  strong-electroweak  coupling  constant  a*  and  the  values  for  the  masses  of  the  electron  and  muon 
neutrinos.  We  note  further  that  since  {g,^}  can  be  empty,  as  the  prime  41  demonstrates  [21],  the  demand 
for  the  existence  of  a  nonempty  set  {g^}  constitutes  an  additional  mathematical  requirement  for 
acceptable  physical  values  [  13]  of  the  prime  P„ . 

In  order  to  explore  the  needed  characteristics  of  primitive  root  inverse  mappings,  we  will  initially 
limit  our  treatment  to  moduli  of  the  form  P  s  l(mod  4)  and  pairs  of  primitive  roots  (g^,  gp)  that  sum  to 
the  modulus  P„ ,  the  condition  that  expresses  supersymmetry  [9,10].  Further,  since  the  Euler  function 
<KP  -  1)  is  always  even  [24-26]  for  any  prime  P  >  3,  and  a  result  of  Gauss  [27]  certifies  that  all 
primitive  roots  of  primes  of  the  genre  P„  h  l(mod  4)  can  be  arranged  in  supersymmetric  pairs  [9,10] 
such  that 


ga  +  gp  =  P«. 


(13) 
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it  follows  from  Lemma  (2)  below  that  the  inverse  state  operation  given  by  Eq.(6)  always  maps  a 
primitive  root  pair  gp)  according  to  the  pattern 


(ga.gp)^(ga>  gp) 

(14) 

with 

&x  ^  8a 

(15) 

and 

go  +  gp=go+gp=  Pa- 

(16) 

Therefore,  primitive  root  pairs  perforce  transform  to  different  primitive  root  pairs  with  both  pairs  in 
conformance  with  Eq.(13),  the  postulated  mass  relationship  expressing  supersymmetry  [9,10,22],  This  is 
an  extension  of  the  Q-pairing  pattern  described  in  an  earlier  study  [9]  of  the  mass  scale.  We  now  prove 
this  important  statement. 

Lemma  (2):  Assume  go  +  gp  =  Pa ,  ga ,  gp  are  primitive  roots  of  Pa ,  and  Pa  prime. 

Let  the  inverses  of  go  and  gp  respectively  be  Ba  and  Bp. 

Then,  (1)  Bo  +  Bp  =  Pa. 

(2)  Bo  and  Bp  are  primitive  roots  of  Pa . 

Proof:  (1)  can  be  proved  by  using  Lemma  (1)  of  an  earlier  study  [9]  whose  proof  was  therein 
given. 

Let  x  =  Pa  -  2,  an  odd  number. 

By  the  previously  demonstrated  lemma  [9],  we  have 

and 

Ba  Bp  ~  Po  • 

To  show  that  both  Ba  and  Bp  are  primitive  roots  of  Po,  we  use  the  relation 

go  •  Ba  s  l(mod  Pa), 

the  defining  statement  of  the  inverse  mapping. 

For  any  y  ^  Po- 1, 

g2-B^  =l(modPa) 

^  ^  ^  (mod  P«) 

^“"y-^modPa) 

o  bIL  s  ~^(mod  Po),  since  ^  s  l(mod  Po). 

P  -Y— 1 

The  residue  of  Bo  is  in  one-to-one  relation  with  the  residue  of  go*  .  Hence,  Bo  is  a  primitive 
root.  The  corresponding  result  for  Bp  follows  mutatis  mutandis. 


QED 


It  follows  from  Lemmas  (1)  and  (2),  if  either  gaor  gp  is  a  member  of  {g^},  the  (ga,  gp)-)>(g'a,  g'p)  inverse 
transformation  is  immediately  calculable  on  the  basis  of  the  divisor  structure  of  -  1  through  Eq.(lO). 

C.  Mass  Parameter  Inverse  Mapping  Characteristics 

Scale  invariance  of  the  representation  of  the  mass  parameters  is  a  fundamental  characteristic  of  the 
power  residue  formulation  [9,10].  Therefore,  the  basic  features  of  the  inverse  state  mapping  described 
above  can  be  folly  illustrated  with  an  example  utilizing  a  prime  modulus  of  small  magnitude.  Table  I 
displays  the  inverse  state  primitive  root  pairings  for  the  prime  modulus  Pa  =  37.  We  observe  that  one 

primitive  root  (ga  =  2)  is  a  divisor  of  Pa  -  1  =  36,  so  that  [gpi^j  =18,  the  result  given  directly  by 

Eq.(lO).  In  order  to  show  the  gross  characteristics  of  the  scaling  of  the  pairing  pattern  on  the  magnitude 
of  the  modulus,  a  second  example  of  the  primitive  root  inverse  mapping  is  given  in  Table  n  for  P„=  73 
and  it  is  seen  that  the  qualitative  features  of  the  transformation  are  unaltered. 

Two  new  conserved  quantities  are  revealed  by  the  primitive  root  inverse  mappings  illustrated  in 
Tables  I  and  H.  The  statements  of  the  conserved  magnitudes  are  presented  in  Table  HI.  It  is  conjectured 
that  these  relationships  correspond  to  the  conservation  of  internal  attributes  of  the  particle  under  the 
inverse  state  operation.  Further  work  [28]  is  underway  to  elucidate  the  potential  physical  interpretations 
of  these  conserved  quantities  in  terms  of  the  subgroup  structure  of  the  units  of  Fp  . 

*  a 

An  explicit  picture  of  the  quartet  of  residues  representing  the  octet  of  physical  particle  masses 
{P,  P,  P«,  Pm,  Pufc  Pin,  (PiD)n,  (Piii)ss}  Can  be  niade  by  further  examination  of  the  example  for  Pa  =  37.  For 
this  illustration,  we  will  take  (g^,  gp)  =  (2,  35),  one  of  the  primitive  root  pairs  identified  in  Table  I.  Fig. 
(1)  shows  the  spectrum  of  these  primitive  roots  in  conjunction  with  the  corresponding  inverse  state  pairs 
(g^,  gp)  =  (19,  18).  These  four  generators  can  now  be  used  jointly  to  organize  both  (a)  the  particle 

mass  parameters  and  (b)  the  corresponding  index  space  for  state  classification. 

Inspection  of  Fig.  (1)  shows  that  the  operators  corresponding  to  the  supersymmetry  O*  and  inverse 
state  Oin  transformations  commute,  so  that 


0i„0„  =  0„0i„. 


(17) 


We  also  have 


0?=Oi=l,  (18) 

the  unity  operator,  and  the  inverses 

(Oi.]-‘=Oi.,  [0.r‘=0..  (19) 

so  that  {Oin,  Om}  constitutes  an  abelian  group. 

Table  IV  comprehensively  illustrates  the  mappings  of  the  residues  that  arise  from  the  four  primitive 
roots  shown  in  Fig.  (1).  The  particle  states,  supersymmetric  partners,  and  corresponding  inverse  states 
that  occur  through  organization  by  the  quartet  of  primitive  roots  are  shown.  In  general,  each  residue 
possesses  four  power  indices  (h;  i  =1,2,3 ,4}  and  four  corresponding  modulus  indices  {Iq,  i  *1, 2,3,4}, 
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parameters  which  define  the  coordinates  of  the  8-dimensional  space  in  terms  of  residue  classes 
[x]p  designations. 

Examination  of  Table  IV  demonstrates  that  the  Q-  and  D-pairing  motifs  developed  in  an  earlier 
study  [9]  remain.  The  extension  of  the  pattern  with  inverse  states  simply  adds  new  relationships  that  the 
four-dimensional  structure  could  not  reveal.  Specifically,  for  power  index  t  odd,  a  supersymmetric  pair 
is  coupled  to  a  corresponding  inverse  supersymmetric  pair  (e.g.  t  =  7,  (20,  17)  (13,  24))  that 

preserves  the  Q-pairing  structure.  All  supersymmetric  primitive  root  pairings  faU  in  this  class,  as 
described  in  Section  n.B.  Likewise,  for  power  index  t  even,  the  degenerate  pattern  exhibited  by  D- 
pairing  [9]  is  transferred  to  the  corresponding  inverse  state  pair  which  also  exhibits  D-pairing.  The 
appropriate  conjugate  panels  [9]  are  always  related  by  a  shift  in  the  power  index  t  by  the  transformation 
t  t  +  (Po  -  l)/2,  the  relation  established  in  the  earlier  work  [9]. 
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INVERSE  STATE  PROPERTIES 

Power  Residue  Systematics  (Pa  = 
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Table  I:  The  systematics  of  the  power  residue  counting  system  represented  by  (ga,  gp.  Pa)  =  (2,  35, 
37).  The  parameter  t  represents  the  respective  integral  power  indices  ta  and  tp.  In  green  are  the 
residues  (Bq)  for  ga  =  2  and  correspondingly  (Bp)  in  red  for  gp  =  35.  Two  patterns  are  apparent;  the 
one  for  t  odd  always  yields  Ba  +  Bp  =  Pa  =  37,  the  value  coincident  with  (ga  +  gp  =  Pa  =37,  the 
condition  governing  supers>'mmetr\'.  The  remaining  entries  for  t  even  always  represent  the  degenerate 
condition  Ba  =  Bp.  Conjugate  panels  are  shifted  t  t  +  (Pa-l)/2  uniformly  in  both  cases  such  that  the 
respective  sums  give  Pa.  Hence,  both  patterns  provide  pairing  which  represents  the  prime  modulus 
P„  =  37.  Every  residue  has  an  associated  inverse  in  a  field.  Shown  below  is  the  subclass 
corresponding  to  the  primitive  root  pairings  for  prime  Pa  =  37  established  by  constructing  units  in  the 
finite  field  Fs?-  The  primitive  roots  are  {2,  5, 13, 15, 17,  18,  19, 20, 22, 24, 32, 35}. 


STATES 

INVERSE  STATES 

(ga»  ga) 

(&;»  ge) 

(2, 35)  - 

- 

(19, 18) 

(5. 32)  - 

- ^ 

(15,22) 

(13,24)  - 

- 

(20, 17) 

The  inverse  operation  in  F37  relates  conjugate  panels  for  t^  distinct  sets  primitive  roots,  thereby 
conferring  a  unique  status  on  these  four  residues.  For  example,  the  inverse  mapping  (1)  satisfies 
supersymmetry  (ga  +gp  =  ga'  +  gp'  =  Pa  ==  37)  and  (2)  the  pair  (ga,  gp)  =  (2,  35)  is  transformed  to 
(8^1+  8^)  (19,  18),  another  primitive  root  pair.  Therefore,  t^  distinct  quartet  manifolds  are  related  by 

the  inverse  mapping  for  primitive  root  states.  Hence,  four  (2,18,19,35)  distinct  masses  are  related 
together  with  a  corresponding  quartet  of  antiparticle  states.  Therefore,  the  existence  of  inverse  states 
give  the  mass  scale  an  extended  foundation  through  the  construction  of  octets  and  a  corresponding 
eight-dimensional  index  space  which  classifies  the  particle  states  with  residue  class  designations.  These 
octets  involve  (A)  particle/antiparticle  (CPT),  (B)  particle/supersymmetric  particle  (Supersymmetry), 
and  (C)  particle/inverse  particle  (Unitary  symmetry)  relations.  Generally,  [ga][gpl~’  =  -l(mod  Pa)  and 
only  t  odd  states  are  involved  in  the  primitive  root  mappings. 


Residue  Mappings  for  Supersymmetric  Primitive  Root  Pairs  (P^  =  37) 
l8aJ37~2,  [gp]37~  35,  =19,  [gpj^7  =18 
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25 

17 

20 

24 

13 

26 

3 

3 

25 

25 

27 

31 

6 

6 

31 

28 

12 

12 

34 

34 

29 

13 

24 

20 

17 

30 

11 

11 

27 

27 

31 

15 

22 

5 

32 

32 

7 

7 

16 

16 

33 

23 

14 

29 

8 

34 

28 

28 

4 

4 

35 

18 

19 

35 

2 

36 

1 

1 

1 

1 

Table  IV:  Pattern  of  residue  classes  generated  by  the  quartet  of  inverse  state-related  primitive  roots 

inverse  ^ 

(got>  gp)  “  (2,35)  ^  (g^,  g^)  =  ( 1 9, 1 8)  involving  the  residues  B,  their  supersymmetric  partners  Bss, 
and  the  corresponding  inverse  state  pairs  Bin  and  (Bi„)ss.  See  Fig.  (1)  for  reference. 


Power  Residue  Systematics  (Pa=73) 

ga=34,  gp=39,  ga  +  gp=Pa 


Primitive  root  pair  its  inverse  power  index  power  index  mapping  sum 


{5,68}  -^{44,29} 

_  5  ^  67 

72 

{11,62}  ->  {20,53} 

_  59  ->  13 

72 

{13,60}  ^  {45,28} 

_  7  ■>  65 

72 

{14,59} ->{47,26} 

_  61  ->  11 

72 

{15,58}  ^{39,34} 

_  35  ->  37 

72 

{31,42}  ^  {33,40} 

_  55  ->  17 

72 

Table  EE:  Systematics  of  the  inverse  states  for  the  prime  Pa  =  73  represented  for  the  primitive  roots 
ga  =  34  and  gp  =  39. 


NEW  STATES  and  NEW  CONSERVATION  LAWS 

•  Inverse  States  Define  a  New  Symmetry  Based  on 
Supersymmetry  Preserving  Primitive  Root  Mapping 

•  (ga,  gp) - ►  (gi,  gp 

ga  gp  ~  ?  ga  gp  “  ?  ga  ^  ga 

•  Inverse  State  Transformations  Define  New  Conservation 

Laws 

•  Power  Index  t 

(ga ,  gp)  - ►  (gi  ,  g'p) 


•  Modulus  Index  k 

kaga  +  Kga  =  0  (mod  Pa) 
kpgp  +  kpgp  =  0  (mod  Pa) 

Presumably  indicate  conservation  of  internal  attributes 

•  One  New  Fundamental  State  Classification  (Symmetry) 

and 

Two  New  Conserved  Quantities 

Table  EH:  Inverse  states  define  a  new  symmetry  based  on  a  supersymmetry  preserving  mapping  that 
can  be  achieved  with  =  l(mod  4). 


Pa  =37 


Pa=  l(mod4) 


[Spls? 


[Sala?  ^ 


=fci37 

18  =  [g,]«  =  [gi,)37 


Fig.  ( 1 ):  Example  of  the  spectrum  of  primitive  roots  governing  the  octet  state  manifolds  for  =  37. 

The  primitive  roots  are  (g^,  g^)  =  (2,  35);  the  corresponding  supersymmetric  inverses  of 
(ga.  8(5)are(g;,gp  =  (19, 18). 


nL  Computational  Findings 
A.  Muon  Neutrino  Mass 

The  results  derived  in  Section  n  above  can  be  applied  to  several  specific  computations.  The  use  of 
Eq.(lO)  petmits  the  determination  of  a  mass  parameter  which  we  prospectively  identify  as  the  muon 
neutrino  niass  m^  .  An  earli^  study  [9]  gave  an  electron  neutrino  mass  m^^  of  ~  0.808  meV,  a  value 

corresponding  to  the  mass  parameter  [9,10] 

g«  =  2q  =  2-3^-7-13-31-37-47-53-59-71-73.83*107-109-113*137-139-149.  (20) 

Direct  evaluation  [9, 10,13]  showed  that  ga  is  a  primitive  root  of  the  modulus  Pa,  the  ritd  property  that 
lets  it  serve  as  a  generator.  Furthermore,  since  ga  is  manifestly  [13]  a  member  of  the  divisor  set  {Bx}  of 
the  integer 


Pa-1  =  2’-3*5*^7-11-13-17-19-23-29-31-374143-47*53-59-61-67-71-73 
.79-83-89-97  10M03  107  109113  127-13  M37.139- 149151, 

Eq.(lO)  of  Lemma  (1)  can  be  used  to  compute  directly  the  inverse  of 

8p  ~  ~  8a  • 


(21) 


(22) 


This  gives 


2-5^1M7-19-23-2941  43-61-67-79-89-97-10M03-127-13  M51 


(23) 


in  which  the  integer  r  is  defined  in  an  earlier  study  [13].  Accordingly,  on  the  basis  of  the  previous  work 
[9,10,13],  we  immediately  recognize  the  specific  magnitude  of  [gp^J  as  a  &ctor  in  the  known  result 

[13] 


qr  =  (2qXr/2)  =  Pa-l. 


(24) 


These  earlier  analyses  [9,10,13]  also  demonstrated  the  important  finding  that  the  fine-structure  constant 
a  could  be  uniquely  represented  as 


a  =  ^  =  (137.03599926808)"^  (25) 

r 

and  that  this  value  is  in  full  agreement  with  the  measured  magnitude  [29,30]. 

In  summary,  Eqs.  (22)  and  (24)  can  be  respectively  restated  as 

(g„lp  +[g|,]p  -0(n>odPa).  (26) 

a  a 
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the  additive  relation  representing  supersymmetry  [9,10]  in  the  field  Fp,  and 

tea]  =  -KmodP  ),  (27) 

Pa  ^  ® 

a  corresponding  multiplicative  condition.  However,  the  unusual  factor  structure  exhibited  by  [gp;^!  in 

Eq.(23).  specifically,  that  [g0t^6{Bx},  requires  the  special  condition  that  [g  ]„  e{B  }  hold,  an 

outcome  legislated  by  Eq.(lO).  Only  a  very  rare  fi-action  of  the  primitive  roots  of  Pa  will  satisfy  this 
exceptional  requirement  and,  without  conformance  to  this  condition,  the  crucial  statements  given  in 
Eqs.(24)  and  (25)  would  not  follow.  For  the  physically  relevant  values  [9,10,13]  of  Pa,  this  firaction  is 
less  than  one  in  ~  1 0^^. 

The  magnitude,  smoothness,  and  factor  structure  of  [gp^l  correspond  to  a  light  propagating  particle 

[13]  whose  properties  would  resemble  those  of  the  electron  neutrino  which  has  been  prospectively 
identified  [9,10]  with  a  mass  parameter  of  ga  and  a  physical  mass  of -0.808  meV.  From  the  symmetry 
expressed  by  Eq.(27),  we  will  provisionally  assign  [gp]^l  as  the  mass  number  of  the  muon  neutrino 

This  gives  the  physical  mass  corresponding  to  the  muon  neutrino  v^as 

m^^  =  27.68  meV ,  (28) 

a  value  that  is  consistent  with  both  the  available  data  and  the  range  expected  for  this  particle  [1 1]. 

B.  Unified  Strong-Electroweak  Coupling  Constant  a* 

The  results  presented  in  Section  in.A  above  lead  directly  to  an  expression  for  the  magnitude  of  the 
unified  strong-electroweak  coupling  constant  a*.  A  previous  study  [9]  derived  the  relationship 

(a*)-^  =  34.26  (29) 

8a 

on  the  basis  of  characteristics  of  the  mass  scale.  However,  the  introduction  of  the  concept  of  inverse 
states,  the  results  given  in  Eqs.(20),  (22),  (23),  and  (26),  and  the  identification  of  m^^  with  the  quantity 

[gp]p^ ,  enable  Eq.(29)  to  be  immediately  recast  in  the  form  of  the  fundamental  mass  ratio 


=  (34.26)- 1  . 


(30) 


The  expression  given  in  Eq.(30)  has  several  attractive  properties.  First,  it  gives  a  magnitude  that 
stands  in  good  correspondence  with  the  anticipated  range  [7,8,16,17]  of  a*.  Further,  it  simultaneously 
relates  the  utterly  basic  quantity  a*  mathematically  to  the  ratio  of  two  primitive  roots  that  are  members 
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of  a  restricted  set  {g*}  of  special  elements  in  the  field  and  physically  to  the  mass  ratio  of  two 
fundamental  particles.  Moreover,  the  integers  [g  ]  and  their  corresponding 

Pa  ^  Pa 


respective  supersymmetric  partners  [g-]  and  [g„]~^ ,  are  all  primitive  roots  by  Lemma  2.  This  enables 


Pa 


Pa 


them  to  serve  collectively  as  the  generators  of  the  mass  scale  and  the  organizers  of  the  eight¬ 
dimensional  index  space  t^t  classifies  all  particle  states.  Fig.  (2),  in  parallel  with  the  model  shown  in 
Fig.  (1),  illustrates  the  rdationships  connecting  these  four  primitive  roots.  Basic  physical  entities  and 
correspondingly  fundamental  mathematical  quantities  unite  to  (a)  construct  the  particle  mass  scale  and 
(b)  specify  the  magnitude  of  the  cardinal  physical  parameter  governing  their  interactions. 
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MASS  SCALE  OCTET  GENERATORS 

Pa  =  l(mod  4) 


(^e)ss»  [Sr^P  ^ 


Fig.  (2);  Illustration  of  the  primitive  root  states  that  organize  the  mass  scale,  construct  the 
eight-dimensional  index  space,  and  specify  the  unified  strong-electroweak  coupling  constant 
a*.  The  magnitude  of  Pa  represents  the  monopole  mass  niQK>,  Eo  denotes  the  energy  unit,  and 
vi  a  boson  neutrino  with  a  mass  equal  to  Eo.  See  references  [9,10,13,22]  for  further  discussion. 


C.  Definition  and  Computation  of  the  Higgs  Mass 

The  analysis  presented  above  demonstrates  how  the  definition  of  inverse  states  extends  the 
previously  developed  four-dimensional  index  space  [9]  to  eight  dimensions,  thweby  enriching  its 
classificational  capacity.  Inspection  of  the  pattern  represented  by  the  inverse  states  immediately 
identifies  a  fundamental  degeneracy.  For  a  particle  P,  it  is  the  condition  of  the  equality  of  the  mass 
parameters  for  the  supersymmetric  and  inverse  Pin  states;  hence, 


B„  =  B 


iii» 


(31) 


a  statement  signaling  a  point  of  symmetry.  Furthermore,  this  degeneracy  only  occurs  for  two  critical 
values  of  the  power  index 


(32) 


From  Eq.(32)  it  is  apparent  that  the  existence  of  the  degeneracy  given  by  Eq.(31)  requires 
Po  =  l(mod  4),  a  condition  made  necessary  earlier  by  several  other  considwations  [9,10,13], 
Accordingly,  this  degeneracy  appears  in  Table  IV  for  t  =  9  and  t  =  27.  In  these  two  cases,  the  residues 
represent  the  equivalent  supersymmetric  pairs  (6,  31)  and  (31,  6).  It  can  be  fijrther  shovm  that  the 
magnitudes  of  these  residues  are  unique,  since  they  are  invariant  over  the  full  set  of  primitive  roots  of 
the  prime  Pa.  It  follows  that  these  two  residues  (e  g.  6  and  31)  depend  solely  on  the  value  of  the  prime 
modulus  Pa,  the  large  scale  parameter  in  the  counting  system,  and  obey  the  quadratic  congruence 

x^  =  -l(modPa).  (33) 


Indeed,  it  can  be  shown  [23]  that  Eq.(33)  has  a  solution  if  and  only  if  Pa  =  l(mod  4),  the  condition 
stipulated  by  Eq.(32).  Fintdly,  these  two  residues  always  have  order  4,  hence,  the  quaitic  congruence 

x^sl(modPa)  (34) 


simultaneously  holds. 

The  unique  status  of  these  residues  that  is  established  by  Eqs.(3 1),  (32),  and  (33),  in  alliance  with  the 
fact  that  they  depend  only  on  the  modulus  P®,  whose  value  is  fully  determined  by  the  large  scale 
parameters  of  the  universe  [13],  nominates  these  supersymmetric  states  as  candidates  for  the  Higgs 
particle.  With  this  assignment,  we  can  write  fi-om  Eq.(33)  the  mass  parameter  of  the  Higgs  system  as 

B’»„=-Kn.odP„).  (35) 


We  now  explore  the  consequences  of  this  prospective  identification. 

Consider  the  definition  of  the  Higgs  mass  presented  in  Fig.  (3).  The  defining  statement  given  by 
Eq.(31)  is  expressed  by  the  simultaneous  fulfillment  of  two  conditions,  one  additive  and  the  other 
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multiplicative.  These  conditions  respectively  relate  (1)  the  supersymmetric  (X,Y)  pair  mass  distribution 
and  (2)  the  inverse  state  definition.  They  are 


[X)p__+tYi,^=0(niodP„) 

(36) 

[X]p^  [Y5*  •  -HmodPa)- 

(37) 

Satisfaction  of  Eqs.(36)  and  (37)  leads  to  the  statement  given  in  Eq.(31),  the  degeneracy  of  the  P„  and 
Pin  states  and  the  specification  of  the  supersymmetric  residue  pair  uniquely  dependent  upon  the 

modulus  Pa- 


D.  Seesaw  Mass  Congruence 

The  results  expressed  in  Eqs.(27^  and  (35)  immediately  lead  to  the  relation 

“  [®Higgs^^  “  0(mod  Pqj) 

and  the  seesaw  statement 

fealp^  ^  ’ 


(38) 

(39) 
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HIGGS  MASS  DEFINITION 


Fig.  (3):  Defining  condition  for  the  supersymmetric  Higgs  particle. 


a  fiinHatnftntal  expression  that  connects  the  electron  neutrino  mass,  the  muon  neutrino  mass,  the 
masses  of  the  supersymmetric  Higgs  pair.  Importantly,  since  Eq.(39)  is  written  as  a  congruence,  it  holds 
simultaneously  for  the  masses  of  both  particles  comprising  the  supersymmetric  Higgs  system. 

An  explicit  evaluation  of  the  Higgs  masses  can  be  made,  since  ga  and  Po  are  known  [9,10,13]  from 
Eqs.(20)  and  (21),  and  the  critical  power  indices  X±  are  given  by  Eq.(32).  fri  order  to  simplify  the 
calculation,  we  can  use  the  smallest  primitive  root  of  Pa,  which  has  the  value  gmin  =  14,  since  the 
residues  corresponding  to  the  supersymmetric  Higgs  pair  are  independent  of  the  primitive  root  and 
depend  only  on  the  modulus  Pa.  Hence,  we  evaluate  the  congruence 


Pg-l 

8mfn  =  ®Higgs(®0^^a) 


(40) 


and  obtain  the  mass  parameter  of  Ae  jffiggs  boson  as 


=  5295067887317048385691240239639344309098972470130691907921115, 


(41) 


an  integer  composed  of  the  factors 


5,  13,  31,  109,  677, 25717,  165235181,  3158418359, 

5645180551, 470014220572174240309.  (42) 

The  mass  parameter  presented  in  Eq.(41)  gives  a  physical  mass  of  ~9.56  x  10‘*  GeV.  The  corresponding 
supersymmetric  Higgs  fermion  is  then  found  directly  by  the  condition  for  supersymmetry  through 


[f  =P  - 
Higgs  a  Higgs 

=146451816243788693417546684723041573 1421582430740004414836186,  (43) 


a  mass  parameter  having  the  fiictor  structure 

2, 3,  29, 37,  344353,  204901219,  14323448690905379, 
225085635424845454780404199  (44) 


and  representing  a  physical  mass  of  ~  2.64  x  lO'*  C5cV.  Consequentiy,  both  members  of  the 
supersymmetric  Higgs  doublet  possess  masses  >10^“  GeV.  These  results  given  by  Eqs.  (41)  and  (43) 
were  confirmed  by  a  calculation  using  the  primitive  root  g«.  Therefore,  since  the  two  neutrino  specie| 
have  a  lower  bound  on  their  masses  (<  30  meV)  and  the  Higgs  particles  are  both  superheavy  (>^ 
GeV)  we  find  that  Eq.(39)  contains  no  representative  of  an  intermediate  scale  (-100  GeV)  mass.  This 
outcome  contrasts  shaiply  with  customary  seesaw  relationships  [31]  developed  originally  from 
Lagrangian  field  theory  [11,32,33],  The  modular  seesaw  given  by  Eq.(39)  expresses  a  fundamentally 
different  physical  conclusion.  Further,  since  the  mathematical  procedures  used  to  compute  the  Higgs 
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mass  with  Eq.(40)  are  entirely  divergence-free,  no  renormalization  is  necessary,  a  second  result  that 
highlights  the  contrast  with  findings  derived  by  customary  quantum  field  theoretical  methods 
[6-8,34,35]. 

The  seesaw  congruence  given  by  Eq.(39)  expresses  a  deep  connection  between  the  concepts  of  space 
and  mass.  In  particular,  Eq.(39)  states  that  the  Higgs  system  specifies  the  set  of  generators 

{g  ,  g  ,  g'S  S"‘}  of  ®o®^®  iWustrated  in  Fig.  (2).  The  mass  parameters  corresponding  to  the 

a  p  a  P 

supersymmetric  Higgs  doublet  cannot  be  generators  themselves,  since  they  satisfy  the  quadratic 
congruence  given  by  Eq.(35)  and  are  thereby  barred  from  status  as  primitive  roots  of  P„.  Further,  since 
the  magnitudes  of  the  Higgs  residues  are  determined  in  Eq.(40)  solely  by  the  value  of  the  modulus  P„, 
the  seesaw  congruence  founds  the  generators  of  the  particle  masses  directly  on  the  large  scale 
characteristics  of  space.  Hence,  the  cosmic  seesaw  creates  a  new  archetype  by  expressing  the  fusion  of 
the  concepts  of  space  and  mass. 


E.  Higgs  Mass  and  the  Man^Body  Problem 

Quantum  field  theoretical  procedures  have  been  applied  to  the  computation  of  the  Higgs  mass 
[6-8,34,35].  Electroweak  theory  requires  quantum  loop  corrections  that  are  intrinsically  divergent  and 
lead  to  an  estimate  of  the  Higgs  mass  that  is  in  the  superhea^  regime  (>  10  GeV).  Basically,  the 
underlying  theory  has  a  many-body  character  and  these  corrections  are  a  manifestation  of  that  feature. 
Hence,  the  calculation  of  the  Higgs  mass  in  quantum  field  theory  is  a  basic  many-body  problem,  an 
outcorne  that  reflects  the  fundamental  nature  of  the  Higgs  mechanism  and  our  concept  of  the  vacuum 

[19].  . 

The  result  expressed  by  Eq.(40)  for  the  mass  parameters  of  the  supersymmetric  Higgs  doublet  is  an 
exact  value  obtained  by  mathematical  procedures  that  do  not  admit  a  perturbative  approach  [13],  With 
the  assumption  that  the  problem  is  properly  formulated,  the  only  result  obtainable  is  the  unique  correct 
solution.  In  the  approach  using  finite  fields  that  leads  to  Eq.(40),  this  occurs  because  of  the  large  number 
of  hidden  constraints  that  are  present.  Basically,  the  solution  is  gridlocked  by  the  combined  influence  of 
the  magnitude  and  the  factor  structure  of  the  integer.  This  exceptional  intolerance  to  computational 
approximation  follows  as  a  consequence  of  the  profound  influence  of  unity,  at  any  number  scale,  on  the 
factor  structure  of  integers.  Therefore,  since  there  exists  no  iieighborhood  ^out  the  exact  result  m  which 
an  approximation  is  meaning&l,  the  procedures  of  perturbation  theory  are  invalid. 

We  conclude  with  a  statement  of  the  obvious  conjecture,  namely,  that  the  counting  procedure 
intrinsic  to  the  theoretical  framework  provided  by  a  finite  field  ^tomtically  and  perforce  yield?  thg 
<|afaf!t  solution  to  a  proneriv  posed  manv-bodv  problem.  Further,  if  this  conjecture  holds  for  the  Hi^s 
particle,  we  may  expect  that  is  valid  for  any  system,  potentially  even  proteins.  In  the  case  of  proteins, 
the  verification  of  this  conjecture  would  allow  a  fundamental  reformulation  of  the  protein  folding 
problem,  a  central  question  of  immense  importance  in  biology  and  medicine. 

HL  Conclusions 

The  construction  of  inverse  states  in  a  finite  field  l>.enables  the  organi^tion  of  the  mass  scale  with 
fundamental  octets  in  an  eight-dimensional  index  space  that  classifies  particle  states  with  residue  class 
designations.  Conformance  with  both  CPT  invariance  and  the  concept  of  supersymmetry  follows  as  a 
direct  consequence  of  this  formulation.  Based  on  two  parameters  (P^  and  gj  that  are  anchored  on  a 
concordance  of  physical  data  embracing  both  large  scale  cosmic  paramrters  and  small  scale  couplings 
[9,10,13],  this  analysis  leads  to  (1)  a  prospective  mass  for  the  muon  neutrino  of  ~27.68  meV,  (2)  a  value 
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of  the  unified  strong-electroweak  cxiupling  constant  a*  =  (34.26)“^  that  is  physically  defined  by  the 
electron  and  muon  neutrino  mass  ratio,  and  (3)  a  seesaw  con^ence  connecting  the  Higgs,  the  electron 
neutrino,  and  the  muon  neutrino  masses.  Specific  evaluation  of  the  masses  of  the  corresponding 
supersymmetric  Higgs  pair  reveals  that  both  particles  are  superheavy  (>10  GeV).  No  renormdization 
of  the  Higgs  masses  is  necessary,  since  this  theoretical  approach  intrinsically  excludes  the  possibility  of 
divergences.  Further,  the  Higgs  fulfills  its  conjectured  role  through  the  seesaw  relation  as  the  particle 
defining  the  origin  of  all  particle  masses,  since  the  mass  parameters  of  the  electron  and  muon  neutrino 
systems,  together  with  their  supersymmetric  partners,  are  the  generators  of  the  mass  scale  and  the 
corresponding  index  space.  Finally,  since  the  computation  of  the  Higgs  masses  is  entirely  determined  by 
the  modulus  P^,  which  is  fully  defined  by  the  large-scale  parameters  of  the  universe  through  the  value 
[36]  of  the  universal  constant  G  and  the  requirement  [9,10,13,37]  for  perfect  flatness  (Q  =  1.0),  the 
cosmic  seesaw  congruence  provides  an  explicit  unifying  statement  for  the  concepts  of  mass  [1,2]  and 
space  [18,19].  As  inseparable  entities  in  this  picture,  these  two  physical  archetypes  coiistitute  an 
identity.  This  outcome  mirrors  the  fundamental  connection  between  space  and  geometry  in  general 
relativity.  _ _ _ 
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Abstract 

A  theoretical  analysis  of  laser-driven  collisional  ejection  of  inner-shell  electrons  is  presented  to 
explain  the  previously  observed  anomalous  kilovolt  L-sheil  x-ray  emission  spectra  from  atomic  Xe 
cluster  targets  excited  by  intense  sub-picosecond  248nm  ultraviolet  radiation  [A  McPherson  et 
al..  Nature  370,  pp.  631-634  (1994)].  For  incident  ponderomotively-driven  electrons 
photoionized  by  strong  above  threshold  ionization,  the  collisional  ejection  mechanism  is  shown  to 
be  highly  /-state  and  significantly  n-state  (i.e.  radially)  selective  for  time  periods  shorter  than  the 
collisional  dephasing  time  of  the  photoionized  electronic  wavefunction.  The  resulting  preference 
for  the  collisional  ejection  of  2p  electrons  by  an  ionized  Ap  state  produces  the  measured 
anomalous  Xe(L)  emission  which  contains  direct  evidence  for  (i)  the  generation  of 
Xe27+(^3f/i0)  and  Xe2*'*’(^3<jP)  ions  exhibiting  inner-shell  population  inversion  and  {it)  a 
coherent  correlated  electron  state  collision  responsible  for  the  production  of  double  2p  vacancies. 
For  longer  time  periods,  the  selectivity  of  this  coherent  impact  ionization  mechanism  is  rapidly 
reduced  by  the  combined  effects  of  intrinsic  quantum  mechanical  spreading  and  dephasii^  -  in 
agreement  with  the  experimentally  observed  and  extremely  strong  pump-laser  wavelength 
dependence  of  the  efficiency  of  inner-shell  (^)  vacancy  production  in  Xe  clusters  excited  in 
underdense  plasmas  [K.  Kondo  etal.,  J.  Phys.  B  30,  pp.  2707-2716  (1997)]. 
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Introduction 
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Several  methods  have  been  previously  evaluated  in  attempts  to  generate  coherent  radiation  in  the 
extreme  ultraviolet  (XUV)  and  x-ray  spectral  regions.  Most  notable  amongst  these  are;  high- 
order  harmonic  generation  [1-3];  discharge-pumped  soft  x-ray  lasers  [4];  laser-pumped 
recombination  [5,6],  electron-impact  ionization  [7,8],  and  hybrid  collision-recombination  [9] 
sources;  x-ray  [10]  and  ponderomotively-generated  Larmor  radiation  [11]  excited  emitters;  and 
Thomson  scattering  sources  [12],  However,  none  of  these  methods  have  been  successful  at 
generating  coherent  radiation  with  a  brightness  in  excess  of  10^  ys*‘(mrad)*^(mm)'^(0. 1%  BW)‘‘  at 
a  wavelength  X,  <  lOA;  that  is,  a  kilovolt  x-ray  brightness  level  sufficient  for  interferometric 
biological  imaging  capable  of  providing  an  atomic-level  resolution  visualization  of  the  molecular 
anatomy  of  cells,  tissues  and  organisms  in  the  natural  state  [13].  Indeed,  to  date,  kilovolt  x-ray 
imaging  applications  have  been  limited  to  static  [14]  and  dynamic  [15]  differential  absorption 
measurements  and  time-resolved  diffi-action  experiments  [16]  using  incoherent  radiation  fi-om 
laser-driven  plasmas.  Moreover,  practical  limitations  in  the  generation  of  high-power  sub-5-fs 
laser  pulses  and  the  spectral  merging  of  harmonic  orders  above  -200  may  preclude  the 
development  of  a  high  brightness  source  of  >lkeV  photons  by  high  harmonic  generation 
techniques,  while  a  recent  analysis  of  linear  and  nonlinear  (i.e.  multiple)  electron-photon 
interactions  [17]  suggests  that  Thomson  scattering  sources  may  never  reach  a  brightness  level 
sufficient  for  a  ground-breaking  interrogation  of  living  organisms. 

A  potentially  very  effective  method  of  producing  an  extremely  bright  source  of  kilovolt  x-ray 
radiation  would  be  to  generate  a  controlled  population  inversion  in  the  inner  shells  of  high-Z 
atoms,  i.e.  produce  specific  species  of  “hollow  atoms”  [18].  In  this  paper,  we  outline  a 
mechanism  for  achieving  precisely  this  type  of  population  inversion  in  laser-driven  plasmas.  The 
mechanism  is  based  on  the  selective  coUisional  ionization  of  inner-shell  electrons  by 
ponderomotively-driven  electrons  initially  photoionized  by  strong  above  threshold  ionization 
(ATI)  [19].  The  selectivity  rdies  on  the  fact  that  the  waveftinction  of  a  photoionized  electron 
distribution  maintains  its  phase  and  geometric  symmetry  for  a  time  shorter  than  the  charateristic 
coilisional  depharing  time  which,  by  comparison  with  electron-electron  scattering  rates  in 
semiconductors  [20,21],  is  of  the  order  of  one  femtosecond  (10->5s).  Consequently,  if  the 
photoionized  wavefimction  collides  with  an  ion  or  atom  within  this  time  period,  the  matrix 
element  associated  with  the  collision  will  be  dominated  by  the  intrinsic  dynamics  of  the  quantum 
mechanical  spreading  of  the  photoionized  atomic  orbital.  In  stark  contrast  to  a  “plane-wave” 
collision,  this  results  in  a  strong  angular  momentum  (/-  and  m-state)  and  a  significant  radial  (/t- 
state)  dependence  in  the  selectivity  of  the  interaction.  Las^-driven  classicaUy  fi’ee  electrons  in  the 
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plasma  [9]  cannot  exhibit  this  selectivity,  since  their  wavefunctions  contain  no  information  about 
the  atomic  orbitals  from  which  they  were  ionized. 

The  essential  features  of  the  selective  collisional  ionization  mechanism  will  be  discussed  in  Section 
III  using  the  interaction  of  ponderomotively-driven  4p  electrons  photoionized  from  Xe  with  the 
bound  state  electrons  of  Xe  ions  as  an  example.  The  results  of  this  theoretical  analysis  are  shown 
to  be  in  good  agreement  with  the  observed  anomalous  Xe(L)  spectra  described  in  Section  II, 
which  were  obtained  under  intense  ultraviolet  (UV)  248nm  laser  irradiation  of  5-20  atom  Xe 
cluster  targets  [22,23].  Specifically,  it  explains  the  strongly  prefered  generation  of  ^  vacancies  in 
Xe  cluster  targets;  a  process  leading  naturally  to  the  observed  evidence  for  (i)  a  laser-driven 
generation  of  Xe  ions  exhibiting  population  inversion  between  the  2p  and  3d  shells  and  (/V)  the 
production  of  multiple  2p  vacancies  [18].  Furthermore,  for  gaseous  targets  such  as  Xe  clusters, 
the  mechanism  displays  a  substantial  pump-laser  wavelength  scaling  favoring  the  use  of  shorter 
wavelength  driving  sources  -  in  agreement  with  the  experimental  determination  of  at  least  a 
scaling  in  the  efficiency  of  Xe  L-shell  (2.5-3. 1  A)  emission  between  248nm  and  800nm  laser 
excitation  [24,25].  The  possible  extension  of  this  technique  to  the  selective  generation  of  l5 
vacancies  in  Kr  for  strong  ~1  A  K-shell  emission  is  explored  in  Section  IV. 


II.  Xe  L'Shell  Emission 

Arguably  the  strongest  experimental  evidence  supporting  the  existence  of  an  anomalous  selective 
laser-driven  collisional  ionization  mechanism  comes  from  the  L-shell  (}d->2p)  hollow-atom 
emission  observed  from  atomic  Xe  clusters  [18]  under  intense  irradiation  (IQl^-lO^^/cm^)  by 
sub-picosecond  ultraviolet  (UV)  terawatt  (ITW  ==  lO^^W)  laser  pulses  at  248nm  generated  by  a 
hybrid  Ti:Sapphire/KrF*  laser  system  [26],  As  described  elsewhere  [18,23,25,27-29],  these  TW 
UV  laser  pulses  are  focused  by  an >73  optic  into  a  target  of  Xe  clusters  (typically  of  5-20  atoms) 
generated  by  a  cooled  pulsed  valve.  A  mica-crystal  von  H4mos  spectograph  equipped  with 
Kodak  RAR  2492  x-ray  film  and  appropriate  filters  was  used  to  collect  the  spectrum  of  the  x-ray 
radiation  emitted  by  the  irradiated  target.  As  shown  in  Figure  1,  the  kilovolt  Xe  L-shell  spectra 
detected  under  these  experimental  conditions  displays  the  characteristic  double  peaked  hollow- 
atom  structure  of  the  3d->7p  transition  [30].  It  is  important  to  note  that  this  Xe(L)  emission 
must  be  generated  by  a  collisional  mechanism  in  which  photoionized  electrons  are  accelerated  by 
the  laser  field  to  a  sufficient  energy  to  eject  an  inner-shell  2p  electron  via  Coulombic  interaction 
This  is  because  the  Xe  2p  state,  which  has  an  ionization  potential  of  ~8keV,  cannot  be  directly 
photoionized  by  above  threshold  ionization  (ATI)  [19]  at  irradiances  below  lO^^w/cm^ 
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However,  a  fully  relativistic  classical  model  of  the  photoionized  electron  trajectories  for  248nm 
irradiation  at  lO'^-lQi^w/cm^  [25]  has  shown  that  the  ponderomotively-driven  electrons  ionized 
front  more  weakly  bound  orbitals  can  return  to  their  cluster  of  origin  with  an  energy  of  10-20keV; 
a  value  sufficient  to  eject  2p  electrons  collisionally  and  thus  produce  the  observed  L-shell 
emission. 

A  striking  feature  of  the  spectra  shown  in  Fig.  1  is  the  strong  L-shell  emission  due  to  single  2p- 
shell  vacancies  from  Xe+9  ions  for  charge  states  of  ^  =  27  to  37  corresponding  to  2p^3d 
transitions  ranging  from  z.  filled  3</-shell  {2{f3(fiir-2p^3d\^)  to  one  containing  only  one  electron 
{^<r-2p^3d)  respectively.  This  observation  is  very  unusual  in  light  of  the  magnitudes  of  the 
cross-sections  evaluated  from  a  single-electron  plane-wave  analysis  for  the  colUsional  ejection  of 
inner-shell  electrons  with  an  ionization  potential  /p  >  250eV  [31,32].  As  shown  in  Figure  2,  a 
plane-wave  collisional  analysis  indicates  that  an  electron  incident  on  a  Xe  ion  with  a  collisional 
energy  E^qw  of  10-20keV  is  ~10  times  more  likely  to  eject  a  3d,  3p,  or  3s  electron  with  ionization 
potentials  of  ~2-3keV  than  a  2p  electron  with  /p  «  8keV.  This  analysis  is  clearly  in  stark 
contradiction  to  the  experimental  observations  of  Fig.  1  on  two  grounds.  First,  it  would  predict 
that  Xe(L)  emission  from  ions  should  be  approximately  an  order  of  magnitude  less  than  that 
from  Xe*^^')  ions,  ultimately  suggesting  that  on  average  all  3d  electrons  should  suffer  collisional 
ejection  before  a  single  2p  electron  is  ejected.  Consequently,  the  experimental  observation  of 
significant  L-shell  emission  from  Xe^^*  with  a  filled  3d  shell  (i.e.  a  2ffi3(fi<r-2p^3(fi^  transition) 
indicates  that  the  collisional  cross-section  for  ejection  of  a  2p  electron  could  be  at  least  an  order 
of  magnitude  larger  than  that  for  the  ejection  of  a  3d  electron  under  248nm  UV  irradiation 
atl0**-10*5W/cm2;  that  is,  precisely  the  inverse  of  the  scaling  associated  with  the  plane-wave 
collisional  analysis.  And  second,  it  suggests  that  a  significant  number  of  3p  and  35  vacancies 
should  also  be  generated  by  the  collisional  mechanism  ejecting  the  2p  electrons.  Yet,  strong 
emission  from  ionic  species  possessing  such  vacancies  has  not  been  observed. 

More  anomalous,  given  the  plane-wave  collisional  analysis,  is  the  observation  of  spectral  features 
associated  with  double  Tp  vacancies  in  the  range  2.60-2.6SA  (Fig.  1)  [18,22,23].  This  result  is 
even  more  astonishing  when  one  considers  the  fact  that  the  lifi^ime  of  a  single  2p  vacancy  is  only 
-l.Sfs  [32,30].  This  means  that  the  generation  of  2p  vacancies  in  the  collisional  ejection 
mechanism  must  be  occuring  at  a  rate  in  excess  of -lOi^s*';  that  is,  effectivdy  producing  multiple 
2p  vacancies  within  one  optical  cycle  of  the  248nm  pump-laso*  radiation  used  m  the  experiment. 
Previous  analyses  [28,33,34]  have  suggested  that  this  could  be  accomplished  by  ordered  or 
cooperative  ponderomotive  electronic  motions  in  which  Z  photoionized  electrons  behave  as  a 
single  quasi-particle  of  charge  Ze  with  an  effective  collective  energy  Z  times  their  individual  mean 
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energy.  For  example,  the  cooperative  action  of  3-4  electrons  under  248nm  UV  irradiation  at 
lO'S-lO’^W/cm^  would  certainly  enhance  the  probability  of  collisional  2p  vacancy  generation  by 
ensuring  that  the  quasi-particle  has  an  effective  collisional  energy  of  ~40keV,  i.e.  near  the  peak  of 
the  2p  cross-section  (Fig.  2).  However,  such  a  mechanism  fails  to  explain  the  apparent  selective 
generation  of  2p  vacancies  with  a  rate  exceeding  the  production  of  other  multiple  inner-shell 
vacancies  (e  g.  2s,  3s,  and  3p  vacancies)  in  the  Xe  ions. 

It  is  also  interesting  to  note  that  a  comparison  experiment  conducted  with  a  TW  femtosecond 
Ti:  sapphire  laser  system  operating  in  the  inJfrared  (IR)  at  SOOnm  did  not  reveal  the  generation  of 
multiple  2p  vacancies  in  the  observed  Xe  L-shell  emission  [24,25].  Moreover,  the  overall  strength 
of  the  L-shell  emission  was  observed  to  be  three  orders  of  magnitude  (~3000x)  weaker  under  IR 
irradiation.  Consequently,  any  plausible  collisional  ionizatioo  mechanism  must  (/)  be  selective  in 
nature  in  order  to  explain  the  anomalous  features  of  the  Xe(L)  emission  outlined  above  under  UV 
irradiation  and  (;/)  be  consistent  with  the  experimental  data  obtained  at  both  pump-laser 
wavelengths.  Evidently,  the  latter  effect  could  arise  directly  from  a  strong  scaling  with  the  laser 
excitation  wavelength,  specifically  the  pump-laser  radiation  period. 


III.  Selective  Collisional  Ionization 
in.  A.  Interaction  Geometry  and  Matrix  Element 

The  cross-section  for  any  collisional  interaction  is  proportional  to  the  modulus  square  of  the 
matrix  element  associated  with  the  interaction  mechanism  [35]  which,  for  an  initial  state  1/) 
coupled  to  a  final  state  |/)  through  the  interaction  potential  V,  is  written  as 

(I) 

In  our  case,  we  are  interested  in  the  collisional  section  of  bound  inner-shell  electrons  by  a 
ponderomotively-driven  photoionized  electron  through  the  action  of  the  Coulomb  potential.  The 
general  interaction  geometry  for  this  collision  is  shown  in  Figure  3.  The  bound  states  are 

defined  by  the  radial  coordinate  tj  centered  on  the  nucleus  of  the  target  atom  or  ion.  The 
displacement  vector  R  describes  the  trajectory  of  the  cento'  (centriod)  of  the  photoionized 
electron  distribution  (or  state)  during  the  motion  determined  by  the  laser's  oscillating 
ponderomotive  potential  and  the  Coulomb  potential  of  the  ionized  atomic  cluster  [25,36].  The 
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internal  radial  coordinate  of  the  photoionized  electron  waveflinction  |  \\i)  is  given  by  r2  so  that  R  + 

r2  defines  the  position  of  the  ionized  electron  with  respect  to  the  target  ion.  Within  this 
framework,  the  initial  state  |/)  can  be  written  as 

(r I /)  =  (|)^  (r,  )vi/(rj , /)exp[/ko  .(R  +  )]  ;  (2) 


that  is,  as  the  product  of  the  incident  ponderomotively-driven  time-dependent  photoionized 
wavefimction  (R+r2;/|\|/),  for  which  Kq  is  the  wavevector  associated  with  its  motion,  and  the 
stationary  inner-shell  electronic  states  of  target  atom  or  ion  |(|>y^.  In  the  limit  of  a  sufficiently  long 
time  after  the  interaction,  the  final  state  |/)  for  the  scattered  and  ejected  inner-shell  electrons  may 

be  well  represented  as  the  product  of  two  spherical  waves  emanating  from  the  origin  of  the 
coordinate  system  centered  on  the  target  atom  (Fig.  3): 


(r|/)  =  exp[/1t,.rj 

'  '  r,  Tj 


(3) 


Here,  kj  and  k2  are  the  wavevectors  of- the  ejected  and  scattered  incident  electrons  after  the 
collision  respectively.  Thus,  since  the  coUisional  ejection  mechanism  is  subject  the  momentum 
conservation  condition  ko  =  k,  +  k2  and  the  vector  fj  -  (R  +  r2)  describes  the  displacement 
between  the  two  electrons  involved  in  the  Coulomb  mteraction,  the  matrix  element  associated 
with  the  interaction  is  of  the  form 


(4) 


In  general,  the  full  exact  calculation  of  this  matrix  element  is  highly  complex  [37,38]. 
Accordingly,  for  the  scope  of  this  paper,  we  will  develop  a  suitable  approximate  analysts  bssfd 
upon  the  fact  that  the  I/r  dependence  of  the  Coulomb  potential  will  strongly  select  coUisioaal 
interactions  between  |  v)  and  for  which  R  *  0,  i.e.  interactions  for  which  the  center  of  the 
photoionized  state  |  y)  passes  close  to  the  nucleus  of  the  target  ion.  Consequently,  the  essential 
physics  behind  the  cross-section  for  inner-shell  coUisional  ejection  is  expected  to  be  obtained  from 
an  evaluation  of  |A^^f;(R  =  0)p.  In  particular  and  notwithstanding  this  approximatioii,  it  is 
immediately  evident  from  equation  (4)  that  the  strength  of  the  matrix  element  wUl  be  strongly 
dependent  upon  the  overlap  between  |  y)  and  that  is,  the  form  of  |  v)  since  the  bound  states 

l^y)  are  known.  SpedficaUy,  the  spherical  symmetry  of  the  incident  photoionized  wavefUnction 


\ 
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|v(/)  with  respect  to  the  target  wavefunctions  should  play  a  major  role  in  determining  the 

collisional  cross-sections  for  particular  inner-shell  ejection  mechanisms.  Hence,  as  in  atom-atom 
and  ion-atom  collisions  [39],  the  interaction  is  expected  to  have  strong  geometrical  features 
associated  with  the  electronic  charge  distributions  of  |  v)/)  and  | <l>y ^ . 

IILB.  Temporo-Spatial  Dynamics  of  Photoionized  Electron  Wavefunctions 

Before  the  collision  of  the  electron  with  the  target  atom  or  ion,  it  must  first  be  photoionized  and 
then  accelerated  by  the  pump  laser  field.  For  laser  irradiances  above  -lO^'^W/cm^,  such  as  those 
of  10i*-10*^W/cm2  required  to  generate  the  observed  Xe(L)  emission  (Fig.  1),  above  threshold 
ionization  (ATI)  is  the  dominant  ionization  process  [19].  Moreover,  under  our  experimental 
conditions,  the  incident  UV  pump-laser  radiation  is  observed  to  undergo  relativistic  self-focusing, 
resulting  in  channelled  beam  propagation  over  more  than  100  Rayleigh  ranges  [40].  Numerical 
simulations  of  the  self-focusing  mechanism  have  indicated  that  significant  self-steepening  of  the 
leading  edge  of  the  laser  pulse  occurs  over  irradiances  of  -'lO^'^-lO^W/cm^  as  it  propagates  down 
the  channel  [41].  As  a  result,  the  conditions  of  photoionization  experienced  by  the  irradiated 
atoms  or  ions  are  expected  to  be  highly  dynamic.  In  particular,  it  is  very  probable  that  the 
photoionization  mechanism  closely  follows  the  "ionization  ignition”  model  of  Rose-Petmck  et  al. 
[42]  where  entire  orbital  shells  may  be  ionized  by  ATI  at  the  peak  laser  field  in  each  half  cycle  of 
the  radiation.  Thus,  for  our  experimental  irradiances  of  10i*-10*^/cm2,  we  may  make  the 
ansatz  that  the  fi’ee  electronic  state  is  initially  a  replica  of  the  photoionized  bound  state  i.e. 

the  bound  state  is  simply  “tipped”  out  of  the  atomic  or  ionic  potential  by  the  perturbing 
ponderomotive  potential  of  the  strong  laser  field  in  ~0.1fs.  Consequently,  immediately  after 
photoionization  at  /  =  0,  we  may  write  within  this  "sudden  approximation”  that 

(rj;/  =  0|\|/>«(r,|<|>,).  (5) 

It  follows  that  the  initial  photoionized  electronic  state  (rj;/  =  0|\p)  possesses  the  angular  and 

radial  characteristics  of  namely  ^(r)T^(6,i|>)  where  R{r)  and  are  the  usual  radial 

eigenfunctions  and  spherical  harmonics  for  a  central  Coulomb  potential  respectively.  As  a  direct 
consequence,  one  would  expect  {r^d  =  0|\|/)  to  exhibit  substantial  p-v/vtt  symmetry  within  the 

irradiance  range  of  lO^^-lO^^/cm^  under  which  the  Xe(L)  emission  of  Fig.  1  was  observed, 
since  the  six-electron  Xe  4/7  state  with  an  ATI  threshold  of  ~2xi0i^/ciii2  is  one  of  the  last  states 
to  be  photoionized  at  these  levels  of  irradiance. 
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The  initially  ionized  state  =  0|\j/)  is  subsequently  accelerated  in  the  ponderomotive  potential 

of  the  oscillating  laser  field  Ecosto/  [25,43],  This  process  is  usually  described  by  the  expansion  of 
I  \|/)  in  terms  of  Volkov  states  |4»^^  [44],  the  eigenstates  of  the  ponderomotive  potential; 

•  (6) 

However,  the  expansion  coefficients  (^n|v)  are  generally  not  explicitly  evaluated  since,  as  we 

will  demonstrate,  the  ionized  electronic  wavefimction  evolves  into  phase-unrelated  Volkov  states 
I*.)  in  a  timescale  significantly  less  than  the  3-4fs  radiation  period  of  the  most  commonly  used 

high-power  IR  laser  systems  [45],  The  resulting  loss  of  information  on  both  the  initial  phase  and 
symmetry  of  |  y)  does  not  allow  any  possible  selectivity  in  collision  process  due  to  the  details  of 

the  spatial  overlap  of  the  returning  electron  wavefimction  (r2;/|v|/)  with  the  bound  states  |<|)y^  of 

the  target  atom  or  ion  to  be  treated.  Furthermore,  to  the  authors’  knowledge,  the  explicit  time 
evolution  of  the  photoionized  state  |v|/)  due  to  intrinsic  quantum  mechanical  spreading  and 
collisional  dephasing  (due  to  internal  scattenng  amongst  the  electrons  in  the  ionized  state  |  vp))  has 
not  been  previously  considered.  The  evaluation  of  these  effects  and  their  affect  on  collisional 
ionization  interactions  is  a  crucial  feature  of  the  analysis  presented  in  this  paper. 

The  intrinsic  quantum  mechanical  spreading  of  the  initial  wavefimction  (i‘2i^  =  ^o|v)  during  the 

acceleration  by  the  laser  field  can  be  readily  treated  within  its  internal  coordinate  system  r2  in  the 
single  non-interacting  electron  limit.  This  is  because  we  expect  the  spreading  of  |  \p)  to  be 

essentially  radial  since  =  0|\|;)  ori^nated  fi-om  a  bound  state  in  a  central  potential.  In 
other  words,  the  symmetry  dictated  by  the  spherical  harmonic  of  the  ionized  state  |4>y^  is 

predominantly  preserved  in  (r2;/|v)  since  no  torque  is  applied.  Thus,  to  a  first-order 
approximation,  the  time  evolution  of  |  q/)  may  be  written  as  [35] 

(  t^ 

(*’2;/|v)  =  (ri(/),0.<l>|v)  ;  i+-y  ,  (?) 

\  ) 

where  r2o  is  the  original  radial  coordinate  and  x  =  ma^jh  is  the  characteristic  time  for  the  radial- 
width  a  of  the  wavefiinction's  initial  probability  distribution  to  increase  to  V2  of  its  initial  value 
For  example,  the  wavefimction  for  an  electron  photoionized  fi’om  the  Xe  4p  state,  which  has  a 
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bound  radius  of  ~0.7A,  will  double  its  width  in  ~0.07fs,  i.e.  in  a  fraction  of  a  femtosecond!  More 
significantly,  when  the  photoionized  Xe  4p  electron  returns  to  collide  with  its  cluster  of  origin 
after  one  0.8fs  period  of  the  248nm  radiation,  |v^)  has  expanded  to  a  radius  of  -HA.  This  is  to 

be  compared  to  the  typical  ~10A  diameter  10-20  atom  Xe  cluster  size  used  as  the  target  material 
in  the  experiments  generating  the  observed  Xe(L)  emission  shown  in  Fig.  1.  Thus,  since  the 
photoionized  electron  tr^ectory  at  irradiances  of  10**-10i^W/cm2  is  known  to  come  back 
through  the  cluster  near  its  center  [25],  the  expanded  wavefiinction  is  optimally  sized  to  interact 
with  all  the  atoms  in  the  Xe  cluster.  This  optimal  geometrical  coupling  condition  for  collisional 
ionization  of  inner-shell  electrons  from  the  Xe  atoms  in  the  cluster  will  be  degraded  at  longer  laser 
excitation  wavelengths.  For  example,  upon  return  to  the  target  atoms  in  the  cluster  under  800nm 
IR  laser  excitation  (optical  period  of  2.7fs),  |  vj/)  would  be  ~3  times  broader  so  that  its  central 
probability  amplitude  (r2  «  0|  \|/)  is  reduced  by  a  factor  of  5-6.  In  turn,  this  should  result  in  more 
than  an  order  of  magnitude  reduction  in  the  cross-section  for  the  collisional  ejection  of  any  bound 
state  from  the  a  Xe  target  atomrion. 

In  addition  to  the  intrinsic  quantum  mechanical  spreading,  the  photoionized  wavefiinction  |  v)/)  will 
also  be  dephased  (i.e.  loose  its  initial  coherence)  by  collisional  scattering.  For  sufBciently  strong 
ATI,  where  all  electrons  in  an  entire  orbital  are  simultaneously  ionized  [42],  the  dominant 
scattering  mechanism  will  be  internal  electron-electron  collisions  within  the  photoionized 
state(r2;/|\|/).  Classically,  this  process  results  in  a  “Coulomb  explosion”  of  the  photoionized 

electron  distribution  [25]  which  converts  of  the  shape  of  the  initial  coherent  wavefiinction 
(r2;/  =  0|\|/)  into  incoherent  spherical  waves  ori^ating,  on  average,  from  the  center  of  the 

ionized  electron  distribution.  In  other  words,  this  process  is  primarily  responsible  for  the 
transition  into  the  regime  where  a  plane-wave  analysis  of  the  collisional  ejection  of  inner-shell 
electrons  becomes  valid. 

We  may  estimate  the  electron-electron  scattering  rate  Yee  for  this  dephasing  medianism  by 
comparison  with  the  large  body  of  experimental  work  on  carrier-carrier  collisions  in 
semiconductors.  In  particular,  the  observations  of  Portella  et  al.  [20]  and  Becker  et  al.  [21] 
indicate  that  the  electron-electron  scattering  rate  in  bulk  GaAs  increases  with  the  cube  root  of  the 
carrier  density,  i.e.  Yee  *  where  v4g,as  *  3.7xl05m/s.  Adjusting  for  the  effective 

electron  mass  of  GaAs  (m*  =  0.067Wo  [46],  where  is  the  free  electron  rest  mass),  we  obtain  a 
value  of  ~1  A/fr  for  Aq  in  a  vacuum.  Thus,  for  the  maximal  simultaneous  photoionization  by  ATI 
of  all  six  electrons  of  the  Xe  4p  state  in  an  initial  volume  with  a  radius  of  ~0.7A,  we  would  expect 
an  initial  scattering  rate  of  around  l-2xl0i  corresponding  to  a  characteristic  electron-electron 
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scattering  time  of  0.5-1  fs.  This  initial  scattering  time  is  roughly  a  factor  of  5  longer  than  the 
characteristic  time  associated  with  the  intrinsic  quantum  mechanical  spreading  of  the  photoionized 
wavefiinction.  Consequently,  the  Ap  wavefunction  photoionized  from  a  solitary  Xe  ion  spreads 
faster  than  it  dephases.  This  means  that  it  is  possible  for  a  significant  fi^ction  of  the  photoionized 
Ap  state  to  remain  coherent  during  the  ponderomotive  acceleration  over  one  radiation  period 
before  returning  to  interact  with  its  parent  Xe  cluster. 

The  fraction/of  the  photoionized  state  (r2;/|M/)  that  remains  coherent  as  a  function  of  time  may 

be  estimated  by  incorporating  the  reduction  in  the  rate  of  internal  electron-electron  collisional 
dephasing  caused  by  the  spreading  of  the  |  \|/)  into  a  rate  equation  of  the  form 


1 

2  Vi 


V. 


J 


(8) 


where  /lo  is  the  initial  electron  density  of  the  photoionized  wavefiinction.  The  solution  to  this 
equation,  subject  to  the  initial  generation  by  ATI  of  a  totally  coherent  state  (i.e./=  1  at  /  =  0),  is 


In  /  = 


t 

-  + 
T 


(9) 


For  the  photoionized  Xe  Ap  orbital  with  six  electrons  initially  in  a  radius  of  --O.TA,  this  analysis 
indicates  that  the  ionized  state  {r2;/|\t/)  will  be  30-40%  coherent  when  it  returns  to  interact  with 

its  cluster  of  origin  after  the  0.8fs  radiation  period  of 248nm  laser  excitation  [25].  As  a  result,  an 
accurate  description  of  the  collisional  ejection  mechanism  generating  the  spectrally  observed  7p 
vacancies  (Fig.  1)  must  incorporate  a  non-plane- wave  component.  Such  considerations  may  even 
be  necessary  for  800nm  IR  irradiation  since  equation  (9)  suggests  that  the  Ap  orbital  is  still  20- 
30%  coherent  after  the  longer  2.7fs  radiation  period,  even  though  (r2;/|\y)  is  now  substantially 

more  spread  out.  However,  this  intra-orbital  analysis  does  not  consider  dephasing  due  to  inter¬ 
orbital  electron-electron  scattering  which  vnll  begin  to  contribute  to  the  decoherenctng  of 
(*'2i^|v)  once  the  expanding  photoionized  states  fi’om  different  atoms  in  the  cluster  begin  to 

overlap.  The  scattering  rate  due  to  this  inter-orbital  interaction  is  estimated  to  be  0.5-lfii‘‘,  which 
suggests  that  the  electron  wavefiinction  may  be  almost  completdy  dephased  upon  returning  to 
interact  with  its  cluster  of  origin  when  the  driving  radiation  period  exceeds  ~2fs,  i.c  for  pump- 
laser  wavelengths  greater  than  -^OOnm. 
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IIJ.C.  Angular  Momentum  Selectivity  in  Collisional  Ejection  Mechanism 

The  above  analysis  of  the  temporo-spatial  dynamics  of  the  photoionized  wavefunction  clearly 
indicates  that  inner-shell  impact  ionization  due  to  incident  Xe  4p  electrons  under  248nm  laser 
excitation  could  be  primarily  affected  by  the  residual  coherence  of  the  spreading  (r2;t|n/)  rather 
than  plane-wave  collisions  resulting  from  a  dephased  (r2;/|H/).  Accordingly,  in  the  rest  of  this 

paper  we  will  examine  the  physics  behind  such  a  coherent  “orbital”-orbitaI  collisional  ionization 
process  [38,39], 


With  a  knowledge  of  the  form  of  the  expanding  and  ponderomotiveiy-driven  photoionized 
wavefunction  (r2;/|vi/)  given  in  equation  (7),  the  essential  physics  behind  the  matrix  element 

associated  with  the  collisional  ejection  of  an  electron  in  the  bound  state  (ri|<t>/)  by  the 

photoionized  electron  state  may  be  readily  determined  for  the  case  when  R  -  0.  In  this  limit, 
equation  (4)  reduces  to 


V2kl-|-2| 


exp[-/k,.(r, 


(10) 


which  may  be  solved  analytically  using  the  expansion  [47] 


L^O 


Ki^^L 


01) 


where  are  spherical  Bessel  functions  of  order  L,  are  Hankel  functions  of  the  first 

kind  of  order  I,  and  k  =  -  kicos^  with  ^  equal  to  the  angle  between  ki  and  the  vector  n  •  r^.  As 
with  other  collisional  interactions  [35,37],  the  matrix  element  is  a  function  of  the  details  of  the 
collision  dynamics,  that  is  ki  through  ic. 

< 

We  will  begin  by  considering  the  collision  physics  associated  with  the  angular  contributioa  to 

A^R  -  0)  in  the  integrand  of  equation  (10).  Here,  it  is  important  to  note  that,  unlike  the 
direction  of  ri  -  r2,  the  direction  of  the  ejected  electron’s  momentum  ftk,  cannot  be  explicitly 

related  to  the  integration  variables  variables  9i,  63,  and  Instead,  the  allowed  directions  of 

k|  are  determined  solely  by  the  momentum  and  energy  conservation  in  the  collision.  As  a  result. 
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the  angle  (3  between  kithe  vector  ri  -  rz  must  be  considered  as  a  random  variable  in  the  evaluation 
of  Ai^R  =  0)  so  that,  for  example,  the  energy  conservation  condition  for  the  ejection  of  an  inner- 
shell  electron  with  an  ionization  potential  Ip  may  be  written  as 


(12) 


since,  on  average,  cos*  P  =  | .  It  is  then  manifest  that  the  form  of  the  angular  integrations 

involved  in  Mji  imply  a  strong  angular  momentum  selectivity  to  the  collisional  ejection  of  inner- 
shell  electrons.  Specifically,  if  one  is  interested  in  the  collisional  ejection  of  an  electron  from  a 
bound  state  |(t>y  ^  with  an  orbital  angular  momentum  quantum  number  /,  then  the  insertion  of 

equation  (11)  into  equation  (10)  ensures  that  only  the  L  =  /  term  of  the  Coulomb  expansion 
contributes  to  (the  6|  integration),  and  thence,  that  the  collisional  cross-section  is  non-zero 
only  if  the  incident  photoionized  electronic  state  (r2;r|vi/)  has  an  angular  momentum  component 

of  Iti  in  the  collision  (the  02  integration).  Alternatively,  if  the  ponderomotively-driven  incident 
wavefunction  (r2;/|\{/)  is  predominantly  a  state  with  a  particular  orbital  angular  momentum 

quantum  number,  then  strong  /-state  selectivity  will  result  in  the  collisional  ejection  mechanism. 


This  type  of  /-state  selectivity  in  the  collisional  ionization  process  is  clearly  in  agreement  with  the 
observation  of  strong  Xe  L-shell  emission  from  Xe“^  ions  under  intense  10**-10*^/cm2  UV 
laser  excitation  (Fig.  1).  In  particular,  unlike  the  plane-wave  collisional  analysis,  it  explains  the 
generation  of  2p  vacancies  in  Xe  ions  which  have  all  thdr  35  and  a  substantial  fraction  of  their  Zd 
electrons  remaining.  This  is  because,  for  laser  irradiances  ^.5xl0'*W/cm2,  sufficiently  strong 
ATI  can  ^ultaneously  photoionize  the  entire  Xe  Ap  sheU  to  generate  an  initial  6-electron  state  of 
the  form 


(13) 


After  approximately  one  period  of  the  driving  laser  radiation,  the  time-evolving  state  (r2>/|v) 

returns  to  interact  with  the  Xe  atoms  in  its  cluster  of  origin  [25].  Under  248nm  UV  laser 
excitation  with  a  0.8fs  radiation  period,  the  30-40%  of  the  returning  wavefunction  (r2;/|  v)  which 
is  still  coherent  with  its  initial  photoionized  Ap  state  (r,;/  =  0|\p)  will  predominantly  coUisionally 

eject  a  bound  electron  from  a  state  with  the  same  spherical  harmonic  i.e.  a  /estate 

electron  with  /  =  1.  The  matrix  elements  with  all  other  bound  states  with  /  *■  1  will  be 
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severely  supressed  since  the  spreading  of  (r2;/|v|/)  has  mostly  preserved  its p-wave  character.  For 
longer  radiation  periods,  such  as  the  w3fs  period  of  near  IR  laser  radiation,  this  /-state  selectivity 
in  the  collisional  ejection  process  is  expected  be  sharply  reduced  by  intra-  and  inter-orbital 
electron-electron  scattering  which  destroy  the  symmetry  of  the  initial  photoionized  state 
(fj;/ =  0|v]/),  thus  rapidly  and  progressively  returning  the  inner-shell  ejection  process  to  the 

customary  plane-wave  regime.  However,  the  coherent  nature  of  the  collisional  interaction  may 
not  be  completely  destroyed,  as  evidenced  by  the  observation  of  weak,  yet  identifiable,  L-shell 
emission  from  Xe  ions  with  charge  states  between  27+  and  29+  (Fig.  1);  i.e.  Xe  ions  with  filled  35 
and  mostly  filled  "id  orbitals. 

III.D.  Multiple  Inner-Shell  Vacancy  Production 

It  is  interesting  to  note  that,  in  addition  to  the  /-state  selectivity,  m-state  selectivity  is  also 
expected  in  the  coherent  inner-shell  collisional  ejection  mechanism.  This  arises  because  the  strong 
electric  field  fi^om  the  exciting  laser  radiation  imposes  a  preferred  axis  to  the  interaction  through 
the  Stark  effect.  As  a  result,  the  angular  integrations  over  <|>i  and  <{>2  ensure  that  a  non-zero  matrix 
element  A^R  =  0)  for  the  collisional  ejection  of  an  electron  fi-om  a  particular  bound  state  |(t>y^ 
exists  only  when  the  incident  photoionized  state  (r2;/|v{/)  has  a  component  of  the  same  azimuthal 
angular  momentum  as  Since  the  Coulomb  interaction  is  spin  invariant,  this  has  an  important 

consequence  for  the  shell-selective  collisional  ejection  process:  namely,  that  it  is  possible  for  the 
two  electrons  in  the  same  /-  and  w-state  of  (r2;/|v)  to  act  collectively  (or  coherently)  [28,33,34] 

to  collisionally  eject  two  inner-shell  electrons  with  the  identical  orbital  angular  momentum  and 
magnetic  quantum  numbers.  Such  a  collision  between  two  correlated  anti-symmetric  spin  states 
of  the  form  “  |-  i .  i))  the  target  and  photoionized  wavefunctions  has  the  same  value 

of  the  matrix  element  as  the  single  electron-electron  collisional  ejection  mechanism  in  the 
absence  of  dephasing  efifects.  Thus,  the  generation  of  double  m-state  vacancies  by  tlus  selective 
coherent  process  is  only  4  times  less  likely  than  the  production  of  single  vacancies  -  there  being  4 
possible  single  electron-electron  collisions  for  each  /-  and  m-state  selective  interaction.  This  ratio 
of  single  to  double  vacancy  production  is  to  be  contrasted  to  a  ratio  of  lO’-lO^  for  the  collision  of 
high  Z  atoms  (Z  >  20)  with  high  energy  quanta  [48].  This  difference  is  expected  since  the 
coherent  ejection  mechanism  is  more  reminiscent  of  the  orbital-oibital  interactions  of  ion-atom 
collisions  [38,39]. 

At  first  glance,  one  would  expect  the  collision  producing  a  double  vacancy  to  be  more  affected  by 
dephasing  effects  than  the  single-vacancy-generating  interaction  since,  for  the  former,  both 


13 


electrons  in  the  particular  m-state  of  (r2;/lv|/)  must  remain  coherent.  However,  to  a  good  first 
approximation,  this  is  not  the  case  within  a  single  photoionized  orbital,  because  the  most  likely 
electron-electron  collision  will  occur  between  electrons  in  the  closest  spatial  proximity;  that  is, 
between  the  spin  pair  in  the  same  m-state  wavefunction.  Clearly,  this  dominant  collisional 
dephasing  process  will  have  an  equal  effect  on  the  rates  of  both  single  and  double  vacancy 
production,  resulting  in  an  approximately  fixed  ratio  of  single  to  double  vacancy  production  of  4. 
For  atomic  cluster  or  solid  targets,  the  overlap  between  expanding  photoionized  wavefunctions 
from  neighboring  atoms  or  ions,  and  consequent  inter-orbital  dephasing  collisions,  will  cause  this 
ratio  to  increase.  However,  the  experimental  evidence  of  Fig.  1  indicates  that  the  ratio  of  the 
peak  emission  strengths  from  Xe  ions  with  single  and  double  2p  vacancies  is  only  ~5,  suggesting 
that  the  inter-orbital  scatteiing  rate  may  be  less  than  O.Sfs*'  and  thus  not  a  strong  effect  for  short- 
period  UV  laser  excitation.  The  observation  of  strong  double  ^  vacancy  production  in  Xe 
clusters  excited  by  248nm  radiation,  therefore,  provides  direct  evidence  for  a  coherent  (or 
correlated)  collisional  ejection  mechanism.  Moreover,  the  lack  of  the  same  spectral  signiture 
under  800nm  IR  laser  excitation,  where  the  longer  radiation  period  only  serves  to  further  “mix” 
and  dephase  the  expanding  ionized  states  {r2;t|v|/),  is  dearly  consistent  with  such  as  mechanism. 

HLK  Radial  Selectivity  in  the  Matrix  Element 

To  fully  account  for  the  experimentally  observed  anomalous  Xe  L-shell  emission  spectra  (Fig.  1), 
the  proposed  coherent  inner-shell  ionization  mechanism  must  also  display  significant  n-state  (or 
radial)  selectivity  in  addition  to  the  strong  angular  momentum  selectivity  discussed  above. 
Specifically,  it  must  explain  the  observed  lack  of  strong  3p  ejection  when  the  incident 
ponderomotively-driven  state  (r2;/|v|/)  originated  from  a  photoionized  4p  orbital  at  UV 

irradiances  of  IQi^-lO^^/cm^.  This  follows  because  the  experimental  designation  and 
calculational  verifications  of  the  charge  states  in  the  observed  Xe(L)  emission  from  2p  vacancies 
support  the  conclusion  that  the  3p  shell  is  intact.  In  other  words,  the  collisional  cross-section  for 
the  ejection  of  a  ^  electron  by  the  proposed  mechanism  should  be  greater  than  that  for  the 
ejection  of  a  3p  electron,  a  result  that  stands  in  contradiction  to  the  conventional  plane-wave 
collisional  anatysis  (Fig.  2).  Since  this  type  of  selectivity  cannot  be  /  or  m  dependent,  it  must  be 
due  to  radial  dependence  of  the  matrix  element  =  0). 

Evaluation  of  the  radial  contribution  to  the  matrix  element  for  R  -  0  requires  knowledge  of 
the  radial  forms  of  both  the  initially  ionized  wavefiinction  (r,;/  =  0|\|/)  and  the  bound  states 

of  the  target  ion.  In  order  to  simplify  this  calculation,  we  will  approximate  these  states  by  suitably 
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scaled  (to  account  for  screening)  hydrogenic  waveflinctions  [35],  Figure  4  compares  the 
computed  normalized  hydrogenic  wavefunctions  rR(r)  (dashed  curves)  with  those  obtained  from 
an  exact  calculation  (solid  curves)  [49]  for  the  3p,  and  4p  radial  orbitals  of  atomic  Xe.  The 
relatively  good  agreement  in  all  three  cases  allows  the  matrix  element  A^R  =  0)  to  be  evaluated 
analytically  to  a  high  precision.  The  slight  disagreement  on  the  form  of  the  outer  lobe  of  the  4p 
orbital  accounts  for  less  than  a  2%  error  in  the  calculation,  since  (/)  the  1/r  dependence  of  Vg 
supresses  its  contribution  relative  to  the  first  two  lobes  and  (//)  the  overlap  of  the  outer  lobe  of  the 
4p  state  with  the  3p  and  2p  orbitals  is  minimal.  Moreover,  the  rapid  radial  expansion  of  the  4p 
orbital  after  ATI  ensures  that  effectively  only  the  inner  two  of  its  well-represented  radial  lobes 
contribute  to  the  matrix  element  after  a  time  /  >  2x  «  0.1 4fs.  In  fact,  for  small  laser-excited  Xe 
clusters,  where  the  photoionized  electrons  return  to  interact  with  their  cluster  of  origin  after  one 
radiation  period  (i.e.  after  a  time  tlx  >  10)  [25],  essentially  only  the  first  radial  lobe  of  the 
expanding  4p  orbital  overlaps  with  both  the  bound  Ip  and  3p  states  for  R  =  0.  In  addition,  it 
should  be  noted  that  the  use  of  the  radial  forms  of  the  2p,  3p,  and  4p  states  for  a  Xe  atom  rather 
than  an  Xe+«  ion  of  the  appropriate  charge  should  not  affect  the  evaluation  of  the  matrix  element 
significantly.  This  is  because  the  basic  shdpe  of  both  the  tightly-bound  inner-shell  2p  and  3p  states 
and  the  innermost  first  two  lobes  of  the  4p  state,  which  dominate  the  collisional  ejection 
mechanism,  should  not  be  significantly  perturbed  for  typical  charge  states  of  ^  «  +30  generated  at 
the  experimental  excitation  irradiances  of  10^*-10*^/cm^. 

The  result  of  evaluating  the  radial  integration  of  the  matrix  element  for  the  collisional  interaction 
in  equation  (8)  using  the  hydrogenic  Tp,  3p,  and  incident  expanding  4p  waveftuictions  is  shown  in 
Figure  5.  The  value  of  |A^R  =  0)p,  which  is  directly  proportional  to  the  collisional  cross-section 
for  the  selective  ejection  mechanism,  is  plotted  as  a  function  of  the  normalized  time  t/x  given  by 
equation(7)  and  x  =  kicosp  for  the  collisional  ejection  of  the  bound  ^  and  3p  states  in  Figs.  S(a) 
and  (b)  respectively.  Three  features  are  readily  apparent.  First,  in  contrast  to  the  plane-wave 
collisional  analysis  (Fig.  2)  [32],  the  cross-section  for  3p  ejection  is  less  than  that  for  the  impact 
ionization  of  a  ^  electron  at  equivalent  values  of  tlx  and  k.  Second,  as  expected,  the  aoss- 
section  for  the  impact  ionization  of  inner-shell  electrons  decreases  in  time  due  to  the  reduction  in 
the  collisional  overlap  integral  caused  by  the  temporal  spreading  of  the  photoionized 
state(r2;/|M/).  Third,  the  cross-section  is  also  a  strong  function  of  k  with  slower,  less  energetic 

collisions  being  favored  for  tlx  ^  10.  Physically,  the  latter  trend  may  be  understood  in  terms  of 
the  collision  time  -  a  less  energetic  collision  allowing  for  a  longer  time  period  for  the 
wavefunctions  representing  the  electron  states  to  interact.  Mathematically,  it  is  a  direct 
consequence  of  the  form  of  the  Bessel  function  product  7,(Kr)A['^(Kr)  in  equation  (11)  which. 
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for  sufficiently  large  values  of  k,  has  a  more  rapid  oscillation  with  r  than  the  2p  and  3p  bound 
state  wavefunctions  (Fig.  4),  thus  reducing  the  value  of  =  0)|2. 

To  quantify  the  first  two  observations,  we  will  examine  in  more  detail  the  interaction  conditions 
leading  to  the  L-shell  spectra  of  Fig.  1  emitted  by  Xe  clusters  under  248nm  UV  and  800nm  IR 
laser  excitation.  Form  the  relativistic  analysis  of  Ref.  25,  a  photoionized  electron  state  {r2;/|vi/> 

will  return  to  its  cluster  of  origin  with  a  colUsion  energy  of-15keV  after  approximately  one 
radiation  period  for  both  pump-laser  wavelengths  at  the  experimental  irradiances  of  10»8. 
10«9w/cm2.  With  this  information,  the  most  probable  average  value  of  k  for  an  inner-sheU 
ejection  process  may  be  calculated  using  equation  (12)  and  assuming  the  equipartition  of  the 
energy  available  after  the  collision  between  the  scattered  and  ejected  electrons.  For  the  ejection 
of  a  Ip  electron  with  an  ionization  potential  /p  ~  8keV.  this  calculation  gives  k  «  21  A’*,  whereas 
for  3p  vacancy  production  we  have  k  *  28A'‘  since  I,  -  2.8keV  in  this  case.  For  these  values  of 
K,  Figure  6(a)  displays  the  temporal  dependence  of  l,V/yj(R  =  0)|2  for  the  ejection  of  both  2p  (solid 

line)  and  Zp  (dotted  line)  electrons  by  an  interaction  with  the  incident  expanded  and  coherent  4p 
state. 


Most  notably,  Fig.  6(a)  clearly  indicates  that  the  cross-section  for  the  collisional  ejection  of  a  ^ 
electron  by  the  ponderomotively-driven  and  expanding  4p  orbital  is  always  at  least  two  orders  of 
magnitude  greater  than  for  the  ejection  of  a  3p  electron.  This  additional  radial  or  /i-sute 
selectivity  is  a  direct  consequence  of  the  shape  of  the  3p  radial  wavefunction  relative  to  that  of  the 
2p  state;  namely,  the  bipolar  nature  of  the  3p  orbital  supresscs  the  matrix  element  associated  with 
the  collisional  interaction.  In  other  words,  the  qection  of  electrons  from  states  with  the 
fundamental  single-peaked  radial  wavefunction  is  strongly  preferred  in  this  selective  coherent 


interaction.  As  shown  by  the  plot  in  Fig.  6(b),  this  is  in  fact  true  for  ail  values  of  x  and,  hence, 
incident  collision  energies  for  the  4p  state.  Moreover,  the  ratio  of  7p  to  3p  collisional  cross- 
sections  for  this  coherent  interaction  is  also  relatively  insen^ve  to  the  exdtation  wavelength  for 
X  >  5A  ‘.  This  result  is  clearly  consistent  with  the  observed  Xe(L)  emission  under  intense  laser 


excitation,  because  virtually  no  3p  vacancy  generation  would  be  expected  even  for  multipie  2p 
vacancy  produption;  an  outcome  in  agreement  with  the  charge  labeling  in  Fig.  1  for  both  the  IR 
and  UV  excitation  wavelengths.  Furthermore,  the  calculated  ratio  between  the  ^  and  3p 
collision  cross-section  represents  at  least  a  three  orders  of  magnitude  reversal  over  the  cross- 
section  ratio  evaluated  from  a  plane- wave  analysis  (Fig.  2)  [32]. 


Figure  6(a)  also  reveals  a  strong  pump-laser  wavelength  dq)endence  to  the  selective  fnlliskmal 
ionization  process;  namely,  that  the  coherent  orbital-orbital  collision  is  much  more  probable  for 
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shorter  pump-laser  wavelengths  (i.e.  short  radiation  periods).  Specifically,  since  x  »  0.04fs  for  the 
Ap  Xe  state  with  an  initial  bound-state  radius  a  of  -0.7A,  the  value  of  IA^(R  =  0)l2  for  the 
ejection  of  a  Ip  electron  is  ~350  times  greater  under  248nm  UV  laser  excitation  (collision  at  Hi  = 
19)  than  for  the  800nm  pump-laser  wavelength  (collision  at  th  =  63).  The  intrinsic  quantum 
mechanical  spreading  of  the  photoionized  Ap  state  causes  a  similar  reduction  in  the  cross-section 
for  3/?  impact  ionization  as  the  excitation  wavelength  increases.  The  inclusion  of  collisional 
dephasing  within  the  incident  photoionized  state  (r2i^|v).  which  destroys  the  selectivity  of  the 

inner-shell  ejection  process,  further  magnifies  this  adverse  dependence  on  the  radiation  period. 
For  an  estimated  average  electron-electron  scattering  rate  of  the  order  of  Ifs'*  (Section  III.B),  we 
expect  a  ratio  in  the  selective  collisional  cross-sections  of  -lO^  favoring  the  use  of  248nm  UV 
over  800nm  ER  laser  excitation.  This  cross-section  ratio  is  in  good  agreement  with  the 
experimental  observation  of  a  strong  pump-laser  wavelength  dependence  in  the  Xe  L-shell 
emission  [24,25]  which  indicates  that  the  shorter-period  UV  radiation  is  -3000  times  more 
effective  at  generating  2p  vacancies  in  S-20  atom  Xe  cluster  targets  than  equivalent  terawatt  IR 
(SOOnm)  radiation.  An  exact  calculation  of  all  collisional  dephasing  rates  in  expanding 
photoionized  states  (r2;/|v|/)  and  a  full  evaluation  of  the  matrix  element  mtegrated  over  all 

values  of  R  in  a  more  rigorous  analysis,  both  of  which  are  beyond  the  scope  of  this  paper,  may 
yield  a  closer  agreement  with  the  experimental  results.  Nevertheless,  the  approximate  treatment 
unambiguously  shows  that  the  wavelength  dependence  is  strong  and  that  the  use  of  sufficiently 
long  pump-laser  wavelengths  will  eliminate  the  selectivity  of  the  interaction. 

It  is  intersting  to  note  that  since  the  “ip  cross-section  is  at  least  one  order  of  magnitude  less  than 
the  “Ip  cross-section  for  timescales  tlx  <  5,  it  may  be  possible  to  circumvent  the  strong  pump-laser 
wavelength  dependence  of  L-shell  emission  by  ensuring  that  the  time  between  the  photoionization 
by  ATI  and  the  collisional  ejection  is  much  less  than  the  period  of  the  driving  radiation.  This 
immediately  implies  that  the  ionized  state  (r2;/|y)  must  interact  within  the  filrst  quarter  of  the 

laser  radiation  period  with  an  atom/ion  that  is  remote  from  its  atom  of  origin  as  could  be  the  case 
in  a  solid  target  material  or  in  large  (>100)  atom  clusters.  However,  for  the  selective  mechanism 
to  be  effective,  the  laser-driven  state  {r2;/|'i/)  cannot  be  significantly  dq)hased  by  collisions  with 

other  atoms/iops  in  the  target  material  prior  to  reaching  the  incident  kinetic  energy  necessary  for 
the  collisional  ejection  of  a  particular  inner-shell  electron.  Such  a  condition  is  expected  to  be 
difficult  to  achieve  unless  the  ponderomotive  potential  is  sufficiently  strong  to  allow  acceleration 
to  the  required  kinetic  energy  over  inter-atomic  distances.  Nonetheless,  precisely  this  type  of 
interaction  has  already  been  observed  in  the  generation  of  20g  vacandes  in  N2  [SO]  and  in  the 
more  recent  detection  of  strong  Ba(L)  emission  from  solid  BaFj  targets  [51]  under  excitation  with 
intense  sub-picosecond  248nm  radiation. 


17 


IV.  Inner-Shell  Population  Inversion 


The  selective  coilisional  ionization  mechanism  described  above  is  clearly  well  suited  to  the 
efficient  generation  of  inner-shell  vacancies  and,  hence,  the  population  inversion  required  for  the 
development  of  coherent  kilovolt  x-ray  radiation  sources.  In  fact,  given  the  quantum  state 
degeneracies  of  the  2p  and  2d  levels,  the  experimental  observation  (Fig.  1)  of  anomalously  strong 
Xe(L)  emission  due  to  2/;  vacancies  from  both  Xe27+(:^3f/*0)  and  is  already 

direct  evidence  for  the  efficient  generation  of  inner-shell-inverted  ionic  species  in  laser-driven 
plasmas.  Whether  or  not  population  inversion  can  be  produced  is  dependent  upon  the  relative 
rates  of  generation  of  the  upper  and  lower  states  of  the  considered  transition;  for  example,  the 
electron  configurations  2p^2d^^  and  2f^2(fi  respectively  for  L-shell  emission  from  Xe^'^^. 
Nonetheless,  the  above  analysis  of  the  experimental  results  (Fig.  1)  indicates  that  a  UV  laser- 
pumped  Xe  L-shell  system  has  great  potential  as  a  coherent  SA  («4keV)  x-ray  radiation  source. 
Furthermore,  the  strong  relativistic  self-focusing  effects  observed  in  Xe  cluster  targets  under 
248nm  pump-laser  irradiation  [52]  may  produce  an  optimal  guided-channel  geometry  [29,40] 
extending  over  several  Rayleigh  ranges  for  both  the  optical  pumping  and  the  x-ray  emission. 

The  extension  of  this  coilisional  process  to  the  generation  of  more  energetic  x-ray  radiation  would 
require  the  use  of  either  higher  Z  elements  for  the  equivalent  L-shell  emission  or  K-shell  emission 
through  the  \^2p^<r-\^2f^  transition  in  elements  with  Z >  30.  While  the  former  is  complicated  by 
the  engineering  of  suitable  metal  cluster  or  other  composite  and  compound  (e  g.  UF^)  targets,  the 
latter  possibility  appears  to  offer  more  practical  flexibility.  A  similar  selective  coilisional 
ionization  process  for  K-shell  emission  would  require  the  photoionization  by  ATI  and  subsequent 
ponderomotive  acceleration  of  the  35  state  to  collisionally  eject  an  electron  from  the  l5  state 
within  one  radiative  period.  Unfortunately,  for  Xe(K)  emission,  a  very  high  laser  irradiance  of 
~3xl020W/cin2  to  ionize  the  35  state  would  be  required.  A  attractive  alternative  target  atomic 
species,  which  could  produce  very  useful  -1 A  K-shell  radiation  [13],  is  Kr  for  which  ATI  of  the 
35  state  requires  a  laser  irradiance  of  '-2xl0i^/cin2  -  a  value  associated  with  the  channeling  of 
248nm  radiation  in  underdense  plasmas  [29]. 

Figures  7(a)  and  (b)  show,  respectively,  the  results  of  evaluating  the  modulus-square  of  the 
explicit  /-state  selective  matrix  element,  i.e.  =  0)p,  for  the  coilisional  ejection  of  bound  1 5 
and  25  electrons  in  Kr  by  ponderomotively-driven  35  electrons  as  a  function  of  k  and  the 
normalized  time  i/x  for  the  expanding  photoionized  orbital.  Again,  the  collision  cross-seaion 
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decreases  in  time,  due  to  the  spreading  ionized  3s  wavefunction,  and  with  increasing  collision 
energy,  i.e.  k.  However,  in  this  case,  it  is  not  immediately  apparent  whether  it  is  more  probable  to 
collisionally  eject  a  Is  or  a  2s  electron  in  the  coherent  interaction.  The  answer  to  this  question  lies 
in  the  details  of  the  collision  dynamics. 

At  a  laser  irradaince  of  ~2xlO*9W/cm2,  the  relativistic  kinetic  analysis  of  Ref  25  indicates  that  an 
electron  photoionized  by  ATI  will  return  to  a  small  (~10  atom)  Kr  cluster  with  an  energy  of 
~20keV  after  approximately  one  radiation  period.  Since  the  ionization  potentials  of  the  Is  and  2s 
states  are  ~17keV  and  ~3.2keV,  the  most  probable  average  values  for  x  are  14  and  33  for  the 
ejected  Is  and  2s  electrons  respectively.  For  these  values  of  x,  the  temporal  dependence  of 
=  0)|^  for  the  ejection  of  both  Is  (solid  line)  and  2s  (dotted  line)  electrons  by  an 

interaction  with  the  ponderomotively-driven  coherent  3s  state  is  shown  in  Figure  8.  It  is 
immediately  apparent  that  again  the  impact  ionization  of  the  inner  single-Iobed  bound  Is  state  is 
favored  over  the  ejection  of  an  electron  from  the  2s  state,  but  only  by  roughly  one  order  of 
magnitude.  The  fact  that  the  magnitudes  of  the  two  matrix  elements  are  much  closer  for  this  s- 
state  interaction  than  equivalent  elements  in  the  pervious  p-state  collisions  is  due  to  the  form  of 
the  Laguerre  polynomials  for  1  =  0;  namely,  that  the  relative  width  of  the  first  radial  lobe  of  the  s- 
states  is  less  than  that  of  the  p-states.  Consequently,  the  overlap  of  the  expanding  photoionized 
3s  wavefunction  with  the  relatively  large  second  lobe  of  the  2s  state  dominates  the  matrix  element 
despite  the  1/r  variation  of  the  Coulomb  potential.  Nevertheless,  the  analysis  suggests  that  it  may 
be  possible  to  generate  eflSciently  through  laser  excitation  Kr  ions  exhibiting  inversion  between 
the  2p  and  Is  states  without  significant  influence  from  a  parasitic  2s<-2p  transition. 

The  efficiency  with  which  such  an  inner-shell-inverted  ionic  species  are  produced  is,  however, 
again  a  strong  function  of  the  pump-laser  wavelength  due  to  the  spreading  of  the 
ponderomotively-driven  3s  state  over  one  radiation  period.  In  this  case,  the  smaller  initial  *-0.5 A 
radius  of  the  Kr  3s  orbital  [49]  ensures  that  the  photoionized  wavefunction  spreads  at  almost 
twice  the  rate  of  a  Xe  4p  orbital  ionized  by  ATI.  As  a  result,  for  Kr  K-shell  emission,  the  ratio  of 
the  Is  coUisional  cross-sections  at  pump-laser  wavelengths  of  248nm  (//t »  40)  and  SOOnm  {t/x  « 
130)  only  favors  the  UV  laser  radiation  by  a  factor  of  ~30  -  an  order  of  magnitude  less  than  for 

c 

the  2p  cross-section  in  the  equivalent  Xe  L-shell  process  (Fig.  6(a)).  However,  the  initial  internal 
dephasing  rate  due  to  electron-electron  scattering  in  the  photoionized  2-electron  Kr  3s  orbital  is 
very  close  to  that  of  the  6-eIectron  Xe  4p  orbital  ionized  by  ATI,  since  the  initial  electron  densities 
are  nearly  identical.  Consequently,  after  inclusion  of  coUisional  dephasing,  the  ratio  between  the 
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two  cross-sections  is  again  expected  to  favor  the  shorter  248nm  pump-laser  wavelength  by  a 
factor  of  at  least  100. 

A  comparison  of  the  magnitudes  of  |A^(R  =  0)p  evaluated  for  the  generation  of  Xe(L)  and  Kr(K) 
emission  by  inner-shell  vacancy  production  (e.g.  Figs.  6(a)  and  8)  also  suggests  that  the  cross- 
section  for  Is  ejection  may  be  substantially  larger  than  that  for  the  impact  ionization  of  a  ^ 
electron  for  this  coherent  interaction.  This  is  predominantly  due  to  the  radial  dependence  of  the 
bound  atomic  wavefunctions  which  suppresses  the  dominant  contribution  to  the  matrix  element  of 
the  Coulomb  interaction  near  the  origin  of  the  coordinate  system  (Fig.  3)  when  1^0.  However, 
this  trend  may  be  reversed  in  a  full  evaluation  of  the  matrix  element  Mfi,  because  the  larger  2p 
bound-state  target  wavefunction  will  contribute  more  to  the  interaction  for  R  0  than  the  smaller 
Is  atomic  orbital.  Moreover,  the  2p  ejection  cross-section  will  be  enhanced  by  an  order  of 
magnitude  over  the  Is  ejection  cross-section  since  there  are  three  more  m-states  (and  electrons)  in 
a  full  p-orbital  than  in  a  s-state.  In  addition,  it  should  be  noted  that  the  relativistic  kinematic 
analysis  of  Ref  25  indicates  that  the  Lorentz  force  in  the  ponderomotive  potential  will  force  the 
driven  photoionized  wavefunction  (r2;/|v)  to  avoid  a  collision  with  the  target  atom(s)  in  gas  or 

atomic  cluster  materials  at  sufficiently  high  laser  irradiances.  Thus,  since  this  effect  is  much  less 
severe  for  shorter  radiation  periods  (due  to  the  reduced  impulse),  the  future  development  of  a 
coherent  A  x-ray  radiation  source  based  on  a  gaseous  target  material  strongly  favors  the  use  of 
high-power  short  wavelength  (UV)  rather  than  long  wavelength  (IR)  pump  lasers.  In  fact,  the 
combined  deleterious  effects  of  the  Lorentz  force  and  the  dephasing  of  (r2;/|\p)  by  electron- 

electron  scattering  may  be  the  major  limiting  factors  for  the  efficient  generation  of  the  inner-shell 
population  inversion  required  for  an  atom-based  coherent  hard  x-ray  radiation  source. 


V.  Conclusions 

The  selective  laser-driven  collisional  ionization  mechanism  outlined  in  this  paper  provides  an 
explanation  for  the  observation  of  spectrally  anomalous  and  unexpectedly  strong  L-sheQ  emission 
from  atomic  Xe  clusters  under  intense  irradiation  (10**-l0*^/cm2)  by  sub-picosecond  UV  TW 
laser  pulses  at  248nm  (Fig.  1).  In  particular,  it  elucidates  the  physics  behind  the  preferential 
generation  of  Xe  7p  vacancies  in  collisions  with  ponderomotively-driven  4p  electrons  that  are 
photoionized  by  the  strong  ATI  expected  in  the  observed  stable  channeled  propagation  of  the  UV 
laser  radiation  under  relativistic  self-focusing  conditions  [29,40,52].  Such  a  specific  coupling  of 
the  incident  laser  energy  into  target  atoms,  at  a  rate  that  may  exceed  IW/atom  [22,23,53],  could 
lead  to  the  inner-shell  population  inversion  required  for  the  development  of  friture  coherent  hard 
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x-ray  sources.  Indeed,  the  experimental  observation  of  -SkeV  L-shell  emission  from  both 
and  Xe28+(2/?53(^  already  provides  direct  evidence  that  laser-driven  inner-shell 
inverted  ions  can  be  efficiently  generated.  We  note  that  the  extension  of  this  mechanism  to 
uranium  would  generate  L-shell  emission  in  the  15-17keV  range.  Analysis  of  a  similar  selective 
collisional  ejection  process  to  produce  l5  vacancies  in  Kr  suggests  that  the  development  of  an 
efficient  K-shell  radiation  source  may  also  be  possible. 

The  analysis  of  this  novel  laser-driven  collisional  ionization  mechanism  also  indicates  an 

extraordinarily  strong  laser  wavelength  dependence  for  gas-phase-like  targets  (e.g.  small  atomic 

clusters  and  dilute  gases)  where  the  photoionized  electrons  may  return  to  interact  with  the 

atom/cluster  of  origin  after  approximately  one  radiative  period  [25].  Both  the  intrinsic  quantum 

mechanical  spreading  of  the  photoionized  wavefunction  and  its  dephasing  by  electron-electron 

scattering  ensure  that  the  efficiency  of  the  selective  collisional  ejection  process  is  reduced  by  a 

factor  of  -lO^  as  the  pump-laser  wavelength  is  increased  by  only  a  factor  of  ~3  from  the  UV 

(248nm)  to  the  IR  (800nm).  In  &ct,  for  sufficiently  long  laser  radiative  periods,  the  dephasing 

causes  the  progression  to  the  plane-wave  scattering  regime  for  which  the  relative  cross-sections 

indicate  that  the  generation  of  inner-shell  population  inversion  for  kilovolt  x-ray  emission  in  high 

Z  atoms  is  ruled  out.  The  analyas  also  shows  that  this  impossibility  is  robust  and  cannot  be 

overcome  by  the  use  of  higher  laser  powers,  such  as  petawatt  Ti:sapphire-  or  Nd-based  IR  laser 

systems  whose  long  radiation  wavelength  would  essentially  avoid  the  special  features  arising  from 

the  coherent  phase  dependence  of  the  interaction.  Form  this  point  of  view,  the  recent 

demonstration  of  an  efficient  third  harmonic  generation  scheme  for  high  power  IR  lasers  appears 

very  attractive  [54].  The  strong  pump-laser  wavelength  dependence  could  be  circumvented,  while 

retaining  the  /-state  and  significant  radial  (/i-state)  selectivity  of  the  collisional  ionization 

mechanism,  by  using  condensed-phase  targets  (e.g.  solid  targets  and  large  atomic  clusters) 

provided  that  the  ionization  coUiaons  occurred  within  a  quarter  of  radiation  period,  i.e.  the 

photoionized  electron  interacted  with  an  atom  other  than  its  atom  of  origin.  However,  the 

selectivity  of  ionization  process  is  likely  to  be  rapidly  impwed  by  scattering  with  ions  and  other 

electrons  in  the  laser-generated  plasma  as  the  photoionized  wavefunction  is  accelerated  to  the 

incident  collisional  energy  necessary  for  inner-shell  impact  ionization. 

« 

A  more  elaborate  and  sophisticated  calculation  of  the  matrix  element  Mfi  for  the  kinematics  of  the 
selective  collisional  ionization  mechanism  should  provide  further  insight  into  the  conditions 
required  for  the  optimally  efficient  generation  of  inner-shell  vacancies.  In  particular,  explicit 
evaluation  of  for  R  ^  0  should  allow  the  limits  of  the  strength  of  the  /-state  selectivity  of  the 
mechanism  to  be  determined.  A  comparison  of  the  cross-section  of  this  selective  mechanism  with 
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respect  to  the  plane-wave  case  is  also  required  to  quantify  experimentally  observed  enhancement 
of  the  former.  In  addition,  a  more  detailed  evaluation  of  the  effect  of  quantum  mechanical 
spreading  on  the  collisional  dephasing  dynamics  internal  to  the  photoionized  waveflinction  is 
necessary.  Nevertheless,  it  is  clear  that  the  selective  collisional  ionization  mechanism  provides  a 
new  mode  of  energy  coupling  to  atoms  which  could  result  in  the  efficient  and  controlled 
generation  of  inner-shell  population  inversion  required  of  potential  coherent  kilovolt  x-ray 
sources.  Furthermore,  the  evident  consistency  between  the  theoretical  analysis  and  the 
experimental  data  strongly  suggests  that  the  pump-laser  wavelength  is  a  substantially  more 
important  parameter  than  the  laser  power  or  intensity. 
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FIGURE  CAPTIONS 


Figure  1 

The  kilovolt  Xe  L-shell  x-ray  spectrum  emitted  by  5-20  atom  Xe  clusters  pumped  by 

TW  ultraviolet  (248nm)  and  infrared  (800nm)  sub-picosecond  laser  pulses.  The  charge-state 
designations  are  the  result  of  an  analysis  based  on  the  work  of  Cowan  [30],  which  also 
indicates  that  the  feature  at  2.6-2.6SA  is  attributable  to  the  generation  of  multiple  2p 
vacancies. 

Figure  2 

The  plane-wave  collision  cross-sections  o  for  the  ejection  of  3i,  3p,  and  “id  electrons  with 
ionization  potentials  of  8keV,  3.3keV,  2.8keV,  and  1.8keV  respectively  as  a  function  of  the 
energy  of  the  incident  colliding  electron  (after  McGuire  [32]).  For  convenience,  the  cross- 
sections  are  normalized  to  the  maximum  cross-section  for  the  ejection  of  a  2p  electron  at  a 
incident  collision  energy  of  ~30keV.  The  shaded  region  indicates  the  10-20keV  energy  of 
photoionized  electrons  returning  to  interact  with  their  cluster  of  origin  for  laser  irradiances  of 
10>8-10>9W/cm2. 

Figure  3 

The  interaction  geometry  for  the  ejection  of  an  electron  from  a  bound  state  {t\  via  a 
Coulombic  collision  with  a  ponderomotively-driven  photoionized  electronic  state  (r2;^|v|/) 
displaced  by  R  from  the  target  atom  (or  ion). 

Figure  4 

Comparison  of  the  exact  (solid  curves)  and  hydrogenic  (dashed  curves)  normalized  radial 
wavefunctions  rR{r)  for  the  2p,  3/>,  and  Ap  orbital  of  Xe.  The  hydrogenic  wavefimctions  are 
obtained  using  effective  nuclear  charges  of  SO,  42,  and  40  respectively. 

Figure  5 

The  calculated  value  of  |A/y;(R  =  0)|  as  a  function  of  the  normalized  expansion  time  tlx  and 
K  =  it,  cosP  for  the  coherent  coUisional  ejection  of  (a)  2p  and  (b)  3p  electrons  in  Xe  by  an 
expanding  photoionized  Ap  state. 
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Figure  6 

a)  The  value  of  [A/yj(R=  0)j  as  a  function  of  the  normalized  expansion  time  tlx  for  the 

collisional  ejection  of  a  Ip  electron  with  k  =  2lA’‘  (solid  line)  and  a  “ip  electron  with  k  = 
28A''  (dashed  line)  in  Xe  by  a  photoionized  Ap  state.  The  vertical  lines  indicate  the  periods 
for  248nm  UV  and  800nm  IR  radiation. 

b)  The  ratio  of  |A/y;(R  =  0)|  for  2p  versus  'ip  inner-shell  electron  ejection  as  a  function  of 
K-k^  cosP  for  a  collision  after  one  0.8fs  248nm  radiation  period  (solid  line)  and  one 
2.7fs  800nm  radiation  period  (dashed  line). 

Figure  7 

The  calculated  value  of  |A/y;(R  =  0)|  as  a  function  of  the  normalized  expansion  time  tlx  and 
K  =  k^  cosP  for  the  coherent  collisional  ejection  of  (a)  I5  and  (b)  25  electrons  in  Kr  by  an 
expanding  photoionized  35  state. 

Figure  8 

The  value  of  |A/y;(R  =  0)|  as  a  hiQction  of  the  normalized  expansion  time  tlx  for  the 

collisional  ejection  of  a  I5  electron  with  k  =  14A*‘  (solid  line)  and  a  25  electron  with  k  -  33A*' 
(dashed  line)  in  Kr  by  a  photoionized  35  state.  The  vertical  lines  indicate  the  periods  for 
248nm  UV  and  SOOnm  IR  radiation. 
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Abstract 

Evidence  of  highly  efficient  (1.2%  yield)  Ba(L)  (2.4  -  2.8A)  and  Gd(L)  (1.7  -  2.1  A)  x- 
ray  emission  produced  from  extensively  ionized  atoms  in  laser-excited  plasmas  in  BaF2  and  Gd 
solid  targets  is  presented.  An  analysis  describing  a  shell-selective,  coherent  collisional  ionization 
mechanism  is  given  in  order  to  explain  the  observed  highly  efficient  hollow  atom  production  in 
solids.  The  L-shell  spectra  were  observed  as  a  result  of  intense  irradiation  (peak  intensities  of 
10’®-10’®W/cm^)  generated  by  a  sub-picosecond  ultraviolet  (248nm)  terawatt  (lO’^W)  laser 
system  driven  by  a  hybrid  Ti:Sapphire/KrF*  laser  system.  A  mica-crystal  von  Hamos 
spectrograph  equipped  with  Kodak  RAR  2492  x-ray  film  and  appropriate  filters  was  used  to 
collect  the  spectra  emitted  by  the  targets.  We  discuss  the  limitations  on  producing  this  type  of 
radiation  imposed  by  the  laser  system  used  in  the  experiments  and  propose  a  method  for 
surpassing  this  limitation  with  little  or  no  additional  effort  or  expense  in  order  to  produce  L-shell 
radiation  on  the  order  of  15-17  keV  from  targets  as  heavy  as  Uranium. 
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I.  Introduction 

Recent  studies  have  reported  the  efficient  production  of  multi-kilovolt  x-ray  radiation 
resulting  from  inner-shell  transitions  in  rare  gas  clusters  [1-6]  as  a  result  of  intense  irradiation 
with  a  subpicosecond  ultraviolet  (248mn)  laser.  In  particular,  the  spectra  of  these  emissions 
were  highly  anomalous;  2p  vacancies  were  being  readily  produced  while  the  outer  subshells 
remained  largely  unperturbed.  This  property  seems  extremely  promising  for  achieving 
population  inversion  in  these  rare  gas  clusters  and  could  lead  to  the  coherent  amplified  emission 
of  these  highly  energetic  x-rays. 

In  an  attempt  to  explain  the  properties  of  this  anomalous  x-ray  generation,  Schroeder  et 
al.  [7,8]  proposed  a  new  collisional  ionization  process  in  which  the  ponderomotively-driven 
electrons  retained  a  substantial  portion  of  their  phase  and  symmetry  (i.e.  coherence)  from  their 
original  atomic  state  for  a  brief  period  of  time  after  ionization.  In  reference  [8],  the  authors 
found  it  necessary  to  consider  the  quantum  mechanical  wavefunctions  of  the  ionized  4/?  electrons 
when  examining  the  ionization  cross  sections  for  the  ejected  2p  electrons.  This  is  due  to  the  fact 
that  although  it  has  been  shown  that  a  conventional  collisional  ionization  process  may  produce 
large  ionization  rates  in  rare  gas  clusters  for  sufficiently  high  intensities  [20],  this  approach  fails 
to  explain  both  the  anomalous  nature  of  the  observed  spectrum  as  well  as  the  observed 
wavelength  scaling  of  the  production  of  the  Xe(L)  radiation  from  clusters  [6].  In  fact,  a 
calculation  of  the  plane  wave  collisional  cross  sections  for  the  specific  electron  orbitals  yields  a 
cross  section  for  the  ionization  of  a  3^  electron  that  is  an  order  of  magnitude  greater  than  that  for 
the  2p  electrons  [7,21].  This  does  not  agree  with  the  observations  published  in  [1,6],  which 
show  large  radiative  yields  from  electronic  states  with  2p  vacancies  where  the  3s  subshell  is 
largely  intact. 

In  reference  [7],  the  distinctive  quantum  mechanical  treatment  of  the  laser-driven 
electrons  resulted  in  the  collisional  ionization  dj-namics  being  altered  in  such  a  way  as  to  yield  a 
significant  w-state  and  /-state  selectivity.  In  the  case  of  Xe(L)  radiation,  the  assertion  was  made 
that  it  was  the  previously  ionized  4p  electrons  which  were  acting  to  generate  the  2p  vacancies. 
However,  the  theory  outlined  in  references  [7,8]  is  based  upon  the  supposition  that  the 
pondermotively  driven  photoionized  electron  subshell  is  able  to  propagate  unperturbed  by 
external  interactions  until  it  returns  to  the  cluster  or  ion  of  origin.  In  the  case  of  a  solid  material, 
this  assumption  is  not  valid  due  to  the  increased  density  of  atoms  in  the  vicinity  of  the  parent  ion. 
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In  this  paper,  we  extend  the  model  of  references  [7,8]  to  hollow  atom  L-shell  spectra 
from  2  solid  targets,  Ba(L)  and  Gd(L).  The  spectra  provide  evidence  of  this  shell-selective 
collisional  ionization  mechanism  in  solid  density  targets.  Two  Gd(L)  spectra  observed  under 
different  laser  intensities  will  serve  as  the  basis  to  establish  an  upper  bound  on  the  hollow  atom 
x-ray  generation  based  on  the  peak  laser  irradiance.  A  method  of  extending  this  limit  with  the 
current  laser  source  will  be  examined. 

TT.  Experimental  Setup 

The  laser  system  used  for  the  generation  of  the  spectra  presented  in  this  paper  is  a 
Ti:Sapphire/KrF*  hybrid  system  which  is  similar  to  that  outlined  in  reference  [9],  with  the 
exception  that  the  infrared  amplification  stage  has  been  replaced  by  a  Ti:  Sapphire  regenerative 
amplifier  [10],  The  system  produces  pulses  with  an  average  energy  of  400nJ  and  pulse  duration 
of  230fs.  An  f72  parabolic  mirror  was  used  to  focus  the  248nm  laser  pulse  onto  the  solid  target 
for  a  resultant  focal  intensity  of  ~10‘’W/cm^.  The  spectra  were  recorded  using  a  mica-crystal 
von  Hamos  spectrograph  in  3*^  order  and  Kodak  RAR  2492  film.  A  100pm  Kapton  (polyimide) 
filter  was  used  in  front  of  the  spectrometer  to  block  lower  energy  (M-shell)  photons.  The  target 
materials  used  were  a  BaFj  optical  flat  (to  generate  the  Ba(L)),  and  a  100pm  thick  Gd  foil 
mounted  on  an  A1  substrate  (to  generate  the  Gd(L)).  Since  A1(K)  radiation  is  below  2keV,  any 
radiation  from  the  A1  substrate  that  could  contaminate  the  measured  spectrum  was  removed  by 
the  Kapton  filter. 

nr.  Experimental  Data 

Figure  1  shows  the  recorded  Ba(L)  spectrum.  The  spectrum  was  produced  under  intense 
irradiation  by  a  pulse  energy  of  320mJ,  corresponding  to  an  average  laser  power  of  1.4TW  and  a 
focal  irradiance  of  ~10‘’W/cm^.  As  can  be  seen  from  the  figure,  the  broad  double-peaked 
structure  characteristic  of  the  3d^2p  transition  manifold  [11]  is  present,  with  charge  states 
ranging  up  to  38+  at  the  peak  of  the  spectrum,  as  labeled  in  Figure  1.  By  comparison  with  the 
Xe(L)  spectrum  described  in  reference  [1],  transitions  ranging  from  (2p'^3d^4-2p^3d^°)  to 
(2p^<-2p*3d)  are  present  in  the  spectrum.  In  addition,  the  Ba  spectrum  displays  the  prominent 
characteristic  transitions  Lai,  Lpi  and  Lp2  from  weakly  ionized  (0  -  10+)  Ba  atoms. 
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To  determine  the  efficiency  of  the  x-ray  yield,  a  calibrated  diamond  photoconductive 
semiconductor  device  (PCD)  was  used  to  record  the  yield  from  the  BaF2  target.  A  trace  from 
this  detector  is  shown  in  Figure  2.  The  temporal  width  of  the  signal  is  limited  by  the  response  of 
the  detector.  The  risetime  of  the  detector  signal  is  on  the  order  of  lOOps,  with  a  <200ps  decay, 
which  is  limited  by  the  response  of  the  detector  and  oscilloscope  (Tektronix  SCD5000).  A  filter 
consisting  of  0.  l|im  thick  A1  on  a  7.6|im  thick  Kapton  substrate  was  used  to  block  radiation  from 
Ba(M).  The  signal  showm  in  Figure  2  corresponds  to  a  total  energy  of  ~4.9mJ,  assuming  a 
uniform  In  distribution.  This  translates  into  a  conversion  efficiency  of  ~1.2%!  The  radiative 
lifetime  from  this  solid  target  is  in  stark  contradiction  to  the  observations  of  Ditmire,  et  al., 
which  showed  a  long  lived  emission  on  the  order  of  ~ Ins  [5,20]  from  Kr  clusters.  The  PCD 
signal  supports  the  conclusion  that  the  Xe(L)  emission  is  prompt;  indeed,  due  to  the  radiative  rate 
of  these  charge  state  configurations,  the  emission  could  be  as  short  as  the  lifetime  of  the  laser 
pulse.  However  further  experiments  with  greater  temporal  resolution  (i.e.  fast  x-ray  streak 
camera  measurements)  are  required  to  prove  this  assertion. 

Shown  in  Figure  3(a)  and  3(b)  are  two  Gd(L)  spectra  taken  under  identical  experimental 
conditions,  with  the  exception  of  the  peak  laser  intensity.  As  in  the  Ba(L)  spectrum.  Figure  3(a) 
shows  a  hollow  atom  spectrum  characterized  by  the  two  broad  spectral  features  with  charge 
states  ranging  up  to  40+,  in  addition  to  the  lines  Lai,  Lpi  and  Lp2  from  weakly  ionized  Gd  atoms. 
The  incident  laser  pulse  energy  was  350mJ,  corresponding  to  an  average  laser  power  of  1.5TW 
(irradiance  of  ~10*^W/cm^).  In  Figure  3(b),  the  hollow  atom  spectrum  is  absent,  leaving  only  the 
characteristic  lines.  Here,  the  laser  pulse  energy  was  200mJ,  a  factor  of  0.6  lower.  As  will  be 
discussed  in  the  next  section,  these  two  spectra  provide  crucial  evidence  to  support  the  theory 
presented  in  this  paper  and  in  reference  [8]. 

IV.  Solid  Target  Shell-Selective  Collisional  Ionization  Theory 

As  indicated  in  section  I,  the  theory  outlined  in  reference  [8]  makes  the  assumption  that 
the  pondermotively  driven  4p  electron  subshell  is  unperturbed  by  external  collisions  until  it 
returns  and  collides  with  the  parent  ion  or  cluster.  The  cross-section  for  this  interaction  is 
proportional  to  the  modulus  square  of  the  matrix  element  associated  with  the  interaction 
mechanism.  It  was  shown  in  reference  [8]  that  Mfi  may  be  expressed  as 
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(1) 


(R  =  0)  =  f f "--|^fcpexp[-  /k,  .(r,  -  )] , 

^1^2  **1  *"2 


where  <|)j(ri)  is  the  wavefunction  of  the  electrons  to  be  collisionally  ionized,  located  at  coordinate 
ri  around  the  atom,  vj/(r2,t)  is  the  wavefunction  of  the  ponderomotively  driven  electrons,  located 
at  coordinate  r2  about  the  center  of  the  photoionized  electron  cloud,  ki  is  the  wave  vector  of  the 
ejected  electron  after  the  collision,  and  R  is  the  displacement  vector  between  the  atom  and  the 
electron  cloud. 

In  reference  [8],  the  value  of  Mfi  was  evaluated  and  compared  for  UV  (248nm)  and  IR 
(800nm)  wavelengths,  based  on  the  assumption  that  the  collisional  ionization  occurred  after  one 
laser  period.  This  assumption  was  reasonable  because  of  two  points.  First,  it  was  demonstrated 
in  reference  [7]  that  it  was  kinematically  possible,  for  sufficiently  high  laser  intensities,  that  a 
previously  ionized  electron  may  collide  with  the  cluster  of  origin  with  more  than  enough  kinetic 
energy  to  ionize  a  Ip  electron  from  Xenon  (~8keV).  Second,  it  was  estimated  that  the  inter- 
cluster  separation  in  the  target  gas  jet  was  large  enough  that  the  electrons  could  be  accelerated 
back  into  the  parent  cluster  unperturbed  by  collisions  with  other  clusters. 

In  a  solid  density  material,  the  assertion  of  a  collisional  ionization  occurring  in  the  parent 
atom  is  not  generally  valid  because  the  inter-atomic  spacing  is  much  less  than  the  electron 
excursion.  Since  the  theory  governing  the  ionization  of  the  electron  wave  packet  is  unaffected 
by  the  target  density,  it  is  still  reasonable  to  assert  that  the  Ap  electrons  retain  their  previous 
quantum  mechanical  radial  and  angular  distribution  for  some  period  of  time  after  ionization.  The 
primary  difference  introduced  by  the  density  of  the  solid  target  is  that  in  this  case,  the  dephasing 
time  could  be  much  shorter,  due  to  the  fact  that  it  is  possible  for  the  electron  wave  packet  to  be 
dephased  by  electron-ion  collisions  in  addition  to  internal  electron-electron  scattering.  However, 
the  crucial  point  to  be  made  here  is  that  in  general,  there  is  no  reason  why  the  2p  vacancies  must 
be  created  in  the  same  ion  from  which  the  4p  electrons  originated.  That  is,  this  enhanced 
coupling  is  not  dependent  on  the  specific  ion  of  origin,  only  the  quantum  mechanical  nature  of 
the  two  electron  wave-packets.  Differences  between  ions,  such  as  actual  ionization  state  will 
have  a  small  effect  on  the  specific  bound  2p  wave  functions.  However,  this  will  only  change  the 
magnitude  of  by  a  modest  amount;  the  physics  governing  the  interaction  remains 

unchanged.  Compared  to  the  gas  cluster  case,  the  coupling  should  be  greater  in  solid  targets  due 
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to  the  fact  that  the  electron  orbital  will  have  had  less  time  to  dephase  due  to  internal  electron- 
electron  collisions.  Also,  the  overlap  of  the  waveflinctions  should  be  greater  because  the 
intrinsic  spreading  of  the  wave  packet  will  not  be  as  progressed. 

To  determine  if  this  shell-selective  collisional  ejection  theory  may  be  applied  to  solid 
targets,  the  only  issue  remaining  to  be  considered  is  whether  the  pondermotively  driven  electrons 
will  have  enough  kinetic  energy  to  eject  a  2p  electron  from  the  target  ion.  In  order  to  answer  this 
question,  tv^’o  pieces  of  information  are  required;  the  kinematics  of  the  pondermotively  driven 
electrons  (i.e.  the  kinetic  energy  as  a  function  of  position  or  time)  and  the  approximate  distance 
between  ions.  The  first  point  listed  above  has  been  calculated  in  reference  [7].  Specifically,  the 
authors  used  a  classical  relativistic  model  to  trace  out  both  the  position  and  kinetic  energy  of  the 
ponderomotively  driven  electrons  as  a  function  of  time.  Figure  4  shows  these  curves  as  a 
function  of  the  248nm  normalized  laser  period  ct/X.  The  calculations  in  Reference  [7]  were 
performed  for  Xe  clusters.  However,  the  numbers  are  sufficiently  similar  to  those  of  Ba  to  be 
applicable;  the  ionization  potentials  differ  only  by  ~5%  and  the  Xe  cluster  radius  rc  for  the 
calculations  was  5.3  A,  while  the  lattice  spacing  in  BaFj  is  ~6A.  For  simplicity,  all  calculations 
performed  here  will  be  done  for  Xe,  although  qualitatively  there  will  be  little  difference  in  the 
results  for  Ba  and  Gd. 

To  address  the  second  point,  one  finds  that  measuring  the  distances  between  ions  in  the 
solid  target  of  the  expanding  plasma,  may  be  extremely  complicated  [12],  In  order  to  quantify 
this  parameter,  detailed  knowledge  of  the  plasma  density  generated  by  the  laser  is  required.  In 
addition,  due  to  the  high  density  of  neighboring  ions  near  the  ion  of  origin  for  the  4p  electrons,  a 
rigorous  Monte-Carlo  simulation  must  be  done  to  approximate  the  average  frequency  of 
collisions.  While  these  details  are  necessary  in  order  to  precisely  estimate  the  magnitude  of  lMf,p 
for  each  solid  target,  they  are  not  required  to  ascertain  whether  the  theory  may  be  applied. 
Continuing  with  the  above  comparison  to  Xe  clusters,  we  approximate  the  distance  between  ions 
in  the  plasma  as  being  on  the  order  of  a  few  times  the  spacing  in  Xe  clusters. 

The  accuracy  of  this  approximation  can  be  determined  by  establishing  realistic  upper  and 
lower  limits  for  the  ion  separation  in  the  laser-driven  plasma.  Clearly,  a  reasonable  upper  limit 
for  the  interaction  distance  is  that  of  the  gaseous  target,  which  would  have  an  inter-cluster 
spacing  of -lOOA  [7].  It  has  already  been  shown  in  reference  [8]  that  the  theory  readily  applies 
to  this  case.  For  a  lower  limit,  a  first  approximation  would  be  the  actual  lattice  spacing  of  the 
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unperturbed  solid.  A  typical  number  for  this  spacing  is  in  the  range  of  SA.  As  can  be  seen  from 
Figures  4(a)  and  (b),  an  electron  will  reach  a  displacement  of  Ir/tc  (5.3A)  after  ~0.03  laser 
periods,  and  at  that  point,  will  not  have  the  required  kinetic  energy  to  ionize  a  2p  electron 
(8keV).  However,  from  examining  Figure  4,  it  is  seen  that  the  average  atom  separation  need 
only  increase  by  dL  factor  of  two  (r/tc  =  2  in  the  figure)  to  allow  the  electrons  sufficient  time  to 
accelerate  to  energies  large  enough  to  ionize  a  2p  electron.  Therefore,  using  an  ionic  spacing  on 
the  order  of  10-20 A  as  an  extremely  conservative  lower  bound  is  not  only  reasonable,  but 
sufficient  to  demonstrate  that  the  theory  may  be  applied  to  the  case  of  solid  targets.  This  lower 
bound  also  gives  an  estimate  of  the  delay  between  arrival  of  the  laser  pulse  and  emission  of  the 
L-shell  radiation,  as  well  as  a  time  scale  for  the  plasma  formation. 

In  addition,  the  argument  may  be  made  that  although  there  is  a  high  density  of  atoms  near 
the  parent  ion,  the  probability  of  the  ponderomotively  driven  electrons  interacting  with  a 
neighboring  ion  is  less  than  unity.  Therefore  it  is  entirely  possible  that  the  electrons  might  not 
interact  with  the  nearest  neighbor  ion,  but  the  next-nearest  neighbor.  Under  this  condition  of  a 
10  -  20 A  separation  in  the  plasma,  clearly  the  electrons  would  have  sufiicient  kinetic  energy  to 
eject  the  2p  electrons  at  the  time  of  collision. 

At  this  point,  a  qualitative  estimate  of  the  effects  on  |Mfif  for  a  given  solid  target  material 
may  be  made.  Figure  5(a)  and  (b)  show  a  plot  of  |Mf,p  for  barium  vs.  t/x,  where  x  is  the 
characteristic  time  for  the  radial  width  of  the  4p  wavefunction’s  initial  probability  distribution  to 
increase  to  V2  of  it’s  initial  value  after  ionization  [8].  As  in  reference  [8],  the  matrix  element  is  a 
function  of  the  parameter  k,  which  is  defined  as  k  =  -kicosp,  with  3  equal  to  the  angle  between 
the  wave  vector  ki  and  the  vector  ri  -  rj.  In  figure  5,  lower  values  of  k  are  used  in  generating  the 
curve,  (k  =  9.9A'^  for  the  2p  case,  and  16.2A'^  for  the  3p  case),  than  was  used  in  the  Xenon  case 
of  reference  [8]  to  reflect  the  fact  that  in  general,  the  electrons  may  have  lower  energy  when 
incident  on  the  target  ion  due  to  the  reduced  time  from  photoionization  to  the  first  collision. 
Qualitatively,  it  may  be  seen  in  figure  5a  that  as  t/x  decreases,  the  value  of  |Mf,p  will  increase, 
and  that  the  value  of  IMep  will  be  greater  for  an  ultraviolet  driving  laser  than  for  an  infrared 
driving  laser.  Also,  for  very  small  ion  separations,  the  actual  time  t  before  collision  for  the 
infrared  case  will  move  closer  to  that  for  the  ultraviolet  case.  In  the  limiting  case  of  the  lowest 
possible  time  before  collision,  the  laser  wavelength  dependence  of  Mfif  would  be  severely 
reduced,  since  the  collision  time  would  become  primarily  a  fimction  of  the  ion  spacing,  and  not 
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the  laser  period.  However,  at  this  time  there  are  no  data  for  the  infrared  case  in  these  solid 
targets  to  confirm  this  assertion. 


V.  Limitations  on  X-Rav  Generation 

Previously  in  section  in.B,  it  was  pointed  out  that  the  L-shell  spectra  for  gadolinium 
serve  to  establish  an  upper  bound  for  the  ability  to  generate  L-shell  hollow  atom  spectra  for  the 
laser  system  used  in  this  study.  Figure  3  (a)  and  (b)  show  two  Gd(L)  spectra  taken  at  laser 
intensities  which  differ  by  roughly  a  factor  of  1.6.  In  case  (a),  the  laser  intensity  is  sufficient  to 
produce  the  hollow  atom  spectrum,  while  in  case  (b),  it  is  conspicuously  absent.  As  is  shown  in 
Table  I,  a  peak  laser  intensity  of  approximately  1.4  x  10*^/cm^  is  required  to  ionize  the  entire 
4p  subshell  by  ATI  theory  [14-16],  as  is  given  by  equation  (3), 


lAii(w/cm')= 


4xlO^(lp(eV))^ 


(3) 


where  Ip  is  the  ionization  potential  of  the  electronic  state  to  be  ionized,  and  Z  is  the  charge  of  the 
residual  ion.  This  implies  that  the  peak  laser  intensity  must  be  greater  than  1.4xl0^^/cm^. 

An  estimate  of  the  laser  focal  intensity  can  be  made  from  the  laser  parameters  of  an 
average  power  of  1.5TW  and  a  measured  ~5|im  focal  diameter  (approximately  twice  diffraction 
limited).  These  numbers  yield  an  approximate  peak  focal  intensity  of  1.5  x  lO^^W/cm^,  which  is 
in  agreement  with  the  above  result.  Similarly,  the  peak  laser  intensity  for  figure  3b  of 
9xlO^*W/cm^  is  less  than  the  ATI  threshold  shown  in  Table  I. 

A  plot  may  be  constructed  showing  the  laser  intensities  required  to  produce  a  hollow- 
atom  L-shell  spectrum  as  a  function  of  atomic  number  (Z).  This  plot  is  shown  in  Figure  6.  The 
data  points  for  the  elements  are  plotted  according  to  the  laser  intensity  required  to  ionize  the 
respective  4p  sub  shells.  Indicated  on  the  lower  half  of  the  plot  is  the  intensity  limit  of  the  laser 
system  described  in  section  E.  According  to  the  plot,  hollow  atom  L-shell  radiation  can  be 
produced  for  heavy  atoms  up  to  Gd,  allowing  for  the  production  of  x-ray  photon  energies  of 
~6keV. 

Although  this  range  of  x-ray  production  can  be  extended  by  increasing  the  laser  intensity 
by  conventional  means,  it  is  more  advantageous  increase  the  laser  intensity  by  focusing  the  laser 
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into  an  underdense  plasma.  The  laser  beam  can  be  made  to  self-focus,  creating  a  channel  in 
which  the  laser  peak  intensity  is  on  the  order  of  lO^^W/cm^  [17],  These  channels  have  been 
readily  observed  in  gas  cluster  targets  [17,18],  Moreover,  recently  the  authors  of  reference  [19] 
used  a  pre-pulse  consisting  of  ~15%  of  the  total  laser  energy  to  create  an  underdense  plasma  in 
front  of  a  solid  target,  allowing  the  main  laser  pulse  to  channel  in  the  plasma.  As  a  result  of  the 
increased  laser  intensity,  the  authors  reported  observing  up  to  2.5MeV  photons  produced  by 
Bremsstrahlung  as  a  result  of  the  increased  kinetic  energy  of  the  ponderomotively  driven 
electrons.  Alternatively,  the  solid  target  could  be  positioned  behind  a  suitable  gas  cluster  target 
(Kr  or  Xe)  in  order  to  form  a  channel  before  the  solid  target.  Utilization  of  either  of  these 
techniques  could  increase  the  laser  intensity  by  2  orders  of  magnitude,  thus  allowing  for  the 
generation  of  hollow  atom  L-shell  spectra  as  energetic  as  15keV  from  a  Uranium  target,  as  is 
shown  in  Figure  6. 

VI.  Conclusion 

In  this  paper  we  have  demonstrated  the  generation  of  L-shell  radiation  with  energies  in 
excess  of  5keV  from  Ba  and  Gd  targets  with  efficiencies  on  the  order  of  1%.  The  mechanism 
attributed  to  this  high  efficiency  is  an  inner-shell  selective  collisional  ejection  process  in  which 
the  laser  energy  is  coupled  into  the  target  at  a  rate  in  excess  of  IW/atom.  The  two  Gd(L)  spectra 
shown  provide  direct  support  for  this  theory  of  collisional  ejection  when  correlated  with  the 
ionization  potentials  for  the  Gd  4p  electrons.  Finally,  this  paper  provides  a  means  of  predicting 
the  ability  to  produce  L-shell  spectra  from  highly  ionized  plasmas  as  a  function  of  laser  intensity, 
and  describes  a  simple  method  for  producing  hard  x-rays  from  elements  spanning  the  periodic 
table  up  to  and  beyond  15keV  U(L)  radiation. 
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Intensity  (a.u.) 


Figure  3b 


KE(cV) 


Material 

Ionization  Potential  Ip 
(eV) 

For  4p  electrons 

Charge  State 

Iati  (W  cm'^  ) 

Ba 

~  670  -  880 

20'"  -  25'' 

4  X  10** 

Gd 

-1100-1400 

28^-33'' 

1.4  X  10*® 

Ta 

-1700-2100 

ZT  -  42^ 

4.4  X  10*® 

Au 

-2300  -  2700 

+ 

00 

1 

+ 

9  X  10*® 

u 

-3700-4350 

56^ -6r 

3.7  X  10^° 

Table  I 


Figure  Captions 


Figure  1. 

High  charge  state  Ba(L)  plasma  spectrum  with  characteristic  lines  from  weakly  ionized 
(<  ~10"^Ba,  produced  by  intense  irradiation  (10^®  W/cm’)  of  a  BaF2  target.  The 
maximum  intensity  peak  results  from  charge  states  as  high  as  of 

Figure  2. 

Diamond  PCD  signal  of  the  L-shell  radiation  from  BaF2.  The  trace  shows  a  rise  time  of 
~150ps,  which  is  the  response  time  for  the  detector.  The  detector  has  a  calibration  of 
160pAAV,  which  corresponds  to  5.6mJA'^  for  the  trace  shown  here,  assuming  a  2:1 
radiation  distribution.  This  yields  a  total  energy  of  ~4.9mJ  of  energy. 

Figure  3 

(a.)  High  charge  state  Gd(L)  plasma  spectrum  with  characteristic  lines  from  weakly 
ionized  (<  -lO"^  Gd  ions,  produced  by  intense  irradiation  (10^^  W/cm^)  of  a  100pm 
thick  Gd  foil. 

(b.)  Gd(L)  characteristic  line  spectrum  only  with  a  background  due  to  Bremsstrahlung 
radiation,  produced  by  peak  laser  intensities  a  factor  of  0.6  less  than  figure  3(a). 

Figure  4. 

(a.)  Electron  displacement  vs.  normalized  laser  period  (cUX)  for  248nm.  Results  are  for 
the  fully  relativistic  numerical  simulation  of  the  kinematic  motion  of  an  electron  for 
an  intensity  of  I  =  lO^^W/cm*  under  the  influence  of  the  electric  and  magnetic  fields 


of  248nm  (UV)  radiation.  Data  is  taken  from  Reference  [7]  by  permission.  Plot  is 
in  units  of  radius  for  a  cluster  of  13  Xe  atoms  frc  »  5.3A).  The  horizontal  line 
indicates  rfrc  =  2,  and  the  vertical  lines  indicate  the  points  in  the  laser  period  where 
the  electron  reaches  displacements  of  rfrc  =  1  and  r/rc  =  2. 

(b.)  Electron  Kinetic  Energy  vs.  normalized  laser  period  (ct/X)  for  248nm.  Horizontal 
line  corresponds  to  8500eV,  the  approximate  ionization  potential  for  a  2p  electron 
in  Ba.  The  vertical  lines  in  both  plots  indicate  the  times  at  which  an  electron 
reaches  r/rc  =  1,  and  2. 

Figure  5. 

Plot  oflMfil^forBa. 

(a.)  The  long  dashed  vertical  lines  correspond  to  248nm  (t/x  =  19)  and  800nm  (t/x  =  63) 
laser  illumination  for  a  gas  cluster  target,  illustrating  the  values  of  for  the  case 
of  the  upper  limit  (~100A)  on  the  ion  spacing. 

(b.)  The  long  dashed  lines  correspond  to  248nm(t/x  =  1)  and  800nm(t/x  =  1.67)  laser 
illumination  for  the  case  where  the  average  ion  spacing  is  twice  the  lattice  spacing, 
illustrating  the  lower  limit  (~5A)  on  the  ion  separation.  As  can  be  seen,  the 
wavelength  dependence  would  be  much  less  severe  in  this  limit. 

Figure  6. 

Plot  of  required  laser  intensities  to  generate  hollow  atom  L-shell  spectrum  for  various 
elements.  The  second  horizontal  axis  shows  the  photon  energy  for  L-shell  radiation  from 
the  respective  elements. 


Table  I. 


Ionization  potentials  and  required  laser  intensities  to  ionize  the  Ap  electrons  by  above 
threshold  ionization  (ATI)  in  various  solid  elements. 


